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Some Roll Characteristics of Cruciform 
Delta Wings at Supersonic Speeds 


ZEGMUND O. BLEVISS* 
Douglas Aircraft Company, Ince. 


ABSTRACT 


Three roll problems associated with a missile-type configuration 
having a body of circular cross section, a cruciform delta wing, 
and some unspecified type of tail are considered. They are (1) 
roll due to aileron deflection of two opposite fins, (2) damping in 
roll, and (3) roll induced by pitch and yaw. Some effects of fin- 
fin, wing-body, and wing-body-tail interference are included. 

Some theoretical solutions are presented for the roll due to 
These results show clearly 
Theoretical 


aileron deflection and damping in roll. 
the effects of fin-fin and wing-body interference. 
calculations for the cross-flow field in the Trefftz-plane behind a 
cruciform wing with aileron deflection are presented for the two 
different assumptions of nondistorted and completely rolled up 
vortex sheets. 

A general qualitative discussion of these problems is presented. 
Effects of high angles of attack are discussed, along with some 
nonlinear, viscous, and gap effects 

Recommendations are made for future theoretical and experi- 


mental research. 


INTRODUCTION 


i PURPOSE OF THIS PAPER is to discuss qualita- 
tively three roll problems associated with a par- 
ticular supersonic missile-type configuration and to 
present some theoretical solutions to a few of these prob- 
lems. The configuration is shown in Fig. 1. It has a 
body of circular cross section, a cruciform delta wing on 
a portion of the body which is practically cylindrical, 
and a tail that is not specified as to type or shape. 

The basic purpose of a cruciform wing is to be able 
to supply lift equally in any desired lateral direction by 
some combination of deflections of the fins. Although 
any wing composed of nonplanar lifting surfaces has 
this property, calculations are simplest for the cruciform 
configuration. The choice of the delta plan form and 
the circular body cross section was dictated by both the 


Presented at the Session on Stability, Control], and Guidance, 
Annual Summer Meeting, I.A.S., Los Angeles, July 12-14, 1950. 

* Aerodynamics Consultant, Santa Monica Plant. Also, Grad- 
uate Student, California Institute of Technology 


simplification of the theoretical problems and, more 
practically, their actual use in some supersonic mis- 
siles. 

The three roll problems considered are (1) roll due to 
aileron deflection of two opposite fins, (2) damping in 
roll, and (3) roll induced by pitch and yaw. Within 
each of these problems it is convenient to classify a sub- 
group as interference problems. The subgroup con- 
sists of fin-fin, wing-body, and wing-body-tail inter- 
ference. Fin-fin and wing-body interference are de- 
fined as the change in any characteristic of a plane wing 
due to the addition of a fin or body, respectively. 
Wing-body-tail interference will be used here to mean 
the change in any characteristic of a wing-body com- 
bination due to the addition of a tail. 

Most of the theoretical results are based on linearized 
nonviscous theory, and, where this is not the case, the 
assumptions are explicitly stated. Furthermore, when- 
ever fins are deflected, the effect of the gap is neglected. 


NOTATION 


It is convenient for illustrative purposes to consider a cruciform 
wing to be oriented so that one pair of fins is horizontal and the 
other pair vertical. 

Wide Delta.—All leading edges supersonic. 

Narrow Delta.—All leading edges subsonic. 

Plane Wing.—Flat plate of zero thickness with delta plan form, 
unless otherwise stated. 

Plane Aileron.—Plane wing with 
the two halves of the plane wing deflected in opposite directions 


aileron deflection—i.e., 
by +6. 

Cruciform Wing. 
norinal to each other and symmetrically arranged about the line of 


Two plane wings with the same plan form 


intersection. 


Cruciform Aileron.—Cruciform wing with aileron deflection of 
The general features are illustrated in Fig. 2. 


Cruciform aileron with wide delta 


the horizontal fins. 
Mixed Cruciform Aileron. 
vertical fins and narrow delta horizontal fins. 


A = body diameter/wing span = 2r/b 


b total span of horizontal fins 
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FIG/ 
MISSILE TYPE CONFIGURATION 
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FIG.2 
CRUCIFORM AILERON 
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FIG.3 
ROLL DUE TO AILERON DEFLECTION AND DAMPING IN ROLL 
FOR WIDE DELTA CRUCIFORM WING 
AILERON DEFLECTION DAMPING IN ROLL 
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Cc = chord of wing extended into the body 
¢’ = chord of elementary strip (Fig. 5) 
Ci = rolling moment coefficient about body axis 


based on area and span of horizontal fins 

extended into the body; the positive roll 

direction is indicated in Fig. 2. 
Ci/(pb/2U) = damping in roll coefficient 


‘ot 
N 
8 

ll 
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d = tan w/tan u (for subsonic leading edges) 
dé = differential width of elementary strip 
dL = differential lift of elementary strip 
dM = differential moment about the axis of the 


elementary strip; positive when counter- 
clockwise facing upstream 
complete elliptic integral of the second kind 


E(V/1—d) = 
with the modulus 1 — d? 


f = tan w/tan w = 1/d (for supersonic leading 
edges) 

M = Mach Number 

m = 1/V M?-1 

p = angular rolling velocity 

q = (1/2)pU? = free-stream dynamic pressure 

r = body radius or radial distance from origin 
in Trefftz-plane (Fig. 10) 

U = free-stream velocity 

uj (7 = 1, 2,3) = perturbation velocities 

n(6) = cross-flow velocity normal to the radius 
from the origin in the Trefftz-plane 
(Fig. 10) 

xj (J = 1, 2, 3) = Cartesian coordinate system (Fig. 2) 

a = angle of attack 

B =1/m = VM?-1 

6 = aileron deflection angle 

0 = angle from the x, axis of the radius from the 
origin in the Trefftz-plane (Fig. 10) 

im = Mach angle 

E = distance between the elementary strip and 
the wing-body juncture 

&° = §/mc’ 

p = free-stream density 

w = half vertex angle of the delta plan form 


THEORETICAL RESULTS 


Cruciform Wing—No Body 

If the leading edges are supersonic, exact solutions 
are readily obtained for roll due to aileron deflection 
and damping in roll. The results are shown in Fig. 3. 


Consider first the roll due to aileron deflection. The 
configuration is conical, and the solution is obtained 
from known conical wing solutions. Since the vertical 
fins have supersonic leading edges, the horizontal fins 
are independent of each other. The presence of the 
vertical fins requires that the sidewash velocities be 
zero in the plane of these fins. Hence, by symmetry, 
each horizontal fin acts like half of a wide delta plane 
wing with the angle of attack of the fin. The results in 
Fig. 3 are replotted from reference 1. Results for the 
plane aileron are plotted for comparison to illustrate the 
effect of fin-fin interference. The rolling moment co- 
efficient for the plane aileron is independent of f and 1s 


given by 


6C./3 = 2/3 (1) 
When the vertical fins are added, the roll carried by the 
horizontal fins is considerably increased above that for 
the plane aileron. However, unfortunately, a counter- 
rolling moment is induced on the vertical fins which 
more than offsets this gain, and the final result is a de- 
aileron effectiveness. The rolling moment 
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CHARACTERISTICS © 


ROLL 


coefficient for the complete cruciform aileron is given by 


: 4 1 — 2f?) 
ss = - >|: f +5] (2) 
6 3x (1 — f?) 


VI-f? 
[he result for four deflected fins is obtained from this 
by superposition. The reduction in rolling moment 
due to fin-fin interference is a maximum of 15 per cent 
when the leading edges are sonic. 

The results for damping in roll are similar. 
tion to this problem is again reduced to the solution of 
a plane wide delta wing, but, in this case, the wing is not 
flat. With the usual assumption that the effect of roll- 
ing can be represented by an apparent linear angle of 
attack distribution, the problem is reduced to that of a 
Since this 


cos* f 


The solu- 


plane wing with linear symmetrical twist. 
wing is not conical, more general wing methods must 
be used. Here, the comparison must be made with two 
plane wings for which the damping in roll is independent 
of f and given by 


BC, = —? oe (3) 
If these two plane wings are combined to form a cruci- 
form wing, the effect is generally the same as for the 
This case cannot be separated conveniently 
the corre- 


aileron. 
into horizontal and vertical fins. Instead, 
sponding rolling moments are called direct damping (the 
damping moment carried by the fin itself) and fin-fin 
induced roll (the counter-roll induced by fin inter- 
Here, again, the result is a decrease in the de- 
The 


action). 
sired characteristic, in this case damping in roll. 


equation for the cruciform wing damping in roll is 


if 
| (4) 


V1 —f 


» 


pC, = _ faye se — df?) + 
—_ ad ho 


7 3r(1 


Cos 


(2 — df? + G6f*) 


Fin-fin interference has less effect in this case. It 
causes a maximum reduction of damping in roll of about 
7 per cent for sonic leading edges. 

When the leading edges are subsonic, the problem is 
much more difficult, and, as yet, no exact solutions are 
known. In this paper an attempt is made to approxi- 
mate the solution for the cruciform aileron in the region 
where the leading edges are not too subsonic (say, 0.75 
<d< 1.0). Actually, the method is nothing more 
than a reasonable guess at the solution from a knowledge 
of some limiting cases combined with under- and over- 
{, the known exact solution for the 


estimates. In Fig. 


plane aileron'is shown. The equation for this curve is 


6C,/é = (2/3)d (5) 


Take Then 
let the vertical fins increase in size from the limiting case 
of no fins (plane aileron) to the other limiting case of 
fins with sonic leading edges (mixed cruciform aileron). 
It seems reasonable that the variation in rolling moment 
during this process will be continuous and monotonic 


a fixed value of d for the horizontal fins. 


F 
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FIG.4 
APPROXIMATION TO ROLL DUE TO AILERON 
DEFLECTION FOR NARROW DELTA 
CRUCIFORM WING 
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FIG.5 
ESTIMATES USING ELEMENTARY STRIP 
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and that the value of rolling moment for the cruciform 
aileron will lie somewhere between the two limits. In 
Fig. 4, the solution for the mixed cruciform aileron is 
shown. This solution is readily obtained from that for 
the plane narrow delta wing’? using symmetry argu- 
ments. The end points atd = Oandd = 1.0 are known 
for the cruciform aileron, and it appears that a straight 
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line connecting these end points should give a reason- 
able approximation at least for the range of d values 


mentioned above. The equation of this line is 


BC,/5 = (16/9m) d (6) 
and it is shown in Fig. 4 as the cruciform aileron ap- 
proximation. When the slender body solution is avail- 
able, the slope will be known at d = 0 and a better ap- 
proximation can be faired in. From these consider- 
ations it appears that fin-fin interference has a detri- 
mental effect for all leading-edge conditions. It is 
believed that this will also be true for damping in roll. 


Cruciform Wing—with Body 


No attempt is made to carry out any exact solutions 
when a body is present. Instead, a simple procedure 
for estimating the wing-body interference is used. 
Using this estimating procedure, results are obtained 
for the three cases treated above. 

It is most convenient for these estimates to use the 
elementary strip’ * 4 (see Fig. 5). The pressures pro- 
duced by this strip may be obtained by interpreting 
the strip as a line of sources of constant strength deter- 
mined by the wing downwash prescribed along that 
line. It is more useful in these applications to inter- 
pret the strip as a lifting wing of differential width dé 
whose length dimension is in the free-stream direction 
and whose angle of attack is constant. Adjacent to the 
wing on both sides is a flat plate at zero angle of attack. 
On one side the plate is always ahead of any disturb- 
ances, while on the other side the leading edge is 
straight and may be subsonic or supersonic depending 
on the wing to which the calculations are to be applied. 
Such a strip is intended to be applied, of course, only to 
wings that are not cambered. 

Now consider this strip as a deflected element in one 
fin of a delta wing-body combination as shown in Fig. 
5. When the strip is close enough to the body, disturb- 
ances will strike the body at the wing-body juncture. 
Part of the disturbances will be diffracted and pass over 
the body downstream of a helical path, and the rest will 
be reflected back onto the fin. The part that passes 
onto the body can contribute nothing to the rolling 
moment about the body axis, since the body cross sec- 
tion is circular. It is this fact that makes a simple esti- 
mation of the wing-body interference possible. The 
disturbances that are reflected back onto the fin increase 


Seri 


(1 + 


ans tmqgc'(&) a(S) J- iesaiiaod | a 
| 2d + EI 


T 


2m*qgc'?(E) a(é) , = 
ae ) va Vd+(1-di-F-FP 


T 
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the lift of the fin and, hence, increase the rolling moment 
carried by the fin. 

It is now apparent how the contribution of a fin to the 
rolling moment may be bracketed by under- and over- 
This is illustrated in Fig. 
supersonic leading edges. It is also 
the leading edges are subsonic. An 
obtained by assuming that no reflection takes place at 
With the elementary strip 


estimates. 5 for the case of 
true, however, if 


underestimate is 


the wing-body juncture. 
this amounts to using only that portion of the pressure 
distribution which is cross-hatched in Fig. 5a. 
stimate is obtained by assuming that complete re- 
flection takes place. This is realized if the body is re- 
placed at the wing-body juncture by a reflecting wall 
that has supersonic leading edges. This is illustrated 
in Fig. 5b, where the region with double cross-hatching 
indicates the added pressures due to reflection. 

supersonic 


An over- 


Using well-known methods of wing 
theory, formulas for the differential lift dL and the dif- 
ferential rolling moment about the axis of the strip 
dM are easily derived. When the disturbances from 
the strip do not strike the body before passing off the 
fin, the solutions for the strip are exact and are given by 


dL = 4mqce'(é) a(&)dé (7a) 
dM =0 (7b) 
for a supersonic leading edge and by 
dL = 4mq Vd c'(&) a(&)dé (Sa) 
dM = —m’q Vd (1 — d) c'*(é) a(€)dé (Sb 


for a subsonic leading edge. When the disturbances 
strike the body the under- and overestimates are given 


by 
4mgc’ (&) a(é) . ‘ ‘ 
dL = = ~~ (r + Esech~! — — cos! £) dé (9a) 
T 
2m°qc'*(E) a(&) = - a 
dM “(V1 — F — EF sech-'£) dé (9b) 
Tv 
and 
dL = 4mqgce'(é) a(&)dé (10a 
tm?gce’*(E) a(€) / — 
dM = (V1-—-2 + 
T 
E? sech—! E— 2E cos—! E) dé (10b 
respectively, if the leading edge is supersonic. If the 


leading edge is subsonic, a single application of Ev- 
vard’s method yields the underestimate 


d) |+ 5 coe | C=O — EN i 
~<a V ad COS rr (4s (lla 
: | E(1 + d) |- 
_— 2d + E(1 — d) 
‘ad(1—d 1 — d) — 2 
Vd ( d) oe af d) has (1 1b) 
Z I l+d f 
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ROLL 


For the overestimate it is necessary to solve the prob- 
lem for a strip in the presence of two subsonic edges. 
The solution requires repeated use of Evvard’s method 
and would lead to rather complicated expressions 
for dL and dM. This solution has not been carried 
out. 

These estimates may be applied only if the disturb- 
ances propagating around the body do not encounter 
other fins; for, if this occurs, rolling moments will be 
induced which the estimates do not account for. Then, 
to apply the estimates, the body must be large enough 
so that the fins do not interact. Obviously, the limiting 
cases where the disturbance from a fin leading edge-body 
juncture just strikes the trailing edge of the next fin will 
be different for the plane and cruciform wings. This 
is illustrated in Fig. 6. For supersonic leading edges, 
the estimates may be used for a cruciform wing if 


A > 2f/(e + 2f) 
and, for a plane wing, if 
A > f/(r + f) 


The corresponding regions of validity for subsonic lead- 
ing edges are 


A > 2/(2 + rd) 
for a cruciform wing and 
A > 1/(1 + wd) 


for a plane wing. 

Now strips can be superposed by integration to give 
the under- and overestimates for any twisted wing. 
If the spanwise distribution of angle of attack is com- 
plicated, the integrations will have to be carried out 
graphically or numerically. 

The method will first be applied to the problem of roll 
due to aileron deflection when the leading edges are 
supersonic. The angle of attack is constant across the 
span, and, by use of Eqs. (9), the underestimate is found 
to be 


BC; = (1 5 {20 4 9A) cos 7! f 
6 om V1—-f? 


w(1 + 2A) + 2f(1 — A | (12) 


In like manner, using Eqs. (10), the overestimate is 

BC 2(1 — A)? cos! f 

i ; Ez + 21 — A)( — 4 ‘) 
0 oT ea f 


The results for wing-body combinations are given in 
Fig. 7 for both the cruciform and plane ailerons for three 
values of the parameter f. For a given value of f, the 
curves for the under- and overestimates are indicated 
by the same type of line. Each pair of curves is 
broken at a value of A indicated by a short solid line. 
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FIG.6 
LIMITS OF FIN-FIN INDUCED EFFECTS 



































(b) PLANE WING 








FIG7 
ROLL DUE TO AILERON DEFLECTION 
FOR WIDE DELTA WING-BODY COMBINATIONS 
CRUGIFORM AILERON PLANE AILERON 
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These are the values of A above which the estimates are 
valid and, for the sake of brevity, will be referred to as 
A,. ForA > A, the curves are given by Eqs. (12) and 
(13). The end points for A = 0 are known exactly from 
Fig. 3. The curves are then completed by drawing in 
curves from the known values at A = 0 to the known 
estimates at A = A,;. These curves are not drawn in 
tangent to the estimate curves, since it is expected that 
the exact solution to this problem will consist of two 
different expressions exhibiting the same type of slope 
discontinuity. Furthermore, the curves are drawn 
convex upward, since this appears to be the most rea- 
sonable shape. In most cases it is rather easy to com- 
plete the curves in a reasonable fashion. 
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The difference between the value for the plane 
aileron at A = 0 and any other value is an indication 
of the effects of fin-fin and wing-body interference. At 
A = 0, any decrease in rolling moment is due entirely 


Damping in roll for supersonic leading edges is treated in the same manner. 
Integrating out, using Eqs. (9) and (10) 


tributed linearly across the span. 


, = = ar * ST — 3f?) (1 2A) 3A2(2 
BCi, -_ mi ~~ =~ le if ( =} ~ + (2 


and the overestimate 


9/ — 2 
ro am {0 — A) [(2 — 3f%) + A(2 
3x(1 — f*) 


BC, - 


The results are shown in Fig. 8 with the exact end points at A 
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to fin-fin interference. 


ference. 
body interference. 


, yields the underestimate 





~ 1 
cos Ss 4 el — pf) (1+ 24+ 3A + 
om: 

f(2 — f?) (1 — A) (1 + 3A) (14) 





— f*) (1 + 2A) + 


A)[(2 — f?) + AQ - apt (15) 


= 0 taken from Fig. 3 


The cruciform aileron with subsonic leading edges can be treated even though the elementary strip overestimate 


is not known. 


BC, a 2(1 i A)? j 


= ——— _¢ [3d(1 + A) 
6 or ] 


— (1 — A)]- 


Using Eqs. (11), the underestimate is given by 


Vd . 3d(1 + A) + ( 
Vd l+d 


tan~! 


~ A 
2 (16) 


The overestimate can be given if the solution for a fin with a subsonic leading edge adjacent to a wall with super- 


sonic leading edges can be obtained. 
overestimate is 


BC, 
= [(1 
; [( 


The results are given in Fig. 9. The end points at A 
those for the cruciform aileron are approximate. 


Cross-Flow in the Trefftz-Plane 


Whenever any of the fins of the cruciform delta wing 
carries lift, a vortex sheet passes downstream from the 
trailing edge of that fin. 
with all fins lifting will have four vortex sheets trailing 
These sheets will de- 


A wing-body combination 


downstream about the body. 





FIG.8 
DAMPING IN ROLL FOR WIDE DELTA WING-BODY COMBINATIONS 
CRUCIFORM WING Two PLANE WINGS 
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This is just the mixed cruciform aileron solution described above. The 


] [4 + A(3x — 4)] (17) 


= (0 are obtained from Fig. 4, where, it will be remembered, 


flect and distort under the influence of the body and the 
spatial distribution of vorticity. This is an extremely 
nonlinear problem whose exact solution is beyond pres- 
ent-day mathematics. Associated with the vortex 
sheets is a cross-flow field that will induce lift on any 


tail fins. This is what is referred to as wing-body-tail 
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ROLL DUE TO AILERON DEFLECTION 
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For 0 < A < A,, the reduc. 
tion in rolling moment is due to both types of inter. 
For A > A), the effects are due solely to wing. 


Here, the angle of attack is dis- 
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ROLL 


interference. These induced tail lifts and rolling 
moments act generally to oppose those of the wing. 

It is obviously important to be able to estimate the 
cross-flow field at the tail. Until recently, most theo- 
retical work on the cross-flow field behind plane wings 
in supersonic flow assumed that the vortex sheet re- 
mained undistorted* ® (see reference 6 for further refer- 
ences), although, in some cases, attempts were made to 
estimate the deflection.‘ It was found that for undis- 
torted vortex sheets the Trefftz-plane (infinitely for 
downstream) values of cross-flow velocities were ob- 
tained for all practical purposes a few chord lengths 
Some elegant Trefftz-plane 
papers*® 


downstream of the wing.® 
methods were developed.*’ 
indicate that, for wings of low aspect ratio such as con- 
sidered in this paper, the trailing vortex sheets are cer- 
tainly rolled up by the time the tail position is reached. 
In reference +, some study is made of a plane wing-tail 
interference problem for the two assumptions of a de- 
flected but nondistorted vortex sheet and the vortex 
sheet replaced by two deflected incompressible vor- 
tices. In reference 9, the effect of a circular body on the 
motion of two incompressible vortices is studied. 


Two. recent 


Little theoretical work has been done on the cross- 
flow field behind wings. However, the 
methods and results for the plane wing problems, which 
are summarized above, act as useful guides for the study 
of the cruciform problems. At first, Trefftz-plane 
solutions were sought. It was soon found that the ele- 
gant Trefftz-plane methods used with plane wings gen- 
erally do not work. An approximate Trefftz-plane 
solution for a wide delta cruciform aileron is given in 
The vortex sheets are assumed to be un- 
The cross-flow field gen- 


cruciform 


reference 1. 
deflected and undistorted. 
erated by the vortex sheets from the horizontal fins 
is replotted in Fig. 10. For comparison, the assump- 
tion that the horizontal vortex sheet is rolled up into 
two undeflected incompressible vortices located at the 
centroids of vorticity from each fin is shown. This com- 
parison is made only for the horizontal fins, because the 
vorticity distribution at the trailing edge of the vertical 
fins is such that the vortex sheet cannot be assumed to 
roll up into two distinct vortices. In Fig. 10, the 
vortex sheet and the two vortices lie along the x; axis, 
the vortices being located at x, = +0.96 mc. The 
velocity components normal to a radial line through the 
origin are shown. These results are presented to illus- 
trate the type of cross-flow calculations that can be 


made for the cruciform wing. 


QUALITATIVE DISCUSSION 


The above theoretical results are based on linearized 
nonviscous theory. Even in the results presented for 
the cross-flow field at the tail where linear theory was 
not used downstream of the wing, linearized nonviscous 
theory was used to determine the needed vorticity distri- 
bution at the wing trailing edges. In actual use, the 
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FIG. 10 
CROSS FLOW IN TREFFTZ-PLANE 
DUE TO HORIZONTAL FINS OF 
WIDE DELTA CRUCIFORM AILERON 
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cruciform wing may be subjected to angles of pitch and 
yaw far above the ranges within which linearized theory 
can be expected to give reasonable results. In these high 
angle-of-attack ranges, the nonlinear and viscous effects 
and the effect of the gap between a deflected fin and the 
These effects are much more 
The follow- 


body will be considerable. 
severe behind the wing than on the wing. 
ing qualitative discussion attempts to indicate the 
difficulties associated with the various roll problems at 
large angles of pitch and yaw and, where possible, the 
qualitative effects on the roll characteristics. 


Aileron Deflection 


Consider first the case of aileron deflection with no 
pitch or yaw. The angles through which the fins are 
deflected are normally within the range of linearized 
theory applications. A certain loss in lift must be 
expected because of the gap between the deflected 
fin and the body through which some pressure equaliza- 
tion takes place. Furthermore, there may be other 
lift losses due to viscous effects such as separation. In 
general, this leads to rolling moments somewhat less 
than predicted by theory. However, it is expected that 
the above theoretical results furnish a good approxi- 
mation to the effects of fin-fin and wing-body inter- 
ference. 

Now if this aileron is given an angle of attack, the 
angle of attack of one fin is decreased, while the angle 
of attack of the other fin is increased. The fin at the 
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high angle of attack carries proportionately more of the 
lift that generates the rolling moment. At the same 
time, the gap and separation losses become more 
severe. The result is probably a decrease of aileron 
effectiveness with angle of attack. 

If, instead of pitching, the aileron is yawed, some roll- 
ing moments must be expected because of lateral asym- 
metry. The body and viscous effects will certainly be 
important factors, and, hence, it is expected that theo- 
retical results for yawed delta wings usually will not be 
applicable. From physical reasoning it does not seem 
possibie to even predict the sign of these additional roll- 
ing moments. 

All the considerations up to this point have dealt with 
the aileron as a stationary flow problem. Actually, of 
course, when the fins are deflected in flight, there is a 
time interval before stationary flow is established. 
Knowledge of the rolling moments in this transient 
flow -region is important in stability and control 
considerations. No theoretical results are avail- 
able. 

If there is a tail on the body, the roll induced by the 
vorticity from the aileron is usually opposite in sign 
to that of the aileron. This effect is generally much 
greater than the effects of fin-fin or wing-body inter- 
ference and results in reduced aileron effectiveness. 
When the body is pitched or yawed or both, the tail 
moves into flow regions where the induced cross-flow 
velocities are lower. This decreases the induced tail 
counter-rolling moments. 
interference effects decrease as the angle between the 
free-stream direction and the body axis is increased. 


In short, the wing-body-tail 


Damping in Roll 


The theoretical results for damping in roll for zero 
pitch and yaw were obtained by replacing a nonsta- 
tionary problem by what is assumed to be an equivalent 
It was assumed that rolling in- 
duced an apparent linear twist. The twisted cruciform 
wing was then solved as a stationary problem. The 
magnitude of the nonstationary effects is not known for 
the cruciform wing, but results for plane wings indicate 
that these effects can be large. 

In any case, gap and separation effects will probably 


stationary problem. 


decrease the actual damping in roll below that predicted 
by the correct linearized solution. 

Another nonstationary problem arises for a rolling 
wing at an angle of attack. The flow will have lateral 
symmetry when the fins are horizontal and vertical or 
when both sets of fins are at 45°. However, for all other 
orientations, the asymmetry of the flow over the lifting 
surfaces will induce rolling moments that vary with the 


orientation. It is not known how important these 
effects are, and no theoretical results are avail- 
able. 


If the body has a tail, the same remarks apply as for 
the cruciform aileron case. 
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Pitch and Yaw Effects 


Consider a cruciform wing that is in combined 
pitch and yaw. If the angles of pitch and yaw are 
equal or if either angle is zero, no rolling moments can 
be induced. For all other orientations, rolling mo- 
ments must be expected because of the asymmetry in 
the flow. Linearized theory cannot predict the roll- 
ing moments,' ! since these moments must be at least 
of second order—i.e., proportional to the product of 
the angle of attack and the angle of yaw. It appears 
that this problem will have to be attacked by experi- 
ment. 

Another effect of pitch or yaw which can be im- 
portant is the following: Consider a cruciform wing 
that has the horizontal fins—say, at incidence. When 
the wing is yawed, lateral asymmetry will induce roll- 
ing moments. 

The tail effects are qualitatively the same as for 
aileron deflection and damping in roll. 


Wing-Body-Tail Interference 


The difficulties associated with the computation of 
wing-body-tail interference are discussed here. 

The first thing that must be known is the vorticity 
produced by the fins. For aileron deflection and damp- 
ing in roll, the vorticity distributions at the trailing 
edges of the fins can be computed by linearized theory 
if the leading edges are supersonic and there is no 
With a body present, an exact linearized solu- 
Any computa- 


body. 
tion would be difficult to compute. 
tions for combined pitch and yaw appear to be out of 
the question. At high angles of pitch and yaw, the 
actual vorticity distributions will undoubtedly differ 
considerably from any linearized computations because 
of the nonlinear, viscous, and gap effects. 

As remarked above, it is now known that the vortex 
sheets are rolled up a short distance downstream of the 
wing. It appears that, if incompressible vortices are 
assumed to start at the trailing edges of the wing, use- 
ful engineering results can be obtained. A vortex is 
placed at the centroid of the vorticity distribution at 
the trailing edge of each fin. The strength of the 
vortex is put equal to the total circulation at the fin 
If the vorticity is all of one sign, a 
However, 


trailing edge. 
single vortex is assumed to trail from the fin. 
there are cases for which the vorticity at the trailing 
This occurs, for example, for a 
For these 


edge changes sign. 
plane aileron or a cruciform wing in roll. 
cases the vortex sheet probably rolls up into two dis- 
tinct vortices of opposite sign and different strengths, 
and the vortex pattern behind a cruciform configuration 
might consist of two vortices trailing from each fin. 
Next, the motion of the vortices as they proceed 
downstream to the tail must be studied. For the cases 
of aileron deflection and damping in roll with the body 
axis in the free-stream direction, the motion of the vor- 
tices as they proceed downstream will be roughly in a 
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direction tangential to the body with some radial dis- 
placement if only a single vortex is assumed to trail from 
Results for these cases can probably be ob- 
If two vortices trail 


each fin. 
tained without excessive labor. 
from each fin, the labor of calculation is greatly in- 
creased. Certain phases of the vortex motion behind a 
cruciform wing-body combination in combined pitch 
and yaw will definitely have to be studied experimen- 
tally. 

Finally, once the position and strength of the vor- 
tices are known at the tail, the induced angle-of-attack 
distributions on the tail fins can be computed. If the 
tail is a wide delta cruciform configuration, the ele- 
mentary strip estimate method may be used to deter- 
mine the induced rolling moment. The integrations 
would be carried out graphically or numerically. For 
other types of tails, estimating procedures would have 
to be developed. 


RECOMMENDATIONS FOR FUTURE RESEARCH 


A number of aircraft companies and laboratories are 
now carrying on theoretical and experimental research 
on cruciform wings. Practically none of this material 
is published, and, for the most part, the details of these 
researches are not known to the author. Hence, some 
of the research recommended here is probably already 
under way. 


Theoretical Research 


Most of the theoretical research on cruciform con- 
figurations appears to be by slender body methods. 
Reference 10 gives the solution for combined pitch and 
yaw and shows that, within slender body theory, no 
rolling moments can be induced. Future slender body 
research should then be confined to the aileron deflection 
and damping in roll problems. Such solutions will be 
useful as limiting cases. Furthermore, once the basic 
incompressible cross-flow problem is solved, solutions 
may be given for more general cruciform wing-body 
combinations than considered here. 

Some attempt should be made to solve cruciform 
delta wings with subsonic leading edges. The simplest 
problem appears to be the case of aileron deflection with 
no body, since it can be attacked by conical flow meth- 
ods. For other cases such as damping in roll and cruci- 
form wing-body combinations, the labor involved in 
obtaining exact solutions may not warrant the effort. 

Some study should be made of the nonstationary 
problems. Results nonstationary plane wing 
problems will be useful. 


from 
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Work should be continued on obtaining vorticity dis- 
tributions at the trailing edges of cruciform wing-body 
combinations. The most useful results for the cross- 
flow field at the tail will probably be obtained by studying 
the motion of incompressible vortices as approximations 
to the rolled-up vortex sheets. Some work along this 


line is in progress. 


Experimental Research 


Some research is now in progress on the vortex motion 
and the cross-flow field behind cruciform wing-body 
combinations. More fundamental research is needed, 
especially for those problems where linearized theory 
can yield no answers. 

The really basic job of experimental research should 
be the systematic hunting down of nonlinear and vis- 
cous effects. 
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Nonstationary Motion of Purely 


Supersonic 


Wings 
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ABSTRACT 


A general theory is presented for the calculation of the total 
forces acting on purely supersonic wings. The method applies 
to wings having an arbitrary downwash distribution (stationary 
or nonstationary) and is valid whenever all of the wing edges are 
supersonic. The general three-dimensional nonstationary prob- 
lem is reduced to an equivalent two-dimensional problem. In 
the case of harmonic oscillations, the aerodynamic coefficients 
are expressed in terms of known or tabulated functions. The 
specific example of an oscillating delta wing is considered, and 
valués of the aerodynamic coefficients for plunging, pitching, and 
rolling oscillations are calculated for two Mach Numbers. 


TABLE OF NOTATION 


z = Cartesian space coordinates 


%, 9,3 

ie = time 

¢g = perturbation velocity potential 
Pp = pressure 

p = density 

c = acoustic velocity 


L = lift 


M = pitching moment 
= rolling moment 
For Fe = coordinates of the port and starboard lead- 
ing edges 
b = semichord 
w = frequency 
&, % 8 = wing coordinates 
M = U/c = Mach Number 
U = free-stream velocity 
vy = whb/U = reduced frequency 
k = vM/p? 
B? = M?—1 
o = sweep angle 
S = wing area 
q = (1/2)poU? = dynamic pressure 
C, = L/qS = lift coefficient 
Cy = M/2bqS = pitching moment coefficient 
l = span 
C, = R/qSl = rolling moment coefficient 
a = angle of attack 
p = rate of roll 
Cds = Case a; plunging oscillation 
()s = Case b; pitching oscillation 
Cle = Case c; rolling oscillation 
(a = Case d; pitching oscillation of a flap 


(1) INTRODUCTION 


W"™ THE ADVENT OF high-performance aircraft 
and the possibility of supersonic flight, the 
problem of flutter and aerodynamic stability of airfoils 
moving at supersonic speeds has become of increasing 
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importance. Consequently, much attention has been 
directed to the problem of an oscillating airfoil in a 
supersonic stream. 

As in the case of the subsonic theory developed by 
Glauert,! von Karman and Sears,? and others, the 
approach to the supersonic nonstationary problem has 
generally been made through the use of the linearized 
equations of motion. It is well known that, if one 
considers a nonviscous nonheat-conducting fluid with 
the assumptions of irrotationality and small dis- 
turbances, the equations that govern the disturbed 
motion of the fluid may be reduced to the wave equa- 
tion. Under these conditions, it is possible to find 
general solutions of the wave equation which may be 
superimposed to satisfy the boundary condition of 
tangential flow on the wing surface. 

Possio*® and Borbely* have both treated the two- 
dimensional airfoil in this These authors 
have obtained solutions for the pressures on an airfoil 
due to an arbitrary chordwise downwash distribution 
that exhibits a harmonic time dependence. These 
solutions are expressed as integrals that cannot be 
evaluated in terms of known functions. However, 
Schwarz® has calculated these integrals by numerical 
methods for a sufficient range of the parameters to 
make the use of this theory practical. 

Garrick and Rubinow® and, recently, Miles,’ fol- 
lowing the early work of Possio and Borbely, have 
treated the two-dimensional theory as applied specifi- 
cally to the flutter problem. These authors have 
obtained analytical expressions for the flutter de- 
rivatives in terms of integrals tabulated by Schwarz. 

The work of Garrick and Rubinow, which includes 
some numerical results, indicates that for a certain 
range of Mach Numbers the aerodynamic damping due 
to torsional oscillations becomes negative; hence, the 
motion is unstable. This point is of paramount in- 
terest to the designer, and the question naturally 
arises as to what further instabilities may arise in the 
case of a finite wing. 

Miles® has carried out an approximate solution of 
the problem of an oscillating delta wing, including 
only terms of first order in frequency. This analysis 
indicates that the instabilities found by Garrick and 
Rubinow may also occur for wings of finite span. 

Recently, Stewart and Li*® have obtained the solu- 
tion for an oscillating rectangular wing. They have 
found that there is a marked change in the nature of 
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Fic. 1. 


the damping of certain modes due to the finite aspect 
ratio. . 

The general solution for an oscillating three-dimen- 
sional wing with supersonic edges* may easily be ex- 
pressed in terms of elementary solutions of the wave 
equation. However, even for relatively simple plan 
forms, the calculation of the pressures involves integrals 
much too difficult to evaluate. Fortunately, know- 
ledge of pressure distributions is not of prime im- 
portance in the application of the theory. On the con- 
trary, interest is centered on the nature of the total 
forces acting on the wing. 

Miles" has pointed out that it is possible to reduce 
the problem of calculating the forces acting on an 
oscillating delta wing to an equivalent two-dimensional 
problem.+ It is the purpose of this paper to present 
a general solution for the arbitrary motion of any 
three-dimensional wing having supersonic edges in 
terms of an equivalent two-dimensional problem. In 
the case of harmonic oscillations, the solution is ex- 
pressed in terms of the familiar two-dimensional 
functions tabulated by Schwarz. This reduction is 
possible only if the trailing edge of the wing is straight 
and perpendicular to the direction of the free-stream 
velocity. 

The specific example of an oscillating delta wing is 
considered. The forces acting on the wing due to roll- 
ing, plunging, and pitching oscillations are calculated 
for two Mach Numbers as a function of frequency. 
GENERAL THEORY 


(II) 

Consider a wing contained in the %, ¥ plane which is 
extremely thin in the 2 direction and is moving in the 
direction of the negative * axis with a supersonic 

* It is common practice to refer to a wing boundary as being 
either subsonic or supersonic according to whether the normal 
component of the free-stream Mach Number is greater or less 
than 1. 

+ The publication in which Miles reported the results of his 
analysis!? was brought to the author’s attention after the com- 
pletion of this manuscript. The present paper is in agreement 
with reference 12 on the results for the wide delta wing. 


PURELY SUPERSONIC WINGS 299 


velocity U (Fig. 1). Let the wing perform any motion 
of small amplitude and small velocity about the plane 
z= 0. 

If the flow is assumed to be irrotational, a velocity 
potential may be introduced such that 


u = 09/0 
v = 0¢/d¥y (1) 
w= do os) 


where u, v, and w are the components of the velocity of 
the fluid in the Z, ¥, and 2 directions, respectively. 

Since the disturbances are assumed to be small, the 
equations of motion for an inviscid nonheat-conducting 
fluid may be linearized to give 


079 OfFy . Og 1 0° 
i Se (2) 
az? * aye t Zz? cl? ; 
O¢/Ol = (Po — P)/pol (3 
c? = (dP/dp)o f " 


where Py and pp are the pressure and density far from 
the wing and c is the acoustic velocity. 

The boundary condition on ¢ is expressed in terms of 
the downwash necessary to produce tangential flow 
over the wing. If 2, = 2(%, J, 7) describes the altitude 
of points on the wing at the time 7, then the vertical 
velocity of the fluid adjacent to the wing will be given 
by 


0¢/02(z, 7, 0, 7) = —w/(z, 9, 2) 


where 
—02,,/Ol 
If the edges of the wing are all supersonic, no dis- 


turbances will occur in the plane of the wing off the 
Furthermore, if the trailing 


w(i, ¥, 2) = 


wing except in its wake. 
edge is supersonic, the pressures on the wing are inde- 
pendent of conditions in its wake. Thus, the pressures 
on the wing are completely specified by Eqs. (2) through 
(5). 

The problem represented by this system of equations 
can, in principle, be solved. However, in certain cases 
(harmonic oscillations, for example), the 
encountered in the calculation of the pressures are so 
complicated that they have thus far resisted evaluation. 


integrals 


If only the total forces acting on the wing need be 
calculated, considerable simplification can be obtained. 
It has already been noted that, because of the 
hyperbolic nature of the differential equation and the 
assumption of supersonic edges, the two sides of the 
wing may be treated independently.{ Further, since 
the wing is essentially a lamina, the downwash on the 
wing will be an even function of 2, and, hence, the over 
pressure will be an odd function of 2. Thus, it can 
be readily seen that the total lift acting on the wing 
will be given by 

t This is equivalent to the statement that there are no dis- 
turbances in the plane of the wing off the wing. 
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L= mf fi [(Po — #"), ‘po |dS 


where the integral is evaluated over the upper surface 
of the wing. Introducing Eq. (3), 


—Ui+2b E+UD d, 
L = 2m f f i F (#, 9, 0,2) dy dé (6) 
Ip(F+UT) 


where ¥, and ¥, are the extreme left- and right-hand 
coordinates, respectively, of the wing, 2) is the maxi- 
mum dimension of the wing in the # direction, and 


¥,(0) = ¥(0) = O 


In an analogous manner, the pitching moment about the 
point ¢ = —U? and the rolling moment about the 


& axis can be written as 


—Ui+2b are Ye 
— 2p f (@ + vn f ~ (&, 9, 0, 2) dv dz 
—ulI ip Ol ; 
—Ui+2b 
= =2m f fore y (Z, 4 y, 


Since it is known that 0¢/O/(#, 9, 0, 7) is zero for 
y< J, and J > %,, it is possible to extend the intervals 
of integration in the ¥ direction without affecting the 
Thus, Eqs. (6), (7), and (8) may be written as 


M= 


0, 2) dy dé (8) 


results. 


—Ui+2b _ 
ra) 
L= 2m f / ca (z, 9, O, 2) dy dt (9) 
OT —o Of ’ 


—Ui+2b 
= 2m f (+ UX 
—Ut 
@ 9 
f- = (z, 9, 0, 2) dy dz 


~Ui+26 
R = — 2p J, 2 


The above equations indicate that, if the following 
two-dimensional functions are introduced, 


y2(#, 2,2) = J Fo(4, F, 2, 2) dy 


the total forces and moments acting on the wing may be 
calculated without explicit knowledge of the solution 
for y, provided that the functions y, and ye can be 
evaluated.* For, if Eqs. (12) and (13) are differen- 
tiated with respect to ? and substituted into Eqs. (9), 
(10), and (11), the expressions for the forces become 


mm = 


(10) 


 @, ¥,0,2)d~ydze (11) 


(12) 


(13) 


—Ui+2b 
Pe 
L =2 ’ w, VU, z 
wf, Di (&, 0, 2) di (14) 
~UT+2b 
M = = 26 f (z + Tee (z, 0, ) dé (15) 


* This possibility was suggested by Dr. P. A. Lagerstrom. 


AERONAUTICAL 


SCIENCES-—MAY, 1951 


- f Ui+2b Ove ‘ss 
—2p — (Z, 
‘0 ui ol 


The functions y, and yY2 have an important physical 
Oy, /O/ represents the lift of a spanwise 


R= 0, t) dz (16) 


interpretation. 
element of the wing, and Oy,/0/ represents the rolling 
moment contributed by a spanwise element of the wing, 


It is now necessary to determine, if possible, the 
differential equations that these functions satisfy and 
the appropriate boundary conditions that must be 
This can best be accomplished 
Consider y, 


imposed upon them. 
by examining the derivatives of y; and yo. 


oY, /dz = / (02¢/dz2) dy 
O°, /02" = / (0°y/02") dy (17) 
OY, /ol?” = / (0°y/Ol*) dy 


Since, from physical considerations, infinite velocities 
and pressures in the flow field are not allowed, ¢ is at 
least a piecewise continuous function of the independent 
Therefore, the integral of Eq. (12) exists in 
have 


variables. 
the ordinary 
tinuities in its first derivatives. 
of Eq. (17) must be considered as Stieltjes integrals. 
The derivatives defined by Eq. (17) will thus have 
meaning whenever the Stieltjes integrals exist. 


sense. However, ¢ discon- 


Therefore, the integrals 


may 


For values of #, 2, 7, such that the second derivatives 
of ¢ > exist, WY satisfies the eque ition 


Oy OW | Oy; f x E Oy 1 oo] 
— = — = AY 
ag? * det ct OF -o Lae? da? car |4 


or, from Eq. (2) 
© a2 
— f “dy (18) 
-2 Oy? ~ 


O*y; 4 Oy, 1 Od Wy _ 
Or O02? 3c? OF? 
Since the sidewash must vanish at infinity, the right- 
hand side of Eq. (18) may be integrated to give 
Oo” 0” 1 oO 

a (19) 
oz? zc? OF 
For certain values of %, 2, and ?, terms of Eq. (19) may 
This is to be expected, since Eq. (19) 
The boundary 


not be defined. 
is a hyperbolic differential equation. 
condition on y; may be obtained in a similar manner. 


Consider 
Oy /0Z = f (0¢/02) d¥ 


If the trailing edge of the wing is straight and _ per- 
pendicular to the # axis and if the leading edges are 


supersonic, Eq. (4) gives 


ra) ? " e 
~ (#, 0, %) = = A w(z, 7, 2) d¥ = —wi(%, 2) (20) 


(0O< €+ Ui< 2b) 





whe 
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b+cT2K+Ct 
/ 
/w,=O h 
: x 
- i 2b onl 
tan 1/U 
Fic. 2. 


Also, 


where w(2, ¥, 4) is the downwash on the wing. 
oy,/d2(#, 0,2) = 0 (#+ Ui< 0) 

The solution to Eqs. (19) and (20) is well known. For 

points on the wing, the solution is 


: £+c(i—71) 
c 
T 0 ¥ —c(i—r) 


The region of integration in the #, 7 plane is illustrated 
in Fig. 2. Now consider yo, 


wi(é, 7) dé dr 
V c(t — 7)? — (x — £)? 
(21) 


y,(z, 0,2) = 


O%Yo/d#? = / 5(0°g/0#*) dF 
O°y2/dz? = / 5(0°p/dz*) dy (22) 
O%y./? = / 3(02p/F?) dF 


Here, the presence of 9 in the integrand introduces no 
convergence difficulties because ¢ and all of its deriva- 


tives are zero outside of a finite J interval. From 
Eq. (2), 
DW. Owe 1 OW = 3to 
awn vs -— ¥ = sf dy (23) 
Ox? on 8 6c* OF —a” OV 
and 
Op O'y i 
p= — = =0 [5,2 Ya, 2,0) 
"OF OF U9 7 


The right-hand side of Eq. (23) may thus be written 


x 9 Y 9 
ra 20 
[> tas= ff > — dy = 
~=" OF” -y" Oye 


Therefore, 
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O*f2/OT? + 072/02" — 1/c?(O*Y2/d??) = 0 (24) 
The boundary condition on Y. becomes 
Oye ; : "Fs : ; , 
— (74,0,%) = — yw(z, J, )dv = — w.(Z, 2) 
Oz ip ‘ ; 
(25) 


(0< + Ui< 2b) 


OY2/02(%,0,2) = 0 (&# + Li< 0) 


The solution for y. for points on the wing may be 
written down immediately as 


w2(é, 7) dé dr 
= 26) 
Ver — r)* — @ — €)* 

The original three-dimensional problem has now been 
reduced to an equivalent two-dimensional problem. 
The methods previously used for two-dimensional 
wings can be applied to obtain the forces acting on the 
wing. 

The solutions represented by Eqs. (21) and (26) are 
general and apply to any downwash distribution, such 
as gust loading, unit step loading, or accelerated flight. 
The special case that will be of interest here is the 
downwash distribution that, at any point on the 
wing, has a harmonic time dependence. This condition 
is 


dy/d2(#, ¥, 0,2) = —w(% + Ul, F)e** 


or, from Eq. (20), 


ov ¥s(8+U72) 
1 - : ” 
— f w(t + Ul, ye’ dy = 
Ppl*t+UD 


— (z, 0,7) = 
Oz 
—w(t + Uije’ (27 


and from Eq. (25), 


¥s(#+UF) 
Op» Vs(F+ oe a - 
~ jw(t + Ul, Fe” dj = 
02 Fp(#+U2) 


(#, 0,7) = 
—w2(F + Ul)e* (28) 


With these boundary conditions, the expressions for 
y, and ye. can be simplified. Consider Eq. (21). 
Using Eq. (27), 

- 2 c oO ttE-9 (E+ Ur)e**" dé dr 
W(2, 0, j= : 

T/0 #—c(i—r) 


Vici — 7)? — (€ — &)? 


(29) 
Introducing the Galilean transformation 
x=#+ Ul, y=795 
s=8 t=i1 
and integrating first along — + Ur = constant, Eq. 


(29) yields the well-known integral relation (see, for 


example, reference 4 or 7) 
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jw ~s ro . It 
e k(x — S) —jkM(x-—s ij 
x(x, 0,2) = - | wild). |e Pe eee Zz all 
B 0 b I 
(30) 
where 
M = U/¢ No 
6B? = M?-1 
k = wh/cp 
and J) is the Bessel function of the first kind of zero 
order. Similarly, wh 
jut 9g 
e R(x — s) , 
~jkM(x—s , 
Yo(x, 0, t) = / ws) 44 k ? = ds (31 ) 
B 0 b wh 
The calculation of y,; and yY. has been reduced to the 
evaluation of a single integral. The calculation of the ,, : : , ; ; 
: 8 8 z 7 The leading edge of the wing is described by the 
forces and moments requires another integration. In . : 
7 : equations 
many cases of interest, w;(x) and we(x) are polynomials 
in x.. When this is true, it will be demonstrated that VY = —x/tano, ys = x/tang (32) 
the expressions for the forces and moments reduce \,. ; —" 
5 j : = Now from Eq. (27) 
to the single integrals tabulated by Schwarz. The : 
= P 5 aed . *x/tan o No 
special case of an oscillating delta wing will now be con- eS a 
. wi(x) = Wo dy = 2wox/ tana (33) by 
sidered. —x/tan o : 
ae : Substituting Eq. (33) into Eq. (30), 
(II) THe DELTA WING : 
a : . 2woe ie J , 
The functions y, and yY» will now be used to calculate V(x, 0, ¢) = sG(x — s) ds (34) ost 
oe : ‘ : Btangd Jo 
the aerodynamic forces acting on an oscillating delta 
wing. The plan form of the wing is an isosceles where 
triangle with a vertex angle of t — 20 (Fig. 3). The P C 
rs = ° “ R(x — S) . na a 
vertex is placed at the origin of the x, y, z coordinate G(x — s) = Jo | | “REG-o/e 8 6(3) 
system. The maximum chord of the wing is 20. b 
Ihree types of motion will be considered: plunging, It will often be convenient to express integrals of an 
pitching, and rolling oscillations. the type appearing in Eq. (34) in slightly different ms 
: : form. <A change of the variable of integration estab- 
Case a: Plunging Oscillation 5 i ait ; 5 
lishes the following relation: 

This motion is characterized by a downwash that is wh 
the same for each point on the wing. Analytically, fg f(s)h(x — s) ds = J. f(x — s)h(s) ds (36) - 
this is expressed by 0 

jot With the aid of the Galilean transformation introduced 
Oy/02(x, y, O, t) = —woe , ; ; a an 
in the previous section, Eq. (14) becomes 
or = flo 
L=2 (Oy, /OF) U (Oy,0/0x)] dx (37) 
w(x, ¥) = Wo = constant ™ 5 [ovr + v1 | ' 
or 


Substituting Eq. (34) and using the relation of Eq. (36), 


ol Wo) _eigee” /[* : 
L= 2h +l ore (x — s)G(s) ds dx 
0 Ot Ox] Btang Jo In 
4 pow oer! 72b % : "2b be 
L= " jw (x — s)G(s) ds dx + l G(s) ds dx 
B an o 0 0 0 0 D: 


Since the area of the wing is S = 4b°/tan o and the dynamic pressure is g = 2)pol”, the lift coefficient becomes 


a 2woe™ jw f f la vf” r isd 38) 
= By2p2|7 (x — s)G(s) ds dx + G(s) ds dx (3! 


or 











yy the 


(32) 


(34) 


(35) 
ils of 
erent 
stab- 


(36) 


uced 


(37) 
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It is interesting to note that the C, is independent of the sweep angle o. It will be found that this is true for 


all of the coefficients considered in this analysis. 
Introducing the quasi-steady value for the lift coefficient C,;, = (4woe?) /BU, Eq. (38) becomes 


C 2 
(=) = - b ah fi (x — s)G(s) ds dx + Op f° x G(s) ds dx (39) 


Now, if the variables of integration are changed by x = 2b and s = 2b, then 


C 
& ) = viv foe — )G(2bn) dn dé + 2 fe fo G(2bn) dn dé (40) 
Lo’ a 


vy =wb/U. Finally, Eq. (40) can be integrated by parts to give 


(Cr/Cr) = jv[To — T: + (1/4)T2] + To — (1/2)T; (41) 


“1 
= 2° | s"G(2bs) ds (42) 


The stalling moment about the vertex of the wing is obtained from Eq. (15), 


tpowoe?™ ; 2b Pe *2b ~ ; 
M=-— jw x (x — s)G(s) ds dx + l x G(s) ds dx (43) 
B tan o 0 0 0 0 


where 


where 


Now, Cy = /qS2b, and the quasi-steady velue of Cy is Cy, = — (Swye?) 38U. Therefore, by integrating 
by parts and using Eq. (42), the final expression becomes 
(Cu/Cup)a = jvlTo — (3/4)T1 + (1/16)73] + (3/4)T. — (3/16)T2 (44) 


Since the wing is symmetrical about the x-axis and the downwash is an even function of y for this mode of 


oscillation, the total rolling moment about the x-axis is zero. 


Case b: Pitching Oscillation w(x, y) = (1 — 2jvp)U ao + jvaol/(x/b) (46) 
In this case, the wing is assumed to be pitching about - “ , : ’ 
é ; o : aa rhe first term on the right-hand side of Eq. (46) 
an axis a distance 2ub downstream of the vertex. The ie ° 
: i corresponds to the downwash used in Case a. There- 
instantaneous angle of attack is given by ; ae ; . 
fore, it is only necessary to obtain a solution for the 
rm : ; 
a = ae” (45) parabolic downwash, 
where a) = constant. Therefore, the z coordinate of w, = (aoUx)/b 


any point on the wing is 
From Eq. (27), 
Ze = —al(x — 2ub) 


. si tela aol x 2aol x“ 
and the downwash necessary to produce tangential wi(x) = dy = 
: —x/tane b b tan o 


flow over the wing is 


Ww Sx dz, and, from Eq. (30), 


= l = —ae™ | jw(x — 2ub) + U 
Oz ot + ox = Ljo( m+ Ol JaUei 
yi(x, 0, t) = al “t — s)*G(s) ds (47) 


or » tan o 


Introducing Eq. (47) into = (14), 


bs 2b x 
L, = poco Ue™* | iw fo / (x — s)°G(s) ds dx + 2U f J (x — s)G(s) ds ax | 


Bb tana 


Dividing by gS and integrating by parts, 


Gaer— . ae oe fen re Se oo 
Cry = 3B | jo(To — 37, +57 - ro+in-insin| 
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A comparison of this expression with Eqs. (41) and (44) leads to the following result: 
Crp = (16aye?/38) [3(Cr/Crp)a — 2(Cu/Curo)al (48) 
Since the differential equations that are being solved are linear, it is permissible to superimpose solutions, 
Therefore, from Eq. (46), 
Cry = (1 — 2jvu)(Uao/wo)Cr, + jvCr, 
For this case, the quasi-steady C, is given by 
(Crp)o = (4e0/B)e*™ = (aoU/wo) (Cry) a 
so that 
(Cr/Crp)o = (1 — 2jvu + 4jv)(Cr/Cry)a — (87v/3)(Cur/ Cay) a (49) 


It is now necessary to calculate the moment about the vertex due to the parabolic downwash. From Eqs. (15) 


and (47), 
4 . oU jut 2b '% 2b x 
M, = —- etree E if vf (x — s)*G(s) ds dx + ou f x [ (x — s)G(s) ds ax | 
Bb tan o 0 0 0 / 0 


By a change of variable and integration by parts, 
‘ maa,  #4,..8.. te ese 
Cup = — B |ie(Te - 37 +37 - T)+in-n+57| 


Since 
(Cup)x=0 = (1 — 2jvp)(Ua0/wo)Curq + jvCuy 
the total moment coefficient about the vertex is 
(Cup)x=0 = — (Saoe* 3B)(1 — 2jvu)(Cu ‘Cup)a + JvCuy 
In the case of pitching motions, it is generally desirable to calculate the pitching moment about the axis of 
rotation. If (Cy,), is the pitching moment coefficient about the axis of rotation, then 
(Cup)u = (Cun)e=0 + uCry 
Since the quasi-steady value for (Cy,,), is 


(Cupp = [4a0(3u — 2) 3B ]e™™ 


then 
(*) eS | 24 — 2jvy) (=) — 3p (=) + jv (=) | (50) 
Cur bu 2— 3u Cup a Cro/ C Mo/p 
where 
Cu a? oe = aw~ -_ ; ” Bas . 
— = wu\ 37y — 47, + - T2 ——= 1) + 47% — 37,4 - Ts (51) 
Cur Dp 2 16 4 


In this case the downwash is again symmetrical about the x-axis, and the rolling moment is therefore zero. 


Case c: Rolling Oscillation 
In this case the motion of the wing and, hence, the downwash will be specified by the rate of roll d. Thus, if 


: jt 4° . . . . . °c 
@ = pe™, then w(x, y) = py. Since the downwash is antisymmetrical with respect to the x-axis, the lift and 


pitching moment will be zero. It is therefore only necessary to calculate the rolling moment. 


From Eq. (28), 
x/tane 
w2(x) = f py? dy = 2px*/3 tan’ a (52) 


x/tang¢g 


Introducing this into Eq. (31), 


¥2(x, 0, t) = (2pe?*, 368 tan? o) / (x — s)8G(s) ds (53) 





an 


bu 


an 


wh 


Sin 


Sol 


Cr 





(48) 


lutions, 


(49) 


qs. (15) 


xis of 


(50) 


(51) 


us, if 
and 


(53) 
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The rolling moment is obtained from Eq. (16): 


i . pest *2b be "2b ms 
R= - pot : iw f / (x — s)®G(s) ds dx + su | | (x — s)*G(s) ds ax | 
3B tan?* Fi. 0 0 0 0 


If the rolling moment coefficient is defined as C, = R/qSIl, where / = 4b/tan o = span of wing, then 


2I pe’ a” 1 7 7] *¢ 
C= < mee aie fi / (& — n)8G(2bn) dn dé + 3 | / (€ — )*G(2bn) dy as| 
oO s 0 0 Zt 0 


Since the quasi-steady value for the C; is 





C, = —lpe™/6BU (54) 
ind using Eqs. (42) and (54), together with the results of Case a and Case b, the final result becomes 
(Ci/Ci)e = 12(Cr/Crp)a — 8(Cur/Cato)a — (Cup/Cuto)p (55) 


Note that for a given span the rolling moment coefficient is independent of the sweep angle. 


Case d: Oscillating Flap 
The flap considered consists of the 


The final case investigated is that of a flap oscillating about its hinge axis. 
The chord of the flap is 


aft portion of a delta wing (Fig. 3). The hinge line is perpendicular to the x-axis. 
(1 — »)2b. The portion of the wing forward of the hinge line remains at zero angle of attack. 
Under the above conditions, the expression for the downwash is 


w(x, y) = 0 (x < 2ub) 
w(x, ¥) = apU[1 — jou + (jrx/b)] (2ub < x < 26) 


From Eq. (27), 


wi(x) = 0 (x < 2ubd) 
2agl Ix. Z jvx 
w(x) = 1 — 2jyp + ° (2ub << x < 2b) 
tan o b 
and, from Eq. (30), 
v(x, 0,t) = 0 (x < Qyd) 
2ayl leit . P jvs es ane 
W(x, 0, t) = 1 — 2j7vun +° sG(x — s)ds (2ub< x < 2b) (96) 
B tan o J 2ub b 


but the lift is given by 


"2b 
LE = 2» f [(O"; Ot) + (Oy; Ox) LU] dx 


and Oy,/Ot = jw. Integrating the second term of the equation for lift, 


2b 
i = 24 J f V1 dx + Uy, (26, 0, 0 | 
2ub 


which becomes, on substituting Eq. (56), 


4 poao Ue . a ' jvs i : jvs\ 
_= -} Jo (1 — 2ivyp +°* ) scx —s)dsdx+U (1 — 2jvp + “— }sG(2b — s) ds 
; 2nd J Mud . b 2ub b 


6 tano 
Since the area of the flap is 


S, = 4b7(1 — yu*)/tan o 


some algebraic manipulation leads to the following expression for the C,: 


; Sae™™ f ; 1 oy 
Cr, = 42jv(1 — pw) [1 + 2jv(1 — uw) (x — s)][1 — w)(x — s) + wJG[20(1 — w)s| ds dx + 
B(1 4 u) | 0 0 


“1 
/ [1 + 27(1 — w)(1 — s)J[A — wl — s) + wJG[200 — us] dsy 
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It is now convenient to adopt the following convention. If F = F(b), then F’ = F(b’), where b’ = (1 —un)b 
This convention makes it possible to use expressions previously obtained but with the reduced frequency »’ 
based on the chord of the flap. 

The expression for the C, can now be reduced in the usual manner using Eq. (42). The result, in terms of 


expressions already obtained, is 


Cr 2 j _ s ef my Basic 7 I ° 7 aan ce) tjv'(1 — pg) Cu | a 
( ), Fe acs (7. ia )|+5¢ lela te 3 ta Se. 


CioJa 1+ al” 
where the quasi-steady lift coefficient of the flap is given by C,, = 4ave™,’B. 
The pitching moment about the hinge axis is 
2b 
= =? f [(Oy,/Ot) + U(Ow,/Ox) |(x — 2nd) dx 


The second term can be integrated by parts to give 


"2b "26 
M = = 2 je J x = 2ub)y; dx + 2b1 ’( l u)W(26, Q, t) — ‘Oi vA ax | 


2ub . 


Using Eq. (56), 


tppay Ve? TL sad ‘ . js\ ., 
M=- jw = 2ub) 1 — 27vp + sG(x — s) ds dx + 2b(1 — p)l x 
B tanoe J ub 2ub b 
"2b _ 2b my ; 
‘ Jvs ; : : ]vs ‘ 
J 1 — 27vyp +° sG(2b — s)ds — l Fi [ 1 — 2j7yn +° )scicw — s) ds dx 
J ub 2 b 2ub J 2ub : b 


Proceeding exactly as in the case of the C_,, the final result is 
=) 6 j | pom ; (; a Rie l ‘ Cy ; jv'(1 — u) (Cu ‘| 
x = Shu 7 / + jv’ OF se Ay + 1 y 4 2jn'u)( : 4° ( = (5S 
(( Mo’ d 2+ ul t 2 s 3 i Cuo/a 6 Cy oJ 
where 


Cy = —[4a0(2 + u)/3B(1 + pn) Je 


Here, the Cy, is based on the chord of the flap. 


TABLE 1 
M = 1.25 
(Cases band d; pw = 0.5) 


ob MM R(o ) 1(C) R(E*) (=) R(E) 1(C) R(e") 1(<) 
Ci.J a Crsf « Cmu/ a Cme/ a Ci] 6 Crest Cto/ bu Cm] ou 


0.6 0.9720 —().0996 0.9665 —0.1375 0).9722 —(). 0862 0.9545 — (0.0727 
1.0 0.9263 —(0.1572 0.9120 —(.2110 0.9276 —(0.1277 0.8832 —0.0919 
1.4 0.8662 —0. 2035 0.8413 —().2611 0.8710 —(). 1489 0.7970 —().0742 
1.8 0.8001 — 0.2364 0.7647 —(. 2851 0.8124 —0. 1460 0.7136 —(0.0206 
2.2 0.7359 —(). 2552 0.6922 —(). 2847 0.7581 —(). 1283 0.6500 0.0659 
2.6 0.6798 —().2612 0.6314 — (0.2654 0.7161 —().0952 0.6138 0.1670 
3.4 0.6044 —(). 2444 0.5576 —(0.1997 0).6727 —(0.0152 0.6210 0.3611 
1.2 0.5726 —0.2119 0.5364 —(). 1400 0.6670 0.0510 0.6827 0.4918 
5.0 0.5613 —(0.1798 0.5356 —0.1052 0.6685 0.0968 0.7271 0.5703 
2kM R(e) 1() R(e ) 1(€ ) r(e ) (6 ) 
Cas Cy/e Ciy/a Ciy/ a Cu,/ a Cm,/ 4 
0.6 0.9867 —().0748 0.9673 —0.0481 0.9556 —(0.0491 
1.0 0.9694 —(). 1212 0.9158 —().0510 0.8866 —().0523 
t4 0.9424 —().1612 0.8526 —().0343 0.8039 —(). 0200 
1.8 0.9097 —(). 1936 0.7903 0. 0086 0.7265 0.0479 
3.2 0.8738 —().2180 0.7391 0.0722 0.6680 0.1426 
2.6 0.8367 —(). 23438 0.7051 0.1476 0. 6364 0.2499 
3.4 0.7656 —(). 2460 0. 6869 0.2970 0.6442 0.4479 
4.2 0.7096 — (0). 2373 0.6878 0.4117 0.6841 0.5821 
5.0 0.6675 —().2195 0.6998 0.5016 0. 6927 0.6800 
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TABLE 2 
M = 2.00 
(Cases band d; » = 0.5) 
CL Cx Cw’ Cm Cz CL Cu (e ) 
9 
2eM a=). (Zé), RC), (=). Re), (ce), R(GE) ' Cuy/ wm 
0.6 0.9890 —().0389 0.9868 —().0539 0.9892 0.0274 0.9839 0.1542 
1.0 0.9707 —0.0618 0.9651 —().0832 0.9718 0.0521 0.9591 0.2685 
1.4 0.9461 —(). 0804 0.9360 —(0. 1036 0.9499 0.0853 0.9300 0.3971 
1.8 0.9181 —(0.0937 0.9033 —0.11385 0.9274 27 0.9040 0.5404 
2.2 0.8897 —0.1018 0.8708 —().1127 0.9081 7 0.8880 0.6954 
2.6 0.8637 —0.1032 0.8420 —(). 1025 0.8951 0.8863 0.8560 
3.4 0.8269 —().0927 0.8049 —(). 0638 0.8919 0.9274 1.1739 
4.2 0.8144 —().0708 0.7985 —(). 0205 0.9117 0.9990 1.4240 
5.0 0.8203 —().0468 0.8135 0.0098 0.9342 1.0508 1.6285 
2b n(&) de) r(&) AS) r(S") (S") 
Ci Jc Ci, ¢ Cis 4 Ciy/ a Cm/a Cm,/ a 
0.6 0.9955 —(0.0296 0.9927 0.1848 0.9902 0.2346 
1.0 0.9879 —().0476 0.9813 0.3119 0.9745 0.3966 
1.4 0.9771 —(0.0634 0.9657 ().4442 0.9546 0.5659 
1.8 0.9638 —(0).0764 0.9493 0.5822 0.9338 0.7431 
2.2 0.9488 —(0. 0863 0.9340 0.7254 0.9155 0.9273 
2.6 0.9330 —().0928 0.9215 0.8726 0.9018 1.1163 
3.4 0.8987 —0.0961 0.9078 1.1716 0.8915 1.4985 
4.2 0.8758 —()_ 0892 0.9052 1.4661 0.8959 1.8684 
5.0 0.8572 —().0765 0.9035 1.7523 0.8980 2.2249 
IV *ALCULATION AND RESULTS l l l 
va. Sa T(k, M) = J. ( 2M) — iJ, ( 2M) 
kM M ; M 


The aerodynamic coefficients obtained in the previous 
section have all been expressed in terms of the functions 
T,(k, \f). A recurrence formula may be obtained for 
T, by integrating Eq. (42) by parts. This formula 


was first given by Borbely;* however, the RTP 
Translation of the paper by Borbely is in error. The 
correct formula is 
: ® "| 5 J (Qk I 1, (2k 4 —_ mo?) 
5 = yA - o(Z = J ZR) 
‘ M~’ 2kM 
2jkM 4 41 — Qn)T7 oe ‘Shall i, (59) 
' =e kM 7 


where Jy) and J; are Bessel functions of the first kind 
of zero and first order, respectively. With this formula 
it is possible to express all of the results for the delta 
wing in terms of the single function 7). Schwarz has 
tabulated the values of two functions, /,;(A, x) and J,(X, 


(60) 


The coefficients of the previous section have been 
evaluated by means of Eqs. (59) and (60) for two values 
of Mach Number, These results are pre- 
sented in Tables | and 2 and are plotted as polar dia- 


25 and 2. 


grams in Figs. 4 to 10. 

In Fig. 7, at J = 1.25, the phase angle of (Cy /Cy,,)» 
is negative for small values of the reduced frequency. 
Since the damping in pitch is proportional to 2 — 3y 
times the sine of the phase angle, the pitching motion 
is unstable in this range. 

It is of interest to determine the Mach Number at 
which the pitching oscillations first become unstable 
for small values of frequency. This may be done by 
obtaining a Taylor expansion of the 7,. Thus, neglect- 
ing terms of order k’, 












































x), which are related to 7) by the equation T, = [2"*°/(n + 1)] — j[2"*°/(m + 2)|RM (61) 
“S M=2.0 —— = M-2.0/ —— 
O ' oO Rie es 

02 a, M=125 — — 02 Ss! |Me125) 
S 3) 

0 06 07 08 09 LO 0 06 O7 ae 09 10 
R(C, /C, ) ~~ ig/ OMk-= R(Cyhy,) 34°_ 107 
a Ly, a 34° - aan J | "501 : ‘ | ~~, _——1— + 

2Mk=, 50 26 18 7 art 2. is 

be he Lire “02 = 7 pe 

0 2 tl 4 6 Is LL — 10 7. ~ 26 _ IS a - LO 

-04 -04 

Fic. 4. Lift coefficient, Case a. Fic. 5. Pitching moment coefficient, Case a 
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Key, 
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| | [ls | on 2Mk= | thes 
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Fic. 8. Rolling moment coefficient, Case c. 
Fic. 6. Lift coefficient, Case b. 
2Mk =5.0 
bad ] J 
05.0 16 MAI25 — ce | J 
1.6F IM=20 | | | ; x05 
M25 — — | 
205 
1.4 ; 14 | 
12k vA l2L | 4 
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3s > | 
LOF = LOH | 
oO 
~N 
J 126 6 
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on 2Mk= 5.0 on 06 18 
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| (34 | ” | 
\ “1.0 4 4 ] 10 
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es \ | | foe oe : 
LJ | \ | | 26 
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fh ee 06 07 og~ ~~ -09, 5 10 
Fic. 7. Pitching moment coefficient, Case b. Fic. 9. Lift coefficient, Case d. 
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By means of Eq. (61), the expression for (Cy/Cryp)o, 
becomes, for small values of frequency, 


ik 
(=) =]+4 = —~ xX 
Cup by M (3p — 2) 
86° — 6B*u* — 38° +5 — 2u) (62) 
The condition that the damping is positive is 


862u — 6B%u2 — 3B? + (3/2) — 2u< 0 


In Fig. 11 the Mach Number at which the damping 
reverses sign is plotted as a function of yw. A similar 
region of instability was found in the two-dimensional 
case calculated by Garrick and Rubinow.® Fig. 11 
shows that the maximum Mach Number for instability 
is V2 and that all instability vanishes for rotations 
about axes aft of the three quarters chord. These 
results agree with thoSe obtained by Miles.* 

Fig. 10 indicates that the damping in the case of the 























flap also becomes negative at M = 1.25. By using 
Eq. (61), the expression for (Cy,/Cy))a becomes 
Cu ) jk’ | ( ‘) | 
—— = 1+ — M? — + 3M — 
on d (2 4 p)M . 2 Z 
(63) 
Ler e20|— 
M=125|— — a 
205 34 
“ | 
: | 
12 Z 
Ss. 026 
oO | 
lol — .| 
\ 
\ 
\ 
08 \ 
18 
\ 4 
2Mk= 7 5.0 \ 
06 ) \ 
= \ 
04 / 10 
126 
| 
0.2 
18 | 
06 0.7 '~ 08 09 1}0 
R(C,,/ ae ‘La 
Fic. 10. Pitching moment coefficient, Case d. 
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Fic. 11. Stability regions for pitching oscillation. 
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Fic. 12. Stability regions for oscillating flap. 


Thus, the condition for positive damping is 
u[M?2 — (7/2)] + 3M? — (9/2) > 0 


The Mach Number for zero damping is plotted as a 
function of wu in Fig. 12. It is seen that, for a flap of 
infinite aspect ratio (u = 1), the Mach Number for 
zero damping is 1/2. This is the two-dimensional 
result obtained by Garrick and Rubinow. 

(V) CONCLUSIONS 

Using the linearized equations for a compressible 
isentropic flow, one can reduce any three-dimensional 
wing problem to an equivalent two-dimensional 
problem, provided that the edges of the wing are super- 
sonic and that the trailing edge is normal to the flow. 
For the usual types of plan forms and downwash 
distributions, the forces acting on the wing can be 
expressed in terms of tabulated functions. 

As an example, the forces acting on a delta wing have 
been calculated for the case in which the downwash 
exhibits a harmonic time dependence. All of the 
aerodynamic coefficients were found to be independent 
of the sweep angle. Since any nonstationary motion 
of a rigid delta wing can be obtained by harmonic 
analysis from the three modes of oscillation considered 
here, it is proved that the force coefficients for such an 
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arbitrary motion are independent of the sweep angle. 
It is also found that the pitching oscillations of a delta 
wing exhibit a characteristic instability for small values 
of the frequency. This is in agreement with the two- 
dimensional result obtained previously. 

Finally, it should be noted that, while it is correct 
to calculate the total loads acting on the wing by strip 
theory, this procedure does not necessarily give the 
correct spanwise loading.* It is interesting, however, 
that, although the load distribution is disturbed by the 
three-dimensional effect, the two-dimensional lift of a 
chordwise strip is preserved. This result is in accord 
with the theorem on the preservation of lift for steadily 
lifting wings given by Dr. P. A. Lagerstrom.!! 
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Symposium on Fluid Dynamics 


The Fourth Symposium on Applied Mathematics of the American Mathematical Society 
(cosponsored by the University of Maryland and the U.S. Naval Ordnance Laboratory) will be 
held at the University of Maryland, College Park, Md., and the U.S.N.O.L., White Oak, Md., 


The general topic of the Symposium is “Fluid Dynamics,’ and sessions are planned on 
Turbulence, Foundations, Compressible and Incompressible Flow. 


I.A.S. members interested in attending may secure additional information by writing to 
M. H. Martin, University of Maryland, College Park, Md. 
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Effects of Pressure Gradient on Stability and 
Skin Friction in Laminar Boundary Layers 
in Compressible Fluids' 


HERSCHEL WEIL? 


General Electric Company 


ABSTRACT 


An extension of the von Karman-Pohlhausen method is carried 
through to find approximate steady-state velocities and tem- 
peratures in the compressible laminar boundary layer with a pres- 
sure gradient in the direction of flow along a two-dimensional 
adiabatic wall. These are used to find the effect of pressure gra- 
dient on skin friction; om momentum thickness; and, by a cri- 
terion due to Lees, on stability to small perturbations. 


At a given point along the airfoil, the boundary layer is stable 
if the Reynolds number Rg based on momentum thickness is less 
than a certain critical value (Ra)cy. min. which is a function 
solely of Mach Number Mand a suitable parameter \ dependent 
on the airfoil shape. A curve family of (Re)cr. min. vs. \ for var- 
ious M’s is given. In addition, the product of skin friction times 
Ro, which depends solely on X, is plotted. Numerical results 
were obtained for flows at M7 = 1.5 and 4 over a 5 per cent thick 
biconvex airfoil which showed that for these the 
favorable effect of the negative pressure gradient on stability is 
not large for M = 4 but is extremely pronounced at the lower 
supersonic Mach Number MM = 1.5, and in both cases there is an 
increase in drag up to about 30 per cent at mid-chord compared 
to the case of a flat plate. The dependence on M and \X of the 
conversion factors between Reynolds Numbers based, respec- 
tively, on distance from the leading edge x, on momentum thick- 
ness 6, and on displacement thickness 6* is exhibited. 


examples 


NOTATION 


Cy = local skin-friction coefficient 

( = specific heat at constant pressure 

1 = enthalpy 

k = [see Eq. (29)] 

K = [see Eq. (10)} 

L = a typical length, such as chord length 

M = Mach Number 

n = distance measured normal to airfoil boundary 
p = pressure 

p = p/p-« 

P, = Prandtl Number 

R = universal gas constant or Reynolds Number 


Reynolds Numbers defined, respectively, with 
x, 6*, and @ as typical length, u,x/», etc. 

distance measured along airfoil boundary from 
the leading edge 

t = [see Eq. (11)] 

y = temperature, absolute, (7 = 7/T)) 
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u = velocity in s direction (tangential) 

7 = velocity in m direction (normal) 

w = u/uy 

x = distance measured along airfoil chord from the 
leading edge 

Y = ratio of specific heat at constant pressure to spe- 
cific heat at constant volume 

6 = boundary-layer thickness 

6* = boundary-layer displacement thickness 

n = n/é 

4 = boundary-layer momentum thickness 

K = [see Eq. (32)] 

r = [see Eq. (16)] 

Mm = viscosity 

p = density 

v = p/p 


A superscribed bar denotes a dimensionless ratio as given by 


Eq. (2). 


Subscripts 

( )o = stagnation value 

(); = “local ambient’’ value—i.e., the value taken just 
outside the boundary layer at the s value in 
question 

( = undisturbed free-stream quantity 

(er. min. = critical minimum 

()w5 = value taken at the solid boundary 

(de = transition 

(I) INTRODUCTION 


rr ORDER TO INVESTIGATE both skin friction and the 
stability to time dependent perturbations of a 
laminar boundary layer in compressible flow, it is 
necessary to obtain the steady-state velocity and tem- 
perature distributions in the boundary layer. An ex- 
tension of the von Karman-Pohlhausen method of in- 
compressible flowtt is given here to find these steady- 
state distributions with sufficient accuracy for stability 
computations for flows over insulated bodies with a 
generally nonvanishing pressure gradient in the direc- 
tion of flow. It is assumed that the Prandtl Number 
is unity and that viscosity is proportional to tempera- 
ture. The velocity and temperature distributions ob- 
tained in this way are then applied to studying the sta- 
bility of laminar boundary layers by the methods of 
Lees and Lin! and to finding skin friction. 

In connection with stability, the formulas are given 
for the computation of a particular Reynolds Number, 


tt This extension is closely related to the work of Dorodnitzyn® 
and Howarth’ as explained in Section (III) of this paper. 
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Fic, 2. 


(Ro)cr. min.» by Lees’? rapid approximate method when 
the steady-state distributions are given in a special 
coordinate system that is used here. The significance 
of (Ro)er. min. lies in the fact that, in laminar boundary 
layers at points where the Reynolds Number R, < 
(Ro) cr. min.» Small time-dependent disturbances orig- 
inating in the layer will be damped out; the bound- 
ary layer will be stable, whereas at points where Ry > 
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(Re)cr. min. these disturbances may build up and cause 
the boundary-layer flow to become turbulent. 


The quantity (Rg)cr. min. computed with the type of 
velocity distributions considered here depends solely 
on the local ambient} Mach Number and a quantity i, 
called the form parameter. For a given local ambient 
Mach Number, J is proportional both to the local am- 
bient velocity gradient in the direction of the undis- 
turbed flow: and to the square of the boundary-layer 
thickness at the point in question. Thus, the informa- 
tion on the dependence of stability on velocity gradient 
or pressure gradient may be essentially contained in a 
family of curves of (Rg)cr. min. VS. \ for various values of 
M. However, for practical purposes, the quantity 
(Rz)er. min. iS Of greater importance than (Rg)er. min, 
since it is (R,)¢;. min. that determines how far behind the 
leading edge the boundary layer is stable—namely, 
from the leading edge to the point where R, becomes 
equal to (R,)er. min. Thus, one is also interested in the 
effect of pressure gradient on the conversion factor 
k = R,/R,?. Unfortunately, k does not depend ex- 
plicitly on M and X only but also depends explicitly on 
position, %, along the airfoil and the local ambient ve- 
locity gradient. Thus, universal curve families of these 
quantities are not possible as with (Rg) ¢-. min... The situ- 
ation is different from the flat plate case where, within 
the approximations used here, & is independent of M 
and x. Similar remarks hold for the quantities 
C/V us Vy and @V u/s related, respectively, to skin 
friction coefficient and momentum thickness, since 
these are shown to be respectively proportional to Vk 
and equal to 1/% ‘k. The quantity C;,R,, however, 
depends solely on X. 


For high-speed flights of noninsulated bodies, where 
heat is conducted from the air to the solid body, 
(Ro)cr. min. iS actually larger than would be predicted 
for the insulated bodies considered here. The bound- 
ary layers in this case are consequently more stable, 
so that the values of (Ry)c-. min. computed by the meth- 
ods given here may be looked on, in a sense, as worst 
possible values. The increase of (Rg)er. min. with de- 
crease of surface temperature is actually a pronounced 
effect. In fact, from the data in Fig. 6b of reference 2 
for flow over a flat plate at Mach Number 0.7 with no 
pressure gradient, (R,)¢r. min. increases from 3.8 X 10° 
to 1.8 X 10° as the wall temperature decreases from 
that of an insulated wall, 7\,,, to 0.77 Tow. 


At supersonic speeds, this effect is even more pro- 
nounced, since, at certain critical temperatures some- 
what lower than adiabatic wall temperature, complete 
stability [(Rz)er. min. = © ] is achieved over a flat plate, 
M > 3. Forexample, at \/ = 4, this critical tempera- 
ture is computed in reference 2 to be 0.86 Ty. 

+ The term “local ambient” is used throughout this report to 
refer to values just outside the boundary layer at an arbitrary 


distance, x, from the leading edge. 
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TABLE 1 
Configuration 5 M, dp/ds PN k (Re Jer. min din: ote 
5 per cent biconvex airfoil 0.0789 1.33 —0.687 0.810 2.48 107 28,400 
0.500 1.5 —0.562 4.72 2.75 Complete stability 
Flat plate rae 1.5 0 0 2.29 55.0 6,930 





(II) RESULTS AND CONCLUSIONS 


(A) Stability 


Curves were computed of (Rg)er. min. VS. X at Mach 
Numbers 0, 1.5, and 4 and are shown in Fig. 1. The 
curve for 1 = 0 was computed from the data in Fig. 2 
of reference 3 by using Eq. (32) of this report, while the 
other two curves were computed by the methods of 
Section (IV).t For M = 1.5, the favorable effect of 
negative pressure gradient is pronounced, (Rg)er. min. 
becoming infinite at \ = 4.38. However, for VW = 4, 
the favorable effect of negative pressure gradient is in- 
significant. 

It is evident from Fig. 1 that, in the supersonic 
range, for a fixed A, this favorable effect of negative 
pressure gradients becomes much less as the Mach 


Number increases—that is, the ratio 


(ita) a. min. (A, M), (Ro )er. sta tO, M) 


for a fixed positive \ decreases as \/ increases. 


Computations were carried out to find the local 
ambient values of Mach Number and of \ along the 
forward half of a two-dimensional biconvex airfoil of 
5 per cent thickness-to-chord ratio and undisturbed 
free-stream Mach Number 1.5. The dA values are given 
by the curve labeled “integral curve” in Fig. 5. Ata 
point 7.58 per cent of the chord length from the leading 
edge, the local ambient Mach Number is 1.33, \ = 0.810, 
and the value (Rg)cr. min. is 107, lying somewhat above 
the corresponding point on the M@ = 1.5 curve of Fig. 
1. At the mid-chord point 7 = 1.5, = 4, and there 
is complete stability. These results are included in 
Table 1. This shows that, near the leading edge, pres- 
sure gradient is not sufficiently negative to bring an 
important increase of (Rg)er. min. over the value for a 
flat plate, but it is evident that the extremely strong 
favorable effect of a sufficiently increased \ is brought 
into play as one nears the mid-chord point. 

Table 1 contains, in addition to the values of \ and 
(Re)cr. min. Mentioned above for the airfoil, the cor- 
responding values of k = R,/R,? and (Rz)er. min. and 
also the values applicable for a flat plate at MM. = 
1.5. The k values were computed by Eq. (29). 

Bear in mind that these results are for an insulated 
(adiabatic) wall and that, at slightly lower wall tem- 





+ The author is indebted to Prof. Lester Lees, of Princeton 
University, for pointing out to him that in reference 9, on which 
this paper is based, the stability curves for M = 1.5 and 1.4 were 
computed incorrectly. This came about because a factor 7 was 
omitted from the equation in that report corresponding to Eq 
(40) in this paper and on which these curves are based. 


peratures than that of an adiabatic wall, finite R.,. min. 
may increase radically as discussed in Section (I). 


In addition, it should be pointed out that a differ- 
ence is to be expected between R,,. »in., the Reynolds 
Number at which instability may occur, and R,,, the 
Reynolds Number at which the originally laminar 
boundary layer may actually become turbulent. At 
present, there is no theory to determine R,,, although 
R,, has been found in a few cases by experiments where 
sufficiently low free-stream turbulence was obtained 
so that transition due to disturbances originating in the 
boundary layer could be observed.‘~* These experi- 
ments indicate that, for adiabatic conditions, (R,),, is 
never less than about 10°. Comparison with the 
(Rz)er. min. Values corresponding to the curves of Fig. 1 
shows that the magnitude of the difference Ro, min. 
R,, is far greater than the increase of Re. min. with favor- 
able pressure gradient for a sufficiently high Mach 
Number such as = 4. This indicates that the pres- 
sure gradient effect on Re,. min. is relatively unimportant 
at extremely high speeds and can safely be neglected. 


Incidental to the stability computations, it was neces- 
sary to compute the data needed for curves of w vs. »— 
that is, velocity divided by local ambient velocity vs. 
distance normal to the airfoil divided by the local 
boundary-layer thickness. Two of these curves are 
shown in Fig. 2—one for flow over a flat plate and the 
other for flow at the mid-chord point of a 5 per cent 
airfoil, both for M = 1.5. 


(B) Friction Drag 


The formulas of Section (III C) show how to com- 
pute not only k but the closely related quantities, 
momentum thickness and skin friction. The numeri- 
cal computations showed an increase in local skin- 
friction coefficient as compared to the case of a flat 
plate at the same free-stream Mach Number. This in- 
crease reached about 30 per cent at the mid-chord 
point for M,, = 1.5 as the results in Table 2 indicate. 


In addition, C,R,, which depends solely on \ (Eq. 
50), is shown graphically as a function of in Fig. 3. 


TABLE 2 


OV m/nx = 


Configuration M, § 1/k Cy Vux/n 
5 per cent bicon- 
vex airfoil 1.33 0.0789 0.635 0.736 
1.5 0.500 0.604 0.975 
3.68 0.0789 0.654 0.704 
4 0.500 0.683 0.883 
Flat plate Oto Otol 0.661 0.661 
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Fic. 3. Skin-friction coefficient times Reynolds Number based 
on momentum thickness plotted vs. form parameter. 
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Principally as a check on the methods used here, it 
was decided to verify some numerical results obtained 
by Howarth’ by a procedure not related to the von 
K4rman-Pohlhausen method for a flow given by % = 
Bz. This is actually the type of flow one gets near the 
forward stagnation point on a blunt-nosed body. Fig. 
6 of reference 7 is, in the present notation, a curve 
of 


= 3 
1+7— um 


1 f 7 
= C;V R01 — Bz?) vs. My, 


It is reproduced in Fig. 4 for the abscissa range 0 to 
0.6 (1 < M, < 0.67), together with the corresponding 
curve computed by the procedures of Section (III) but 
using fourth-order polynomials. 


THEORETICAL BASIS FOR THE STEADY-STATE 
FLOw COMPUTATIONS 


(IIT) 


(A) General Discussion 


The method used here to determine the steady-state 
velocity distribution in the boundary layer is an ex- 
tension to compressible flow of the von Karman-Pohl- 
hausen method used in incompressible flow problems. 
This method consists of assuming that the velocity u 
in the integrated momentum equation is a suitable 
polynomial in the perpendicular distance n from the 
solid body. The polynomial involves an undetermined 
parameter, A(s) occurring linearly in the coefficients. 
By this procedure the integrated momentum equation 
is reduced to a first-order ordinary differential equation 


1951 


for } which may then be integrated under a suitable 
boundary condition. 

It will be shown in Sections (III B) and (III C) that 
the extension to compressible flow is made feasible by as- 
suming P, = K = | and by introducing a certain change 
from #7 to a new variable 7. This transformation, Eqs. 
(11) and (12), was introduced by Dorodnitzyn® simul- 
taneously with an additional but nonessential transfor- 
mation in the s direction which is actually a hindrance for 
numerical application of the theory. Later, Howarth’ 
applied a transformation in the normal direction only, as 
in this report, his new coordinate being a function of s 
times the / of this report. His reference quantities are 
not the same as the ones used here, and his JX is defined 
in a slightly different way. Since Howarth used only a 
transformation in the normal direction, there is no 
basic reason to prefer the transformation used here to 
his transformation, and his \ and differential equation 
for \ would serve as well. However, the work reported 
here was largely completed before reference 7 came to 
the writer’s attention. 

Both these authors considered fourth-order poly- 
nomial approximations for “« which permit the boundary 
conditions on u and its first two derivatives with respect 
to n to be satisfied. However, for the stability prob- 
lem, the second derivatives of u with respect to n 
are extremely important, and, in order to obtain these 
with some accuracy, a sixth-order polynomial is neces- 
sary. One can then impose, on the third derivative of 
u, the proper boundary conditions at the wall and at 
the outer edge of the boundary layer so that at least the 
correct rate of change of the second derivative can be 
ensured near these boundaries. It turns out that the 
transformation (11) is essential not only in obtaining the 
ordinary differential equation for \ but also to enable one 
to satisfy the boundary condition on the third deriv- 
ative at the solid boundary. 

In the following, the derivation of the differential 
equation for ) is different from the equation in refer- 
ence 7 or 8, where the transformations are applied to 
the nonintegrated momentum equation expressed in 
terms of a stream function, and thereby it is attempted 
































1.4 —— 
I 
1 2,4 
ed a oe 
(§=4x) 
1.0 ee 
O 0.2 0.4 0.6 


m, {i+ 2m 


—o— METHOD OF SECTION ILL 
—— HOWARTH 
Fic. 4. 











wl 
th 
pa 
th 
ho 
an 


of 


on 


Th 
the 
cor 


PCy 


ma 


Wh 


and 


itable 


) that 
by as- 
hange 
, Eqs, 
simul- 
nstor- 
ice for 
rarth’ 
ly, as 
1 of 5 
eS are 
fined 
mly a 
is no 
‘re to 
ation 
orted 
ne to 


poly- 
dary 
spect 
prob- 
to n 
these 
eces- 
ve of 
id at 
t the 
in be 
t the 
g the 
> one 
eriv- 


ntial 
efer- 
d to 
d in 
pted 





' 


; 














LAMINAR 


to make the reason for the change from n to + more 
evident than it is in these references. 


(B) Transformation of the Integrated Momentum 
Equation to an Equation for 


The well-known momentum equation in boundary- 


laver flow is 





( ou 1 =) Op 4 a) ( 2) (1) 
u oj =-—- — — 
4 Os on Os on “ on 
In terms of the dimensionless variables 
- u - 5 n ; 6 (2) 
a = = §=-, =—, =— 2 
V9, ae L 


Eq. (1) may be integrated for compressible flow to 


m Ou Pa) Sp wwe 
( a *) = — / tN (m, — ut) dn + 
po V2it On]; = 0. OF So po 


Ou 3 
a (2 —— aan (3)t 
OS 0 \po Po 


where the integration is carried out at a fixed §. In 
the incompressible case, p is constant so that, when a 
particular form of % such as a polynomial is assumed, 
the integrals may be immediately evaluated. Here, 
however, the variable p appears under the integral sign, 
and one must know the behavior of p as a function 
of n. 


The assumption of P, = | takes care of this, for then 
I 


one can show that 


p = \7/(y — 1) - 
= (1 — #,*)""*’ D/(1 — #) (4) 


Po 


This is a consequence of the fact that, with P, = 1, 
the momentum equation [Eq. (1)] and the equation of 
conservation of energy in the boundary layer" 


, ( oT 4 2) Op 1 (“) 4 ra) (« or) 
“u— yv—J=u — — 
Os on Os ” on on On 
(5) 
may be combined to yield 
( ra) r °)( T+ ) °| re) ( r4 “) 
u +yv— = = — IG 
' Os On “ 2 Os 2 on " 2 
(6) 


An integral of this corresponding to an adiabatic wall 
[for which (O7/0n),-) = O] iscpT + (u?/2) = iy = 
constant, which may be rewritten in the following form: 


Tr = 7,(1 — #) (7) 
When this is combined with the equation of state 
pb = pRT (8) 
and Bernoulli's equation, written in the form 
p = poll — a2)” -” (9) 


t See reference 9 for details of this. 
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Eq. (4) follows. 
of perfect gases, p = constant p’. 
In addition, one can assume that 


Eq. (9) makes use of the assumption 


w= pw(1 - u2)* (10) 


where AK is a numerical exponent that will later be 
specialized to K = 1, 

Now, however, if Eq. (4) is substituted in Eq. (3) 
and fourth- or sixth-order polynomials in n are assumed 
for a, the integrals cannot be readily evaluated as they 
would be in Eq. (3) with p = constant because there 
will be terms with factors (1 — #?)~! in the integrals. 
This suggests the introduction of the new variable 


" n — 2\¥/(y¥ -— 
p (1 — 7) 
ft = f — dan = f _— - dn 
0 Po 0 1 — n* 


In addition, one can make the range of integration 
independent of 5 by using, instead of ¢, 





(11) 


r = t/t(6) (12) 


Eq. (3) takes the form 
1d [t(6) |? ; - 
nl H() (/, a2 dr — ty f ii ar) + [t(5)]? X 
2 ds 0 0 
d 1 d . ny diy 
72 1 a ae u dr — —-— x 
E f sila l i ds I aaa 1 — a,” ds 


ai (1 — iy?) ~ vy dit 
(1 — #) dr| + — 7 a = 0 
0 V 2 iL dr r= 0 


(13) 


from which it is evident that, when @ is a polynomial 
in 7, the integrals may be readily evaluated. 

The conditions that the velocity 7 must satisfy are, 
first of all, 


On = O7*n a orn) = 


_ | (14) 
Or, = | Or*|, = | Or’, = | 


u(1) = ty, 
on the usual assumption of a flow outside the boundary 
layer which is independent of y. Because the fluid 
does not slip along the solid boundary, 


u(0) = O (15) 


Two more boundary conditions may be determined at 
7 = 0 by considering the limit as y — 0 of the momen- 
tum equation and also the limit of its y derivative. 
These two limits give, respectively, the conditions 


(0°%/0r?)|, -9 = — mr 
{ (4) PL V 2% 


L Vin \ ie ae 
= V0 — ae Def Gs = 
f__ — 19) P pol. 


12 | 
uous (1 — 02)?7 — PY ds 


(16) 


and 


(0°u /dr?) = 2(K — 1) (On 0r)4|, oa 


One can check that near a stagnation point these boun- 
dary conditions reduce to those in the incompressible 





r=0 
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case. Eqs. (14) to (17) form a set of equations which 
determines the coefficients of the polynomial i(r). 
However, Eq. (17) makes the set nonlinear and there- 
fore not readily solved unless K = 1, which is the rea- 
son for making this assumption. When it is made, the 
set of conditions on # and, consequently, the polyno- 
mial for + take the same form as found by Schlichting 
and Ulrich"! in the incompressible case. 


ua = my [f(r) + dg(r)] (18) 
where 
f = 2r — 5r* + Gr’ — 27° ) 
1 1 3 8) 
gain petri et ie | (19) 


It is interesting to note that the boundary condition 
corresponding to Eq. (17) expressed in terms of deriv- 
atives with respect to 7 is 


(d*i/dn*)|, -9 = 2K (dii/dn)*q = 0 (20) 


Thus, if # were assumed to be a polynomial in 7 
rather than in 7, linear equations for the polynomial 
coefficients would be possible only when K = 0, so 
that » would be everywhere constant, which is less 
realistic than the assumption u/u = T/T) made 
above. 

With the help of the definition [Eq. (16)] of A, Eq. 
(13) may now be written as an equation for \ which is to 
be solved subject to a proper boundary condition at 
some point s = sp. This equation may be manipulated 
to the form 


dd/ds = U,A(A) + UedAs(r) (21) 


th," Qty’ 
U= 4,4 — “(24 2 :) 
ty 1 — uy,’ y¥—Il1/ ¢ 
U2 (thy’/ty) oa [iiyty’ ‘(1 = ii") | 


while the functions A; and A» are the same as in the in- 


compressible case, 


where 


II 





A\(A) = 4A(A)/[A(A) + BO)] 
Ao(X) _ \ (23) 
—2{ac(n) + 204d) + [fF - AC —e)b-o} ( 

A(x) + B(A) J 
where 


1 ) 
A(A) = f (f? + 2Afg + A\*g? — f — Ag) dr 
0 | 
1 | 
B(A) = 2d [ (2\g? — g + 2fg) dr 24) 
J 0 


1 
C(x) = I (f + dg) dr | 
0 


When the sixth-order polynomials [Eqs. (19)] are 
used, A, B, and C may be determined numerically, 


and 


1951] 


a Q5d2 + 28n3 .% ) 
“AS = (905 — 25k — 350" | 
(25) 
180,180 — 27,422 + 1,048d2 + 28n8 | 
te se : Sane - aomae — 


Ae  deaiaiae pe eS 
4,925 — 142.5\ — 35d? J 





Note that U; and U; depend only on local ambient 
values and are actually independent of the form of 
function used for a#(r), although A; and A» depend in- 
directly on #(7). Asacheck on Eq. (22), one can verify 
that, in the neighborhood of a stagnation point where 
i, = 0, Eq. (21) reduces to the corresponding equation 
for incompressible flow. 


(C) Momentum Thickness, Skin Friction, and the 


Relations Between Ro and R., and Re and R5*t 


The simple approximate relations 


uy ‘ Ux oe 
6 = C; = as 0.6607 ’ 
Vix Vv} 


used, for example, in reference 2 no longer hold when 
the pressure gradient does not vanish. However, 
when # is approximated by a polynomial, as in the fore- 
going sections, expressions for @ and Cy, as well as the 
connection between R, and R,, are not difficult to 


R, = 2.25R,? (26) 


compute. 
By the definitions of the Reynolds Numbers, 


R, = R,(0/x) (27) 


while the momentum thickness 6 is defined by 


l . 
deca f pu(u, — u) dn (28) 
pill,” 0 


With the help of the formulas of the preceding section 
[with K = 1 in Eq. (10)], one can find, by direct sub- 
stitution in Eqs. (27) and (28), that 


u 1 ) 
oe? — =-, R, = kR,? 
Vix k | 
-~ -— 7 . (29) 
ru 
k (%, d) = ~ | 


(1 — m2)mAA2A) | 


where A is given by Eq. (24). 
In a similar manner, the displacement thickness 


] 6 
f (piu; — pu) dn (30) 
pil J 0 


+ By generalizing the von Karman-Pohlhausen approach for 
use with compressible flow in a different manner from this paper, 
Oswatisch and Wieghardt!? obtain the quantity 6*/@ as a func- 


tion of 


6* = 


> 2/; ) : 
* = (e,,9 /it,,) (du;/ds) 


by a numerical method based on the use of the Hartree velocity 
profiles for incompressible flow. They also obtain an expression 
relating dé*/ds to 6*/@ and the free-stream velocity distribution. 
This expression is applied to the study of stability in the inter- 
action between the boundary layer and the free-stream flow, al- 
though this study is limited to the flat-plate case. References 
to other related German work are given in this report. 
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may be written as 


. Lt(8) E — Uy" (A + C) “fs ; | 
Pw ing-ay’* @-aa-* 
(31) 
where A and C are defined by Eq. (24). Then, 
R;* 5* At? + C-1 
‘c=- = =- (32) 
R, 6 A (1 = ul) 
Finally, the skin-friction coefficient 
Ge = 2 uo (Ou on)|,, = 0, piu; (33) 


may be put in the form 


x 1 1g 
Cy wm =~ DAN) Vk (? + ‘£) (34) 
V1 1 r=0 


dr oo 
which, for the sixth-order polynomials [Eqs. (19) ], 


becomes 
C,VR, = 4A(A) VR[L + (A/10)] (35) 


where k is given by Eq. (29). Although C,VR; de- 

pends on % through the appearance of k in Eq. (35), the 

quantity C,R, depends only on A. When k in Eq. (35) 

is replaced by R,/R,’, there results 

C,R, = 4A(d) [1 + (A/10)] (36) 

Eqs. (29) and (35) may be compared with Eq. (26) 

when there is no pressure gradient so that 7%,’ = %," = 
4 =0. Inthis case one has, from Eqs. (21) and (25), 

(d/dx) (\/iy') = 36.56/%,(1 — 1%?) (37) 

in which the right side is a constant. Integrating this 

for \/m,’ and substituting in Eqs. (29) and (35) yields 


uy : fur — 
0 = C, = 0.667, 
Vx : VY} 


(D) Relation Between ui and M 


R, = 2.29 R,? = (88) 


It is convenient to be able to convert readily 
between #7 and J. These two quantities are related 
by 

M,? = 21," a 1) (1 aes iy") (39) 
which follows by substituting Eq. (7) and the perfect 


gas relationship 
G = [v/(y — JR (40) 
into the expression 


M,? = u?/yRT (41) 


(E) Computations of \ for the Airfoil 


The results tabulated in Tables 1 and 2 for flow over 
a thin biconvex airfoil depend on the determination of 
\(s), the proper integral of Eq. (21). The computation 


t See reference 9 for details. 
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Fic. 5. Isoclines and integral curve for integration of Eq. (21). 
of the coefficients (22) in this equation was based on an 
approximate formula for u derived by Tollmien, '* 


ty = U. (1 + ce + Ge’) (42) 


where ¢ is the slope of the airfoil with respect to the 
direction of flight and c,; and c2 are the following func- 
tions of Mach Number: 


VM.2—1 

= \ 1) Met +2] (43) 
. = —~1)M.4+2 
Oss } 


qQ = —] 


The boundary condition that \ be zero at the leading 
edge was used, and the equation was integrated by the 
method of isoclines. The isocline field and integral 


curve are shown in Fig. 5. 
(IV) STaBILiry; COMPUTATION OF (Ro)ecr. min. 
This computation is done by applying formulas de- 
rived in reference 2. In the notation of this report, 
these formulas are 


6 T (co) | 
(Ro) er. min. ™ - ‘a * (44) 
Tw Cot V | M,*(1 —_ Co)? 
where Cp is the value for c for which 
m{l1 — 22 (c)|]V(c) = 0.580 (45) 
Here, 
dw (T)? 1 (dw/d 
dn I/w U(dw/dn)* dn vi —_? 
= (c) = (n./c) (dw/dn)w — | (47) 


and ¢c is a number, 0 < ¢c < 1. In Eq. (47), n, is de- 
fined by the relation c = w (n,). 
On the basis of Eq. (7) 


be written as 


the temperature ratio 7 may 


? 
’ 


T = = Tw = (48) 
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fcr dr fi g? dt X 10° 
0 0 


TABLE 3 


5] 


Site dt X 108 4f dg I as 

0 dr dr dr? dr? 
0 2.00 0.200 0 —1.00 
0.108 1.98 0.104 —(0.486 —O0.899 
1.72 1.67 —0.028 —2.65 — (0.387 
$04 1.00 —0.056 —3.75 0.0625 
5.23 0.206 —0.017 —2.10 0.1600 
5.41 0 0 0 0 


[f and g are the sixth-order polynomials, Eq. (19) |] 


T ¥ g 
O 0 0 0 0 
0.1 0.200 0.0151 0.001383 0.00878 
0.3 0.573 0.0209 0.03438 0.0911 
0.5 0.844 0.0109 0.139 0.147 
0.75 0.986 0.0013 0.357 0.151 
1.0 1.00 0 0.605 0.153 
Thus, 
V(c) == (dw dn) w +4 


-. | — 9” + ona, ()'] (49) 
U(dw/dn)* L° dn? dn Bea 


Since our information on the steady-state flow is 
given in terms of 7, it is necessary to express dw/dn and 
d’w/dyn*? in terms of derivatives with respect to r. 
Based on Eqs. (11) and (12), one finds 

dw ou 6 (1 — &,*) 


= . = (50) 
dn dr t(6) a, (1 — a?) 


and 
d*w ie EB 19 (*) tl | x 
dy? dr? Or/ 1 — # 


| 6 | Q—a7/c- 
—— (51) 
t( 6) uy(1 — a7)? 


To find 6/t(6). the formula 


6 l 
2 - 9 ( 1) x 
t(6) (1 — #2)" 


1 
[ — at f (f? + 2fg + are’) | ar (52) 
0 


is valid where the integration is carried out at a fixed 
5 and the Eqs. (19) are used for f and g. 

Since, for fixed 5, w is known as a function of 7, in 
order to find n, it is necessary to know 7 as a function 
of r. This may be found by 

t(6) 


- : x 
. 6 1 — a27o- 


E _ at f (f? + 2dfg + re) | dr (53) 
0 


For computational purposes one notes that a con- 
siderable part of the computations involved in Eqs. (45) 
through (53) is independent of the value of 5. Conse- 
quently, the numerical values of f and g and the neces- 
sary integrals and derivatives are tabulated for vari- 
ous 7 values in Table 3. These functions are given 


graphically in reference 9, where, in addition, a com 
puting procedure is given. 
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The Performance of Axial-Flow 
Compressors As Affected by Single-Stage 
Characteristics 


SEYMOUR M. BOGDONOFF* 


Princeton Unwersity 


SUMMARY 


A series of calculations were carried out to obtain the per 
formance of a six-, an eight-, and a 13-stage compressor, each 
designed for the same tip speed, to handle the same mass flow, 
and to develop a pressure ratio of 3.2. The simplified analysis 
was based on experimental cascade and single-stage studies. 
The effects of boundary-layer growth and distortion of the 
velocity profiles, which occur at the rear of multistage axial-flow 


compressors, were negletted. 


(1) INTRODUCTION 


a SIMPLE BLADE ELEMENT THEORY of axial-flow 
compressor design has been widely used, both in 
actual compressor design and in analysis, to determine 
the effects of many design parameters. The theory, 
however, gives no indication of the physical limitations 
of the calculations. The maximum pressure ratio per 
stage, maximum turning angles, or Mach Numbers 
must be determined from experimental tests. From 
many tests of multistage compressors, a “‘rule of 
thumb’’ has appeared which states that, for highest 
efficiency, widest range, and best stalling characteris 
tics, the number of stages should be large i.e., the 
pressure ratio per stage should be small. It is im 
portant to note, however, that most of these tests were 
conducted on machines designed without the benefit 
of many recent cascade and single-stage studies. It 
has been shown conclusively in recent rotor tests'~® 
that the efficiency of properly designed blades can be 
kept high even with extremely high loadings. This 
permits the design of efficient compressors that can 
deliver a desired pressure ratio in fewer stages than was 
previously possible. 

It is the purpose of this paper to examine the effect 
of single-stage characteristics on the overall per- 
formance of multistage axial-flow compressors. The 
single-stage variables to be studied are the blade 
loading or pressure rise per stage and the efficiency 
variation with mass flow. The overall multistage 
performance will be examined particularly with ref- 
erence to the range of operation with good efficiency 
and the stalling characteristics. The analysis is, of 
necessity, a simplified one. Since information is not 
available to evaluate the effects at the rear of a multi- 





Received November 22, 1949. Revised and received October 
13, 1950. 
* Assistant Professor of Aeronautical Engineering. 


stage compressor, the compressors to be studied are 
assumed to be made up of a series of single stages. 
The performance of each of these stages can be cal- 
culated from theory and actual single-stage tests with 
a relatively high degree of accuracy. The analysis 
will, therefore, neglect the effects of multistaging con- 
nected with the trailing vortexes (distortion of the 
profiles downstream) and _ boundary-layer 
growth. The calculations are based on experimentally 
determined efficiency curves, stalling data, and two- 
dimensional cascade data. It is believed that the 
results obtained give a good qualitative picture; the 
actual quantitative results must wait for a more 
complete understanding of the many multistage effects. 


velocity 


(Il) GENERAL PROCEDURE OF CALCULATION 


The general procedure was as follows: Three com- 
pressors were designed with six, eight, and 13 stages to 
develop a pressure ratio of about 3.2 and all handling 
the same mass flow. The rotational speed at the mean 
blade height was taken as 684 ft. per sec. and the de- 
sign axial velocity as 450 ft. per sec. for all cases. The 
only variable was the loading of each stage——i.e., the 
pressure rise, lift coefficient, or turning angle. 
axial and rotational velocities are in the operating range 
of many commercial axial-flow compressors in use 
today and were chosen only to give numerical results. 
All stages in each design had the same vector diagram 
at the mean blade height. This vector diagram was 
i.e., the stage was of the 50 per cent 
The design 


These 


symmetrical 
reaction type at the mean blade height. 
stagger angles 8 (defined as the angle between the 
entering air and the perpendicular to the blade row) 
and the turning angles 6 (defined as the deflection of the 
air passing through the blade row) are shown in Table 1 
for the three compressors. The product of lift coeffi- 
cient and solidity, ¢C;, is also included in this table. 

These designs were completed to the point of actual 
blade-shape design and were found to be within the 
region covered by references 2 and 4. The mass flow 
TABLE | 


Design Conditions at the Mean Blade Height 


Compressor B, Deg. 6, Deg oCl 
6-stage 49.17 29.02 1.47 
8-stage 46.5 21.08 1.11 
13-stage 43.15 12.97 0.70 
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through each compressor was then varied by steps from 
the flow that would stall the compressor to the flow that 
would cause part of the compressor to act as a turbine. 
The turning angles were predicted from the angle of 
attack and blade section using cascade data.4.? For 
each stage, the efficiency was determined from the 
data in references 1 and 3. This procedure is neces- 
sarily a step-by-step one, since each blade row is 
treated separately. From the final blade row, the 
overall pressure ratio and efficiency were calculated. 
This procedure was repeated with different rotational 
speeds and a second single-stage efficiency curve. 


(III) ASSUMPTIONS AND APPROXIMATIONS 


The assumptions and approximations included in 
the calculations are indicated below. A _ detailed 
discussion of the effects of each of these on the final 
results is given later in the paper. 

(1) The mean blade height was held constant for all 
stages. The flow path then looked like that shown 
in the sketch below. 

















(2) The slope of the d#/da curve was taken to be a 
constant value of 0.9, where @ is the turning angle 
through the blade row and a is the angle of attack. 

(3) A single composite curve of efficiency vs. mass 
flow for a single stage was developed using twice the 
losses shown for a rotor in references 1 and 3, Fig. 1. 

(4) The effects of Reynolds Number and Mach 
Number were neglected. 

(5) No effect of the change of boundary-layer 
thickness on the walls was included. 

(6) Calculations of each blade row were based on 
conditions entering the row. 

(7) Constant entrance pressure and temperature 
(2,116 Ibs. per sq.ft. and 560°R.) were assumed. 


(IV) DETAILED CALCULATIONS 


From the original design of each compressor, the 
passage areas and blade angles were known (corres- 
ponding to design mass flow and performance). The 
axial velocity was constant at 450 ft. per sec. through 
each machine at design conditions. To find the effect 
of varying the mass flow, a new mass flow was chosen, 


1951 


and the axial velocity entering the guide vanes was 
changed accordingly. The turning angle did not 
change, since the entering flow is still axial, but the 
absolute velocity has changed. The pressure and 
temperature after the guide vane was calculated from 


To = To’ + (V2'2/2C,) 


and 
Po! /Po = [T2!/To]}!/7 -— 9" 
where 
T = temperature, °R. 
VY = absolute velocity, ft. per sec. 
P = pressure, Ibs. per sq.ft. 
C, = specific heat of air at constant pressure 
n efficiency (assumed 100 per cent for guide 


vanes), stage efficiency determined from 
single-stage data 
y = ratio of specific heats 


()o = stagnation conditions ahead of compressor 
()2 = flow after guide vane 
()’ = new flow caused by change in entrance mass 


flow 
The density can be computed from 
P = pRT 


where p = density and R = gas constant. With the 
passage area known from the original design and with 
the density calculated as above, the continuity equation 
can be used to calculate the axial velocity leaving the 
blade row. 

The absolute tangential velocity is known from the 
axial velocity and the turning angle, so the angle 
relative to the rotating blade row can be obtained. 
The new angle of attack of the rotor blade is then 
known, and the change in turning angle is calculated 
from d6/da = 0.9. The new leaving velocity can then 
be calculated, and the pressure and temperature can be 
found as above. The new density is found, and the 
axial velocity is calculated again from the known mass 
flow, density, and area. This procedure was repeated 
for each blade row. The overall pressure ratio was 
taken as the ratio of exit static pressure to inlet static 
pressure. The efficiency was calculated from » = 
AT »/ AT act, Where AT, was the value obtained 
adiabatically from the pressure ratio. 

These calculations were carried out at a constant 
blade radius (actually one-half way between the root 
and tip of the blade) and therefore are independent of 
the flow type through the compressor. The results 
are, of course, an approximation to the actual flow, 
since the flows at the root and tip sections of the blade 
are not directly calculated. They are included, 
however, to some extent by the use of an efficiency 
curve obtained from actual tests of the three- 
dimensional blade row designed for free-vortex flow. 
This efficiency curve must be modified for flow types 
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$04 TABLE 2 : 
S P Q p Range Analysis, Curve A 
oO 4 
2 85 fi Range, % Design 
° “a é Mass Flow 
“ 7 6-stage compressor 41 
za 78 J 8-stage compressor 28 
= ff 13-stage compressor 14 
z 707 
2 —— EFFICIENCY CURVE A 
_ ——EFFICIENCY CURVE B 
g machine with highly loaded stages. The eight-stage 
- ; , ee 
4.8, S machine has a range of good operating characteristics 
°o o . ‘ . . e 
al ° of twice that of the 13-stage machine, while the six- 
s stage machine has a range of almost three times that of 
« 4.04 ~ w « . “ lh 
- + F the 13-stage machine. (See Table 2.) 
Oo 
& 3.64 - - E 
e rhe range is defined as the per cent change in mass 
4 3.24 ra cae ° 
§ 4 flow from stalled conditions to the point where the 
o 284 5 pressure ratio dropped to 2.2. At this point for all 
a . . 
ao compressors, the last stage was operating as a turbine. 
For any deviation from the design point, the fewer the 
7 number of stages, the higher the efficiency and the 
us . . . lower the change in pressure ratio. This would, of 
-30 -20 -10 ° 10 


PERCENT DEVIATION FROM DESIGN MASS FLOW 


FIG. 4 EFFICIENCY AND PRESSURE RATIO 
CHARACTERISTICS OF THE THIRTEEN STAGE COMPRESSOR 


other than free vortex. In general, the results seem to 
be similar if the blade shapes are properly designed, 
as shown by a comparison of the results of references 
1 and 5. 

Stalling of the compressor was arbitrarily taken when 
one or more blade rows reached the stalled condition. 
This point was found to be a constant 27 per cent below 
design axial velocity if the design condition was taken 
to be the point of maximum efficiency instead of the 
point used in references 1 and 3. Therefore, with a 
design axial velocity of 450 ft. per sec., the compressor 
was considered stalled if the velocity entering any blade 
row dropped to 328.5 ft. per sec. 


(V) RESULTS AND DISCUSSION 


(A) Performance Curves and Range 


The overall performance of each compressor is 


summarized in Figs. 2-4. It is immediately apparent 
that a change in the number of stages has a strong 
effect on the shape of the pressure ratio and efficiency 
curves. Increasing the number of stages steepens the 
curves, although all three compressors have the same 
pressure ratio, about 3.2, and efficiency, 88.5 per cent, 
at the design conditions. The pressure ratio for the 
six-stage compressor varies between 2.2 and 3.35 in its 
operating range, while the 13-stage machine varies 
from 2.2 to almost 4.5. This steepening of the curves 
associated with increasing number of stages was found 
to be true also for staged centrifugal compressors 
(Pfau). 

The most striking point, however, is the large 
increase in operating mass-flow range shown by the 


course, only hold up to the case where further increases 
in blade loading gave decreased performance of the 
single-stage machine. The range shown for the six- 
stage machine is considerably more than that shown 
by any axial-flow compressor in operation today. It 
has been shown by Howell that the outlet air angle has 
a considerable effect on multistage performance, but the 
results obtained cannot be explained on this basis 
alone. It should be noted that the steeper charac- 
teristics, with the associated greater pressure ratio 
range, may be more desirable for some applications such 
as gas turbines. 
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FIG.5 AXIAL VELOCITY DISTRIBUTION THROUGH 
THE SIX STAGE COMPRESSOR 
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PERFORMANCE OF AXIAL-FLOW COMPRESSORS 


To show the cause of the change in shape of per- 
formance curves with the number of stages, the axial 
velocity distributions through each compressor for 
various mass flows are shown in Figs. 5-7 (all plotted 
to the same scale; M, designates the design mass 
flow). From these curves, it can be seen that, for a 
given change in mass flow, the change in axial velocity 
through the compressor increases with the number of 
stages. This effect is due to a combination of two 
factors: (1) the ‘‘reheat factor’ or the inherent higher 
efficiency of the more highly loaded stages and (2) a 
given change in mass flow gives a smaller percentage 
change in pressure ratio for the mcre highly loaded 
blade rows. 


The first of these, the ‘reheat factor,’’ can best be 
understood by comparing one stage of the six-stage 
machine to two stages of the 13-stage machine. Both 
of them will give about the same pressure ratio. From 
the rotor tests of references 1 and 3, the maximum 
efficiency for both sets of blade rows is the same. The 
one stage of the six-stage machine may have an effi- 
ciency of 90 per cent; therefore, the overall efficiency 


of the unit we are considering (one stage) is 90 per cent. 
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FIG.6 AXIAL VELOCITY DISTRIBUTION THROUGH 
THE EIGHT STAGE COMPRESSOR 
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FIG. 7 AXIAL VELOCITY DISTRIBUTION THROUGH THE 
THIRTEEN STAGE COMPRESSOR 
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The two-stage unit, however, is made up of two 90 
per cent efficient stages, and, since the second stage 
must operate with a higher inlet temperature, the 
overall efficiency of the two-stage unit is less than 
90 percent. This ‘‘reheat factor’’ is small and, in these 
calculations, seems to have a secondary effect. 


The important factor in the change of shape of the 
performance curves is the percentage change in pressure 
ratio or axial velocity which occurs for a given change in 
mass flow. If the three compressors are compared 
from this point of view, it is obvious that the machine 
with the most lightly loaded blades will be most critical. 
The same change in angle of attack gives a much 
greater change in pressure ratio (as a percentage) for 
the stage with a light loading than for one that is 
highly loaded. In other words, for a given change in 
angle of attack (caused by a change in mass flow) the 
blade row with the highest turning angle is least 
affected. As an example, compare the axial velocities 
entering the second stage of the six-stage machine and 
the third stage of the 13-stage machine when the mass 
flow is 0.85 of design. In the original designs, the 
pressures developed at these points were about equal, 
and the axial velocity was 450 ft. per sec. In the six- 
stage machine, the axial velocity is now about 5 per cent 
less than the entering velocity, while the axial velocity 
in the 13-stage machine has dropped to 7.8 per cent 
less than entering velocity. This corresponds to an 
efficiency of 83.1 per cent for the following stage of the 
six-stage unit and an efficiency of only 82 per cent for 
the following stage of the 13-stage unit. 


The performance of the six- and eight-stage com- 
pressors was also calculated for 0.6 and 0.8 design 
r.p.m., although the accuracy of the calculations 
decreases with decreasing rotational speed. The 
results of these calculations, combined with the per- 
formance curves at the design r.p.m., are shown in 
Figs. 8 and 9. The steeper characteristics of the 
eight-stage machine are still noticeable. At 0.8 de- 
sign r.p.m., the six-stage machine has a mass flow range 
of about 35 per cent more than the eight-stage machine 
(compared with 50 per cent more at the design r.p.m.). 
Below 0.8 design r.p.m., the conditions are so far from 
design that there is little effect of number of stage. 
It is interesting to note that the constant efficiency 
contours of Figs. 8 and 9 are of similar shape. There- 
fore, the mass-flow range of any machine depends 
primarily on the steepness of the pressure ratio curve; 
the flatter the curve, the greater the range. 


(B) The Effect of Single-Stage Efficiency Variation 


The entire analysis, so far, has been based on one 
single-stage efficiency curve (designated Curve A). 
This curve was obtained by assuming that a stage 
would have twice the losses found in the rotor tests of 
references 1] and 3. This assumption may be con- 
sidered conservative, since the losses in a stator are 
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not so large as those in a rotor. However, this effi. 
ciency variation with mass flow may be overly opti- 
mistic when applied to blading at the rear of a com- 
pressor where the flows are not uniform and the 
boundary-layer effects are large. This also may be the 
case at low Reynolds Numbers or high Mach Numbers, 
To study the effect on the material thus far presented 
of changing this efficiency curve, a second curve 
(designated Curve B) was derived. This curve had 
the same peak efficiency as the first curve but with 
losses twice that of Curve A off the design point. A 
third case, that of completely neglecting the efficiency 
variation (assuming a constant 90 per cent efficiency), 
was also studied. The three curves are presented in 
Fig. 1. 

Spot checks were carried out on the six- and 13-stage 
compressors, and a complete study was made of the 
eight-stage compressor with efficiency Curves B and 
the constant 90 per cent curve. The effect on the over- 
all performance caused by using Curve B is shown by 
the dashed lines in Figs. 2-4. The general effect is to 
flatten out the pressure curves and peak the efficiency 
curves, but the results presented in Section (A) are 
still valid. It is important to note that the greater 
the number of stages, the more sensitive the com- 
pressor is to changes in the efficiency curve. 

The effect of using Curve B on the internal flows is 
shown as dashed lines in Figs. 5-7. The axial velocity 
in all cases is increased at the exit of the compressor 
and in most cases causes a reversal in slope of the 
velocity curve. Since the axial velocity through the 
whole compressor is increased to some degree, stalling 
will occur at a lower mass flow than in the original 
case. The range is therefore changed and is compared 
in Table 3. There is still a significant increase in 
range with fewer stages, but the increase is not so large 
as shown in Section (A) using efficiency Curve A. 

The effect of neglecting the efficiency variation 
altogether means that the blade row is assumed to be 
able to operate over a wide range of angles of attack 
with no change in efficiency. The effect of such an 
assumption on the eight-stage compressor at 0.85 
design mass flow is compared in Table 4 with the 


TABLE 3 
Range Analysis, Curve B 
Range, % Design 
Mass Flow 


6-stage compressor 37 
8-stage compressor 29 
13-stage compressor 17 
TABLE 4 
The Ejight-Stage Compressor, 0.85 Design Mass Flow 


Design R.P.M. 


Pressure Ratio Efficiency 


Efficiency Curve A 3.6 80.6 
Efficiency Curve B 3.3 76.6 
88.2 


Constant 90 per cent efficiency 
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results obtained with efficiency Curves A and B. 
The effect on the axial velocity distribution is shown 
for the same cases in Fig. 10. 


(C) Stalling Characteristics 


In the previous sections of this paper, stalling of the 
compressor was taken as the point where the axial 
velocity entering any blade row dropped to the value 
that stalled the blade row in the tests of references 
1 and 3. This point may or may not be the point at 
which the whole compressor stalls. It is generally 
believed that a multistage compressor stalls at the rear 
first. This conclusion is reached from the following 
simple picture of the flow through the compressor 
when the axial velocity entering the first stage is below 
the design velocity. The blade angle of attack will 
then be higher than the design angle, and the first 
stage will develop a pressure ratio that is also higher 
than the design. Since the flow area is fixed and the 
density is higher than design, the axial velocity entering 
the next stage will be lower than that entering the 
first stage. This leads to a still higher angle of attack 
and pressure ratio for the second stage, and, as a 
consequence, conditions entering the following stage 
are still further from design conditions. If this crude 
picture is carried to the end, it is at once obvious that 
the conditions at the exit of the compressor are furthest 
from design conditions and so will stall first. 

If, however, this off-design flow is examined a little 
more carefully than above, the general results change 
completely. When a stage operates off its design 
point, the efficiency decreases and, to some degree, 
counteracts the increase in pressure rise caused by 
the increased angle of attack of the blades. This is 
clearly seen in the axial velocity profiles shown in 
Figs. 5, 6, 7, and 10. 

It is apparent that only for the 13-stage compressor, 
with efficiency Curve A, is it possible for the last stages 
to stall first. In this one case it seems probable that 
as the mass flow is decreased the last stage stalls and 
that, with further decreases, more and more stages at 
the rear will stall. At some point the increased re- 
sistance of the stalled stages will completely stall the 
compressor. However, if the actual single-stage effi- 
ciency curve is similar to Curve B, then the above no 
longer holds true, since the axial velocity is almost 
constant through the compressor. The stalling of one 
blade row will, in all probability, stall the entire 
compressor. 

This complete stalling will also hold true for the six- 
and eight-stage compressors with efficiency Curve A. 
With efficiency Curve B, the critical stages are definitely 
at the front of the compressor. The question that 
has yet to be answered is whether a blade row operating 
near its stall point will be stalled by the fluctuations 
in velocity caused by a preceding stalled blade row. 
Actually, this same question must also be answered 
for the case of a following stalled blade row. The 
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answers to these questions are beyond the scope of 
this paper and must be determined by actual experi- 
The answer will necessarily depend on how 
Below the 


ments. 
far from stall the blade row is operating. 
design r.p.m., for example, the axial velocity increases 
rapidly through the compressor, Fig. 11, and so it is 
likely that, although the first stage or two may be 
stalled, the axial velocities at the rear of the compressor 
will be high enough for the rear blade rows to operate 
unstalled. 

If it is desirable to have a ‘‘soft stall’’—that is, have 
the last stages stall first to give warning of approaching 
complete stall—it seems possible to design blading 
to give this effect. If the original designs are used as 
examples, blade sections are alike for each blade row of 
each compressor, since all vector diagrams are similar 
and symmetrical. From the flat axial velocity dis- 
tributions, all stages will stall almost simultaneously. 
If, however, the same blade section is not used for each 
row, some determination of the stalling is left with the 
designer. By using blades of higher camber than 
optimum at the front of the compressor and operating 
at lower angles of attack, the compressor will be forced 
to stall at the exit. The blade rows at the front will 
still be some distance from stall when conditions at the 
rear cause the unmodified blade rows at the rear to 
stall. 
ciency, but, considering the flat efficiency curve of a 
well designed blade row,' * it may not be excessive. 
Of course, the high mass-flow end of the performance 
curve will also be affected, since the blade rows at the 
front of the compressor will be operating further from 
their optimum point than those in the original design. 


This overcambering will cause a drop in effi- 


(D) Effects of Assumptions and Approximations 


(1) The choice of constant mean blade height was 
chosen only for ease of computation. The vector 
diagrams at the design conditions then become similar 
for all stages. This choice of geometry has no general 
effect on the conclusions reached. 

(2) The choice of a constant value of d0/da = 0.9 
was also based on ease of computation. However, 
the results of reference 4 show that this is a good average 
value in the operating range of blades with a solidity 
of about 1.0. The approximation becomes bad near 
the stall and turbine-action points. One check on 
the eight-stage compressor using an experimentally 
determined turning angle curve showed no serious 
deviations from the results presented. 

(3) The efficiency Curve A was derived from tests 
of four blades of various loadings of references 1 and 3. 
The inclusion of results with efficiency Curve B and a 
constant 90 per cent efficiency allows the evaluation of 
any single-stage efficiency curve by direct comparison. 
The effect of having a lower peak efficiency is to accent 
the sensitivity of the compressor with the more 


numerous stages. The axial velocity at the exit for 
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off-design conditions will always be higher for the 
lower stage efficiencies. 

(4) It was assumed that the blades were operating 
above a Reynolds Number of about 200,000 and far 
enough below the critical speed so that compressibility 
effects were not present. For the conditions studied, 
the Mach Number was considerably below critical 
conditions.’ The effects of both Reynolds Number 
and Mach Number can be evaluated, if single-stage 
tests are available, by simply substituting the corrected 
efficiency and turning-angle results. 

(5) No effect of the change in boundary-layer thick- 
ness on the walls was included because little informa- 
tion is available on the subject. These effects will be 
greatest at the rear of the machine, tending always to 
increase the axial velocities. There is no experimental 
basis for saying whether or not the effect will be greatest 
for the six- or 13-stage compressor. However, since 
both develop about the same pressure ratio and the 
13-stage machine will be considerably longer than the 
six-stage machine, it seems probable that the effect 
will be greatest on the compressor with the greatest 
number of stages. 

(6) Again, for ease of computation, performance of 
each blade row was based on the entering axial velocity. 
To be absolutely correct, the mean of the entering and 
leaving velocities should be used.* ® For all cases at 
the design r.p.m., below the design mass flow, the 
effect of using the entering velocity is negligible because 
the velocity changes are small. Above the design mass 
flow for the eight- and 13-stage compressors, the effects 
are noticeable but small, tending to steepen the curves 
(increase the exit velocity). The same effect is found 
to a much greater degree for the 0.6 and 0.8 design 
r.p.m. cases. Spot checks indicate that the general 
results are still valid, but for better accuracy (if needed) 
the calculations should be done with the mean velocities. 

(7) The entrance pressure and temperature were 
chosen to simulate normal sea-level operation. The 
choice of these factors has no effect on the general 
results, since all three basic designs operate at the same 
conditions. 

(8) The effect of the important multistage effects of 
boundary-layer growth, trailing vortexes, and distortion 
of the velocity profiles cannot be evaluated by any 
means known to the author except by experiments. It 
is believed that these effects will undoubtedly change 
the actual figures that were obtained, but the general 
results will still be valid. The choice of axial and 
rotational velocities used in the calculations also have 
no significant effects as long as the blades operate in a 
region where Reynolds Number and Mach Number 
effects can be neglected. 
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(VI) ConcLUSIONS 


A simplified performance analysis of multistage 
axial-flow compressors was used to compare the per- 
formance of a six-, an eight-, and a 13-stage compressor, 
-ach designed to handle the same mass flow and develop 
a pressure ratio of 3.2. 

The following conclusions were reached: 

(1) The compressor with the fewest stages (highest 
loaded stages) has the widest mass-flow range. The 
six-stage compressor has a range of mass flow of 41 
per cent of design mass flow; the eight-stage, 28 
per cent; and the 13-stage, only 14 per cent. 

(2) The compressor with the fewest stages has the 
flattest pressure ratio and efficiency variation with mass 
The more numerous the stages, the steeper the 
pressure ratio 


flow. 
characteristics—1.e., the greater the 
range. 

(3) The flat characteristics are due primarily to the 
fact that, for a given change in mass flow, the per- 
centage change in developed pressure ratio or axial 
velocity is smallest for highly loaded blades. 

(4) The calculations show that, for most configura- 
tions, stalling of a multistage compressor occurs 
throughout the compressor at the same instant. This 
is in contradiction to the usual belief that stalling occurs 
first at the rear and progresses forward. 

(5) A method is suggested for forcing the compressor 
to stall at the rear first. This may be a way to elimi- 
nate the sudden sharp stall usually associated with 


multistage compressors. 
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Coupled Free Vibrations of a Swept Wing 


A. E. ENGELBRECHT* 


North American Aviation, Inc. 


ABSTRACT 


This paper presents a method by which the natural frequen- 
cies of the coupled free vibrations of a swept wing can be approxi- 
mately determined 

The continuous elastic structure of the actual swept wing 
is approximated by a composite wing consisting of a finite num- 
ber of sections. Each section contains a concentrated mass 
located at the outboard end of the section and an inboard weight- 
less elastic bar of either straight or circular arc plan form. The 
choice as to whether the straight or circular are bar is used in any 
section is determined by whichever form more nearly coincides 
with the locus of elastic centers of the section. 

The usual assumptions of thin bar theory are made. 
other than inertial loads are neglected, and the elements of the 


Loads 


sections are assumed to be in harmonic motion. 

A set of six equations is developed for each section. Each 
equation relates a shear, bending moment, slope, deflection, 
torque, or angle of twist at the inboard end of a section to the 
shear, bending moment, slope, deflection, torque, and angle of 
The six relationships are then com- 
A trial 


frequency is then substituted into each of the section matrices, 


twist at the outboard end. 
bined into a single matrix equation for each section. 


and by successive matrix multiplication a single matrix for the 
entire wing semispan is obtained. If the trial frequency chosen 
is a natural frequency, the boundary conditions will be satisfied. 
If the boundary conditions are not satisfied, additional trial fre- 


quencies must be tried. Methods for determining proper fre- 
quency choices are shown, together with a numerical example in 


the paper. 
INTRODUCTION 


hy OF THE RECENT TRENDS in the aerodynamic 
design of high-speed airplanes is toward the use 
of the swept wing. This paper presents a method by 
which the natural frequencies of the coupled bending 
and torsional free vibrations of a swept wing can be 
approximately determined. The continuous elastic 
wing is replaced by a composite wing containing a finite 
number of masses connected by weightless supporting 


structure. 


(A) Methods of Solution 


The methods usually used to determine the natural 
frequencies and deflection patterns of composite wing 
replacements for actual straight wings fall into two 
broad classes. The one class contains the trial deflec- 
tion methods in which an initial arbitrary bending and 
torsional deflection pattern is assumed, and the other 
class contains the trial frequency methods in which an 
initial arbitrary frequency is assumed. 

Two methods of the first class that have been suc- 
cessfully used on straight wings are the method of nu- 
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merical integration'~* °® and the matrix iteration 
method.** '* The method of the second class follows 
the ideas developed by Holzer* * ' and employs an 
arbitrarily chosen frequency. Operationally, the meth- 
ods of these two classes are designed to converge the 
initial arbitrary choice of the deflection pattern in the 
trial deflection methods and frequency in the trial 
frequency method to its true value in the composite 
system. The method developed in this paper is an 
extension of the trial frequency method, which will 
determine the coupled free vibrations of a swept wing. 


(B) Procedure 


The composite wing used to approximate the actual 
swept wing consists of a finite number of sections where 
each section contains a concentrated mass at the out- 
board end and an inboard weightless bar of either 
straight or circular are plan form. 

A set of six equations relating the boundary condi- 
tions at one end of each section to the boundary con- 
ditions at the other end of the section is developed. 
These equations relate the shear, bending moment, 
slope, deflection, torque, or angle of twist at the in- 
board end of the section to the shear, bending moment, 
slope, deflection, torque, and angle of twist at the out- 
board end. The six relationships are then combined 
into a single matrix equation for each section. 

After the matrix equations have been written for each 
section of the wing semispan, a trial value of the fre- 
quency is selected, and its value is substituted into the 
section matrices. Then the matrix representing the 
tip section is premultiplied by the matrix representing 
the section just inboard of the tip. The product of 
these two is then premultiplied by the matrix repre- 
senting the next inboard section. This successive 
premultiplication is continued until a single matrix 
equation for the entire wing semispan is obtained. 

The matrix equation for the entire wing semispan is 
then reduced to a set of six equations. Each equation 
relates one of the above-mentioned boundary condi- 
tions at the center of the wing span to the boundary 
conditions at the tip. 

At the tip of the wing the shear, bending moment, 
and torque are equal to zero. At the midspan of the 
wing for symmetric free vibrations the shear, slope, 
and torque must equal zero. At the midspan of the 
wing for antisymmetric free vibrations the bending 
moment, deflection, and angle of twist must equal zero. 

If the trial value of the frequency chosen were a 
natural symmetric or antisymmetric frequency, the 
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6 = circular arc bar section central angle, rad. 
s = distance along elastic axis, measured from inboard end 
of section, in. 
elastic e = distance from elastic axis to center of gravity of section, 
in. 
h = vertical displacement of section center of gravity, in. 
L = length of bar between masses of straight bar section, in, 
w = frequency of coupled vibration, rad. per sec 
t = time, sec. 
g = acceleration of gravity, in. per sec.? 
fuselage ¢ K = frequency factor, w?/g, rad.? per in. 
n = subscript denoting section number 
m = subscript denoting mass 
b = subscript denoting bar 








Fic. 1. 


Vibrating system. 
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Fic. 2. Coordinate system. 














Dynamic free body of concentrated mass. 


Fic. 3. 


applicable boundary conditions at the tip and the center 
of the wing span would be satisfied. If the boundary 
conditions are not satisfied, a new trial frequency must 
be selected and the procedure repeated. <A procedure 
that also helps in choosing subsequent trial frequencies 
will be described in the next section. 


SYMBOLS 


W = weight of section of wing, lbs. 
I = bending moment of inertia of area of section of wing 
with respect to chordwise axis, normal to elastic axis, 


in.‘ 

E = bending modulus of elasticity of beam, lbs. per in.? 

G = shear modulus of elasticity of beam, lbs. per in.? 

J = polar moment of inertia of area of section of wing with 
respect to the elastic axis, in.* 

S = static unbalance about wing elastic axis, in.]bs. 

I = weight polar moment of inertia of wing section with re- 


spect to the elastic axis, lb.-in.? 
R = radius of curvature of bar section, in. 


DYNAMICAL ANALYSIS OF SWEPT WINGS 


Consider a continuous elastic swept aircraft wing 
whose section properties and mass distribution are not 
uniform. The continuous elastic system is approxi- 
mated by a system containing a finite number of masses 
connected by weightless elastic bars. The wing is 
thus divided into a finite number of sections, where 
each section contains a concentrated mass and its in- 
board weightless bar. A typical division is shown in 
Fig. 1. 

A weightless bar of either straight or circular arc 
plan form is used in any particular section, depending 
on which plan form most nearly duplicates the locus of 
elastic centers of the actual wing. Angular continuity 
is maintained between sections. 

An orthogonal coordinate system is used in the der- 
ivation of the section equations. The origin of the 
coordinate system is set at the inboard end of the sec- 
tion with positive directions as shown in Fig. 2. 

The oscillation of the swept wing is considered as a 
combination of bending motion in the z direction and a 
torsional motion in a plane normal to the centerline 
of the bar. No chordwise motion in the x direction is 
considered, since chordwise stiffness of an airplane 
wing is much greater than the spanwise bending or 
torsion. 

There are six boundary conditions inherent in this 
type of motion; these involve shear, bending moment, 
slope, deflection, torque, and angle of twist. A matrix 
equation relating these boundary conditions at each 
end of a vibrating section is developed. 


(A) Analysis of Circular Arc Section 


Consider the concentrated mass detached from the 
bar, Fig. 3. Six relations are obtained from the free 
body of the concentrated mass. Assume that the mass 
is in harmonic motion so that the vertical displacement 


h = hsin wt 


Then if one considers only the inertial loads, the sum 
mation of forces in the vertical direction reduces to (see 
Fig. 3) 

WO? (h sin wf) 


w = Vn 


m n 2 of 9 


s 














or 


or 


wl 


Al 


fr 
at 





ard end 


section, 


in. 
ion, in, 


wing 
e not 
IF OXi- 
lasses 
ng is 
vhere 
ts in- 
vn in 


r are 
iding 
us of 
luity 


der- 
the 
sec- 


aS a 
nda 
rline 
mn is 
lane 


y 
4 or 


this 
ent, 
trix 
ach 


the 
ree 
ass 
ent 


see 








i 











COUPLED FREE 


or 
oVn = mVan t+ (W/2)hw? 
Since 
h=z+e 
oVn = mVn + (w?/g) (We + WeB) 
or 
oVan = mVn + AWe + KSB (1) 

where 

K = frequency factor for coupled motion, w*/g 

S = static unbalance, We 


Also, from the free body, 


‘ foV.| 1 0 0 
MM, 0 1 0 
ze |=10 0 1 
bn 00 0 
oly 9 0 0 
frBn J} LO O O 











Next consider a circular are bar portion of a wing. 
The subscript ) has been omitted from Eqs. (8) through 
(33), because these equations are needed only in the 
derivation of Eq. (34). 

The shear, bending moment, and torsional relations 
for the bar portion of the section are found from the 
free body of the bar, shown in Fig. 4. These relations 


are 
V» 4 * Ve (8) 
M,-1 = V,R sin (0 — 6) + M, cos (@ — 0) — 
T, sin (0 — 6) (9) 
T,. = V,R[1 — cos (@ — 6)] + 
M,, sin (@ — 0) + T, cos (@ — 0) (10) 


The three remaining relationships necessary to re- 
late the boundary conditions at one end of the bar to 
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Fic. 4. Free body of bar. 
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M, = aM, (2) 

mn = alte” (3) 

bon _ mon (4) 

oBn = nie (5) 


and 
ol n = * = (W, geh 


By the same analysis as was used to obtain Eq. (1), 


this becomes 
ol = 


mln + KS(m2n) + KIp(mBn) (6) 


where J, = weight moment of inertia. 
Eqs. (1) through (6) are represented in matrix form 
by the following equation: 


KW 0 KS] [aa | 
0 0 O aM 
0 0 O nf 
] 0 0 aie (7) 
KS 1 KIg| | mTn 
oS © 135 Labs J 











those at the other end are found from the geometry of 
the bar when it is deflected and twisted from its orig- 
inal position. The usual thin bar assumptions are 
made. For small deflections the deformed shape of the 
bar is determined by the displacement of the section 
centroid and the rotation about the bar centerline. 
Under these conditions, the deflection and twist are 
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Fic. 6. Pure bending. 

considered separately, and the final change in curvature 
and twist is obtained by adding the individual com- 
ponents, !? 


Consider the bar restrained in such a manner that 


only an angular displacement 6 is obtained. The 
curvature, after twist, as shown in Fig. 5, is 
1/R, = — sinB/R & —B/R (11) 
and the change in angle of twist per unit length is 
@ = dB/ds (12) 


The bar is next considered restrained in such a 
manner that only a pure deflection is allowed. As pre- 
viously stated, no motion in the chordwise plane is 
considered. From this geometric consideration, see 
Fig. 6, two more relations involving R; and ¢ are ob- 
tained. 

1/R, = d*z/ds? (13) 
@ = —(1/R) (dz/ds) (14) 
since —dB = dz/R = — dds. 

Eqs. (8) through (14), together with the bar formu- 
las, 

EI(1 R;) = M,, a (15) 

GJ = Ty (16) 

are sufficient to determine the relationships between the 
boundary conditions at each end of the curved bar. 


The curvature for combined bending and torsion is 
obtained by adding the contributions of Eqs. (11) and 
(13). 


1/R, = (d*z/ds*) — (8/R) (17) 

The angle of twist is given by the sum of Eqs. (12) 
and (14). 

@ = —[(1/R) (dz/ds)|] + (dB/ds) (18) 


When the initial shape of the bar section is a circular 
arc, Eqs. (17) and (18) are written in terms of the cen- 
tral angle, ds = Rd6,. Then, 


1/R, = (1/R?) [(d?z/d6,”) — RB) (19) 
@ = (1/R*) [—(dz/d6,) + R (dB/d0,)] ~— (20) 
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Substitution of Eqs. (8), (9), (10), and (20) into Eqs, 
(15) and (16) will yield 
dz R? 


— RB = ~ 


VR in (A 6,;) + 
d0,2 nen : 


M,, cos (0 4;) T, sin (0 — 6)] (21) 


and 


dz dp oe R? 
dé, do, GJ 
M,, sin (0 


{RII - #,)] + 


cos (6 
6:) + T, cos (0 — 6)} (22) 


Integration of Eq. (22) with respect to 4; and solving 
for RB will give ; 


RB = 2+ (R?2/GJ) } V,R[0 + sin (6 — 0)] + 
M,, cos (06 — 6;) T, sin (0 — 6)} + C, (23) 


By substituting Eq. (23) into Eq. (21), the following 
differential equation is obtained: 


(d?z/d0,2) — z = V,R|(A + B) sin (6 — 6,) + Bé,] 4 
M,(A + B) cos (@ — 6) 
T,(A + B) sin (6 — 6) + C, (24) 


where A = R?/EI and B = R?/GJ. 
The solution of Eq. (24) is effected by the operator 
method or variation of parameters and yields the fol 


lowing expression for 2: 


z = V,R[—(1/2) (A + B) sin (@ — 6) — BO,] — 
(1/2) M,(A + B) cos (6 — 0) + 
(1/2)T,(A + B) sin (6 — 6) —- CG + 

Cz cosh 6, + C3 sinh @,; (25) 


Replacing 2 of Eq. (23) with the value of 2 given by 
Eq. (25) and reducing, the value of R@ is 
RB = —(1/2)V,R(A — B) sin (@ — &) — 


(1/2) M,(A — B) cos (@ — 6) + (1/2)T, X 
(A — B) sin (@ — 6) + C. cosh 6 + C3 sinh 6 (26) 
The expression for z’ is obtained by differentiation 
of Eq. (25) with respect tos. Since 
dz/ds = (1/R) (dz/d@,) 


the expression for 2’ is as follows: 
, |i 
sg’ = VJ, 5 (A + B) cos (8 — 4) — B — 


A+B\. _ (A+B 
M, (- OR ) sin (@—96,) —T, (= ) x 
C; # " 
cos (@ — 6;:) + - sinh @; + 7 cosh @; (27) 


The unknown constants C;, Co, and C; are evaluated 
in terms of the boundary conditions at the outboard 
end of the bar. 


(when 6; = 6) (28) 


/ 
= ¢ e= 2 -_ 
2 = Sn, os 6? B aa B,, 
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When these conditions are imposed on Eqs. (25), Ci = —V,RB@— M,B — 2, + 8,R ) 
MR) « 97 >» “sno relati are ni . | 
96), and (27), the following relations are obtained: C, = (1/2)V,R(A — B) sinh 6 + 


(1/2)M,(A — B) cosh 6 — 2,’R sinh @ — 
(1/2)T, (A + B) sinh 6+ 8,Rcosh 6 - (30) 


\ 


s = —V,RBO@ — (1/2)M,(A + B) —-GQ,+ 


C, cosh 6 + C3 sinh @ 
C; = —(1/2)V,R(A — B) cosh 6 — 


(1/2)M,(A — B) sinh 6 + z,’Rcosh 6+ | 


Re,’ = (1/2) V,R(A — B) — 


vs oe - r (29) 
(1/2) T,(A + B) + Cy sinh 6 + C; cosh 6 (1/2)T,(A + B) cosh 6 — 8,R sinh 6 | 
RB, = —(1/2)M,(A — B) + Cy cosh 6 sd | Substituting the values of the constants C;, C2, and C; 
Cs sinh 6 ) of Eqs. (30) back into Eqs. (25), (26), and (27) and at 
the same time setting 6, = 0, the following bar equa- 
Solving these equations for C), C2, and C3 gives tions are obtained: 


, an .. ; Je 
tri = Val (sinh @ — sin @) (sinh 6 + sin 6 — 20)] + 


(cosh @ + cos 6 -2)| + z,’ (—R sinh 0) + 2, + 


By B 
2 M,, E (cosh 6 — cos 0) — - 


a ew B 
r.| . (sinh 6 — sin @) — = (sinh @ — sin 0)| + B, | R(cosh @ — » | (31) 


wn 


5 A B 
a.’ = Fi = = (cosh 8 — cos 0) + - (cosh 6 + cos 6 — 2) | + 


A B 
M, | - OR (sinh 6 + sin 8) + OR (sinh 6 — sin 0) | + 2,’ (cosh 0) + 


~ 


nit B 
; E (cosh 6 — cos 0) + = (cosh @ — cos 0)| + 8, (—sinh @) (32) 


and 


pp .. ; 
; (sinh 6 — sin 6) | + 


=? oe ; 
Boi = Vy = (sinh 6 — sin @) 


a. B 
M, E (cosh 6 — cos #) — OR (cosh 6 — cos | + z,’(—sinh 6) + 


‘i | - OR (sinh 6 — sin #) — OR (sinh @ + sin 0)| + 8, (cosh 6) (33) 


The matrix equation [Eq. (34)], relating the boundary conditions at each end of the bar, is obtained by com- 
bining Eqs. (8), (9), (10), (31), (32), and (33). 

The matrix equation relating the boundary conditions at each end of the section containing the concentrated 
mass and its inboard circular are bar is obtained by replacing the right-hand column matrix of Eq. (34) with the 
right side of Eq. (7). This product is given by Eq. (35). See page 334 for Eqs. (34) and (35). 


(B) Reduction of Circular Arc Section to Straight Section 


The vibrating system, which has been selected to approximate the continuous structure of a swept wing, is repre- 
sented by sections containing either the circular arc bar and a mass or a straight bar anda mass. Eq. (35) relates 
the various boundary conditions at the outboard end of the section to the boundary conditions at the inboard end. 

The matrix equations relating the boundary conditions at the outboard end of a straight bar section to those 
conditions at the inboard end can be obtained by applying a limiting process to Eq. (35), the circular bar equa- 
tion. Eq. (35) is first rewritten by replacing R with R = L/6, where L is the length of the section measured along 
the centerline of the bar. This will reduce the equation to a function of @ only. Passing to the limit, as 6 — 0 
(R— o), the matrix equation of a circular are bar section becomes the equation for a straight bar section. The 
six relationships thus derived agree with those used by Targoff'! and Myklestad® for the determination of the 
vibration characteristics of a straight wing. 

As an example of the limit procedure, consider the part of the matrix equation [Eq. (35)] which relates »2,1’ 
with the boundary conditions at the outboard end of the section. Substitute R = L @ and expand the equation 
to obtain the relation found at the top of page 335: 
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To handle an actual concentrated mass, such as tip 
tanks, bomb racks, and nacelles, Eq. (39) is used with 
L=0. 

Eq. (35) or (39), depending on whether the section 
contains a circular arc or straight bar, relates the bound- 
ary conditions at the ends of a section. Next, a single 
matrix equation for the composite wing, made up of 
all the sections, will be developed to relate each bound- 
ary condition at the wing center span to the bound- 
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: L? (cosh @ — cos 6 L? (cosh @ + cos 6 — 2 
52 n “" =~ m J n _— " o o + — a ae ae + 
. E] 26° G. 26? 


L /sinh 6 + sin 4 L /sinh 6 — sin @ , 
m M, = ” 269 ations 6 + ee + m=n (cosh 6) + 


GJ 20 


cosh 6 — cos 6 L? /cosh 6 + cos 6 — 2 
rie ia 8 = nen = + 
20? GJ 26? 
KS | ‘fe (= 0 — cos *) 4 5, (= 6 — cos @ i} 4.7 L /cosh @ — cos ’) 
~ LEI 20 GJ 20 ) , 20 + 
cosh @ — cos 6 L? /cosh @ + cos 6 — 2 
eee a ee ane = + 
26? GJ 26? 


K i cosh @— cos @ rr. cosh @ — cos 6 ; mm 
EI ; 20 ae nae GJ 99 Ss — sinh Of (36) 








The following limits are necessary : 
. cosh 6 — cos 6 ] . cosh 6+ cos @d — 2 . sinh @ + sin @ 
lim —— —=-, lim —————_ ~ =0, lim ————— = 
6—>0 20? 2 0—>0 20? 0—>0 26 : 
. (37) 
. sinh 6 — sin 6 : 
lin. —-——————- = 0, _ lim [(cosh 6 — cos @)/20] = 0 
o—0 206 6—>0 
Substituting the values obtained for the limits in Eq. (37) into Eq. (36) yields 
i L KWL? KSL? 
2, = 7 Ve — <= M +" ad wil ce oe ee re) (38) 
2EI El 2EI 2EI *" 


The same procedure is used on the remaining portions of Eq. (35). Thus, the straight bar section matrix is 











- , | 0) 0 KW 0 KS L oe a 

»M 1 L | 0 KWL 0 Bw aM 
a | sphial saihs 

a L? | L KWL? ‘ KSL? = 

jie 2K] | KI 2EI 2K] sais 

= (39) 

3 4 @ , 3 —_ 3 

7 L L 2 (: =r) a KSL ia 

G/e] 2k] 6FI 6k 

a a 0 0) 0 KS l KI; Po 

, = () () 0 = KSL . ( _ *) = % 
GJ GJ GJ 

i JL JL | 


ary conditions at the tip. The composite wing is 
the replacement wing, composed of a finite number of 
sections, used to approximate the actual wing when 
the vibration analysis is made. 


(C) The Composite Wing 


Each section is assigned a section number starting 
with one (1) at the center of the wing span and continu- 
ing to m at the tip. The matrix equation, either Eq. 
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(35) or (39), which is used for the tip section is written 
symbolically as 
[C,-1] = [D,.] [Cr] (40) 


The equation form of the station just inboard of the tip 
is 


[C,-2] = [D1] [Cra] (41) 
Eq. (40) is set into Eq. (41) to give 
[C,-2] = [D,-1]. [D2] [C,,] (42) 


After ” such substitutions, 
[Co] = [D,] [D2] . 


A trial value of the frequency factor A is then chosen 
and is substituted into the matrix equations for the vari- 
ous sections {noted symbolically as [D,] to [D,,]} in Eq. 
(43). The indicated matrix multiplication [D;][D2]... 
[D,,] is carried out, and the resulting matrix equation is 
reduced to its set of six linear equations. Each equa- 
tion relates a shear, bending moment, slope, deflection, 
torque, or angle of twist at the center of the wing span 
to the shear, bending moment, slope, deflection, 
torque, and angle of twist at the tip. 

The boundary conditions at the wing tip are 


.. (Da) [Co] (43) 


V,= M,=T, =0 (44) 


Additional boundary conditions for symmetric vibra- 
tions are 


Vo = Zo’ == as = 0 (45) 
and for antisymmetric vibrations 
My = & = By = 0 (46) 


Set z, = 1, and for the symmetric case, for any particu- 
lar choice of K, solve for z,’ and 6, from the equa- 
tions 2’ = Oand 7) = 0. Substitute the values found 
for z,’ and 8, into the Vo) equation. The shear thus 
obtained is plotted against the frequency, and, where- 
ever the resulting curve crosses the zero shear axis, a 
natural symmetric frequency occurs. For the anti- 
symmetric case, set z, = 1, 8) = 2% = O and solve for 
My. Plot the values obtained for the bending moment 
against the frequency, and, wherever the resulting 
curve crosses the zero bending moment axis, a natural 
antisymmetric frequency occurs. 

Using this procedure, z, is always set equal to 1. 
If, for a particular value of K, z, is actually equal to 0 
for the symmetric or antisymmetric mode, the value of 
Vo or My will be driven to ~. Thus, for two values 
of K in this neighborhood, the curve may be drawn to 
cross the axis when actually an infinite discontinuity 
should occur. A little experience with curve shapes 
will indicate when such a possibility exists or an alter- 
nate procedure can be used. 

The alternate procedure is based on the existence 
theorem for the solution of a set of homogeneous linear 


1951 


equations. Set the tip boundary conditions (44) into 
Eqs. (45) and (46). Then for each case calculate the 
value of the determinant of coefficients, and plot these 
values against the frequency. Wherever the resulting 
curve crosses the zero determinant value axis, a natural 
frequency occurs. 


ILLUSTRATED EXAMPLE 


A preliminary sectioning of the wing is made. The 
total weight and center of gravity of each of these pre- 
liminary sections are determined. Also, the weight 
moment of inertia with respect to the elastic axis is 
determined for each of the preliminary sections. One- 
half the fuselage properties is assigned to a zero length 
section at the fuselage centerline and labeled Section 
1. The total weight and moment of inertia of the 
second preliminary section are then considered concen- 
trated at the spanwise location of the center of gravity 
of that section. The same thing is done in the remain- 
ing preliminary sections. The actual Section 2 used 
in the calculation is composed of the spanwise structure 
from the fuselage centerline to the location of the first 
outboard center of gravity. This process is continued 
to the tip, and the actual sections then used run from 
one mass just through the next outboard mass. 

The average bending area moment of inertia and 
torsional polar moment of inertia of the actual sections 
used are assigned to the weightless bars in the actual 
sections. The static unbalance S for each section is 
determined by the weight W and the distance e from 
the center of gravity to the assumed elastic axis. The 
radius of curvature R and the central angle @ of any 
circular arc bar section is determined by the plan form 
of the circular arc chosen to approximate the locus of 
elastic centers of the section. The length Z of any 
straight bar section is the distance along the elastic 
axis between the spanwise location of the mass of the 
section and the next inboard mass, measured positive 
from root to tip. Numerical values of the section 
properties of a typical wing are listed in Table 1. 

The matrix equations for Wing Sections 2, 3, and 4 
are obtained by substituting the properties listed in 
Table 1 into Eq. (35). The matrix equation for the 
zero length section (Section 1) and the equations for the 
straight bar sections (Sections 5, 6, 7, and 8) are ob- 
tained by using Eq. (39). In order to get a more uni- 
form range of numerical values, the following substt- 
tutions are made: 


z = 3(10—*) } 
zs’ = 2'(10~°) (47) 
B = B(10-*) \ 

The matrix equations for Sections 3 and 8, for K = 20, 


are shown in Table 2. 

When the matrix equation for each section has been 
evaluated, a trial value for the frequency factor K is 
substituted into the matrix equation for each section. 
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Basic Wing Data 
































| Ww I | J L e 
| Section | (Lbs.) (In.*) (In.*) (In.) In.) 
| 1 5,226 776.7 1,650 0 15 
2 503 558.3 1,342 
3 | 413 475.7 1,056 
4 | 350 | 378.7 809 
5 | 112 208.7 517 37 0 
6 68 97.1 | 198 37 0 
7 68 | 58.3 | 99 | 37 0 
8 52 34.0 44 41 8 
E = 10.3 X 108 lb./in.? 
G = 3.5 X 10 lb./in.* 
S | Ig R 0 
Section (Lb.-In.) | (Lb.-In.?) (In.) Rad.) 
1 78,390 64,112,000 eT 
2 } 2,012 110,336 110 0.28 
3 4,543 130,808 175 0.14 
4 —350 71,400 175 0.21 
5 ‘ 0 40,000 
6 | 0 23,800 
7 0 18,240 
8 416 7,740 
TABLE 2 
Matrix Equation for Section 3, K = 20 
TV." =. | 0 0 | 0.008260 | 0 0.090860 V3 
M2 24.4195 0.990220 0 0.189026 | —0.13954 | 1.85370 M; 
3,’ —0.061085 —0.004978 1.00982 | —0.000431 0.000814 | —0.14388 2,’ 
Zo = 0.503171 0.061085 | —24.581 1.00355 | —0.00669 1.74672 23 
T2 1.71150 0.139540 0 0.104108 | 0.99022 2.74608 T3 
| Bs J —(.000937 —0.000114 —0. 14046 —0.000612 —0.00665 0.992348 B; 
Matrix Equation for Section 8, K = 20 
Vr] a. 0 0 | 0.001040 | 0 | 0.008320 Vs 
M, 41 1 0 0.042640 | 0 0.341120 Ms 
Zz, | = —2.40148 —0.117141 | 1 —0.002498 | 0 —0.01998 33’ 
27 32.8202 2.40148 —41 1.03414 0 0.273060 Zs 
T; 0 0 0 0.008320 1 0. 154800 Ts 
| Br 0 0 0 —0.002216 —0. 26625 0.958780 Bs 
TABLE 3 
Wing Matrix for K = 20 
Vo . — 49.4739 0.716660 6.46987 Vs 
Mo —342.682 5.11064 : 22.2424 Ms 
Zz’ |=]. 4.27681 —0.064939 Bate —0.78055 Z,’ 
pA ‘ ; —399 .423 5.84252 al 42.9245 Zs 
To — 1,695.08 20.9030 i 1,261.78 Ts 
Bo —0.728831 0.007562 a 0.919914 Bs 


Then the matrix of Station 8 is premultiplied by the 
matrix of Station 7. This product is premultiplied 
by the matrix of Eq. (6), and the process is continued 
until the matrix for the entire wing semispan has been 
calculated. In this manner each boundary condition 
at the center of the wing span is related to the boundary 
conditions at the tip. The first, second, and fifth 
columns of the section products are not calculated, 
since Vs = Ms = Ts = 0 (boundary conditions at the 
tip). Table 3 shows the matrix of the wing for K = 


20. 


The matrix equation for the wing is then reduced to 
a set of six equations. The three equations used to 
calculate the unbalance for the symmetric case are 


Vo = —49.4739%s’ + 0.716662, + 6.46987 Bs 
3)’ = +4.2768123’ — 0.0649392; — 0.7805498Bs 


Ty) = —1,695.08%s’ + 20.90302, + 1,261.78Bs 


Set Zs = 1, 3)’ = Ty = O and solve for 2s’ and Bs. Sub- 
stitute these values into the Vo equation and find Vp = 
—0.047566. Had Vo = 0, the value K = 20 would 
have been the frequency factor of one of the natural 
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symmetric frequencies of coupled wing vibration. 
Since Vo # 0 for K = 20, more values of K must be 
chosen and the unbalance graphed as shown on Fig. 7. 
The natural symmetric frequencies occur where Vp = 
0. 

The conditions used to calculate the unbalance for 
the antisymmetric case for K = 20 are 


My = —342.68223’ + 5.110642, + 22.2424; 
% = —399.423%.’ + 5.842522 + 42.9245Bs 
By = —0.72883123’ + 0.0075622. + 0.919914 Bs 


Set 23 = 1, % = Bo = O and solve for 3’ and Bs. Sub- 
stitute these values into the Mp» equation and find My = 
0.044425. Plot the unbalanced Mo against the fre- 
quency factor K as on Fig. 7. 

The determination of the unbalanced Vp and Mo for 


K = 20 illustrates the procedure used for any particu- 
lar choice of K. No natural frequency was found for 
K = 20, since the boundary conditions were not satis- 


fied. Therefore, more trial frequencies must be tried 
in order to locate the natural frequencies. Plotting the 
values of the unbalanced V») and M) as soon as they are 
obtained will aid in making more intelligent choices 
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of future trial frequency factors. Fig. 7 shows that 
Vo = O at approximately K = 5.85, 57, and 124 and 
that My = O at approximately K = 21.5 and 122. 
Therefore, the symmetric natural frequencies occur at 
nearly 454, 1,418, and 2,090 cycles per min. and the 
antisymmetric frequencies occur at nearly 870 and 2,070 
cycles per min. By successive choices of K, the fre- 
quencies can be determined to any desired degree of 
accuracy. The deflection pattern or any of the other 
four properties can be determined directly from their 
equations. 

The values of the determinant of coefficients, de- 
scribed as an alternate procedure to determine the nat- 
ural frequencies, are also plotted on Fig. 7. 
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Analysis 


of ie te ich Plates 


of the Elastic and Plastic Stability 


by the Method of 


Split Rigidities—I’ 


P. P. BISLAARDt 
Cornell University 


SUMMARY 


The author gives a detailed but simple derivation of his for- 
mula of the critical thrust for antisymmetrical buckling of sand- 


wich plates 


Pom tee Ci ry 


in which Py is the joint critical thrust of the single faces (Case 
0), P; is the critical thrust without taking account of the shear 
deformation of the core (Case 1), and P» obtains if only shear 
deformations of the core would occur (Case 2). In computing 
P, and P», the flexural rigidity of the faces is not taken into ac- 
count. 

The method is applied to twelve different cases of loading and 
boundary conditions, the coefficients for 
buckling in the plastic domain are also given and presented in 
graphs. A detailed explanation is included showing how the re- 
sulting formulas must be used to find the critical load P,,. The 
critical load P.rw for wrinkling of the faces is also considered. 


for which reduction 


INTRODUCTION 


— RIGOROUS DETERMINATION of buckling loads of 
sandwich plates is extremely difficult because of 
the involved interaction of the core and the face plates. 
An approximate solution of satisfactory accuracy can be 
obtained by subdividing the complex interaction into 
its component parts and computing for each of them 
their individual critical loads. An interaction formula 
that combines these individual critical loads then gives 
the buckling load of the assembly. 

The basic principle of this approach is first illustrated 
on a simple example. A simply supported homo- 
geneous bar (Fig. la) buckles in a half sine wave, 
whereby in the elastic domain in any cross section the 
internal moment J; is equal to the external moment 
M, = P,w, in which P, is the Euler load and w is the 
deflection of the bar at the cross section involved. The 
internal moment 7; depends only on the curvature of 
the bar, so that with bending in a half sine wave, as it 
does with buckling, the inner moment will always be 
equal to P,w, independent of the fact whether the bar 


Presented at the Structures Session, Eighteenth Annual Meet- 
ing, I.A.S., New York, January 23-26, 1950. Revised and re- 
ceived December 15, 1950. 

* On request of the Editor, the author split his original paper 
in two parts, of which this first part gives the derivation of the 
general formula and the necessary information for the structural 
analyst to use it. The second part will contain those theoretical 
derivations that are not available elsewhere. 
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is subjected to an axial force or not; hence, 
M, = Pew (1) 
As a simple application of this conclusion, the re- 


sulting deflection w of a bar, with an initial inelastic 
deflection w,, after application of an axial load P, will 
be calculated (Fig. 1b). Assuming the initial deflec- 
tion w; as a half sine wave, the thrust P will cause an 
Since, before P was applied, the 
internal the moment is 
caused by the deflection w, alone. 
half sine wave, it will cause an internal moment V7; 
P,w,, which likewise varies as a half sine wave. The 
external moment 


M, = P(w; + w,) 


excess deflection w,. 
moment was zero, internal 


If w, also follows a 


= Pw 


also will vary according to a half sine wave, so that at 
any point equilibrium exists if! 


M, = M, 
or 
P,w, = Pw (2) 
or 
(Px/P)w, = (Pz/P) (w — 
or 
w= \(Pe/P)/((Pe 1]}w, (3) 
DERIVATION OF THE GENERAL EQUATION FOR THE 


CRITICAL THRUST OF SANDWICH PLATES 


The same principle is now applied to the computa- 
tion of the critical thrust of sandwich plates, consisting 
for example of metal faces and a core of light material 
of small rigidity, of which Fig. 2 a cross section. 
First, the cylindrical buckling of a plate with free un- 
loaded edges is considered (Fig. 3). 
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external and internal moments for a 
sandwich column. 


Fic. 3. Comparing 


With buckling of homogeneous plates, only deforma- 
tion by bending needs to be taken into account. With 
sandwich plates, however, the modulus of rigidity of 
the core is relatively so small that deformations by 
shear in the core must also be taken into account. The 
shear deformation of the faces may be neglected, since 
its effect is negligible with respect to that of the shear 
in the core. Moreover, it follows from a more exact 
calculation* * that the strains of the core in a direction 
perpendicular to the middle plane may be neglected in 
practical cases of plate buckling (half wave length- 
thickness ratio larger than, say, 5). This does not 
apply, however, to wrinkling of the faces, which will be 
considered later. 

The rigidity of the structure is now split into two 
parts in such a way that the following two cases are 
considered: In Case 1, the sandwich plate has its actual 
bending rigidity, but the core is assumed to be infinitely 
rigid against shear. In Case 2, the core is assumed to 
have its actual shearing rigidity, but the sandwich 
plate is assumed to be infinitely rigid against bending. 
As, however, in the latter case, the shear deformation 
of the core would not be possible if the proper rigidity 
of the faces, being that against bending with respect 
to their proper middle plane, also would be infinite, 
this proper bending rigidity of the faces is assumed to 
be zero, so that the deformations according to Cases 1 
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and 2 may proceed independently of each other. The 
effect of the bending rigidity of the faces is taken into 
account separately afterward in Case 0. The critical 
loads per unit width for Cases 0, 1, and 2 are denoted 
by Po, Pi, and P2, respectively. 

Accordingly, first a reduced sandwich plate—i.e., 
with faces of zero proper rigidity—is considered. If it 
buckles, its total deflection w will consist of a deflec- 
tion w, from deformation according to Case 1 and a de- 
flection w2 from Case 2 (Fig. 3a). Its reduced load P, 
will be smaller than the real critical load P,, of the com- 
plete sandwich. 

Now it is assumed that only a deformation from 
Case 1 occurs, in which case the critical load will be 
P, (Fig. 3b). Since only a deflection w; occurs, the 
condition that the internal and external moments must 
be equal leads to 


or 
M, = Py, (4) 


according to Eq. (1). 

Assume, subsequently, that, while the internal mo- 
ment and shear remain the same, the rigidity of the 
core decreases from infinite to its real value. Conse- 
quently, the core is now deformed by the shearing 
stresses. This will increase the deflection by an 
amount we (dotted in Fig. 3b), so that it acquires its 
final value w = w,; + Ww. Since the internal moment 
remains constant and thus equals P,w,, according to 
Eq. (4), the compressive force required to maintain 
the equilibrium will decrease now to its reduced value 
P,. It has to satisfy the condition of equality of in- 
ternal and external moments 


M, = M, 
or 
Pw = 9 0 (5) 
yielding 
w, = (P,/P:) w (6) 


If the deflections w; and w» and, hence, w, and w have 
the same form, Eqs. (5) and (6) will be valid for any 
point of the plate, so that indeed at any point the 
equilibrium is satisfied. Since in any section the shear 
then will be proportional to P and w, according to Eq. 
(5) the shear remained constant during the extra de- 
flection w,.. Comparing Eqs. (2) and (5), it is clear 
that the deflection w, here plays the same role as the 
initial deflection w; in case of the bar, Fig. 1, since both 
w. and w;, éncrease the external but not the internal 
moment. 

It is next assumed that only a deformation according 
to Case 2 occurs (Fig. 3c) and, hence, that the faces 
are infinitely rigid against tension or compression and 
that only a shear deformation of the core takes place. 
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STABILITY OF 


Since the critical thrust is P2, it follows from /; = 
MV, that 


M,; _ 5Ws (7 ) 


Now let, with constant internal moment and shear, 
the rigidity of the faces decrease to its actual value, by 
which the deflection increases by an amount w,; (dotted 


in Fig. 3c). Inasimilar way, it follows then from 1/7; = 


M, that 
Powe = P.w (Ss) 
yielding 


Ws. = (?, Ps) w (9) 


As the plate is here again in its reduced state, the 
value P, is the same as in Eqs. (5) and (6). If in Figs. 
3b and 3c the deflections w, and w» of the first step are 
assumed to have the same ratio as in the buckled re- 
duced plate in Fig. 3a, then the values w in Eqs. (6) 
and (9) are equal, and w; + we = w. Hence, it fol- 
lows from Eqs. (6) and (9) that 


[((1/P1) + (1/P2) |P,w 


W = W, + We = 


or 
FP, = (Py* + Pe") (10) 


In order to take account of the proper rigidity of the 
faces, it is noted that, if their joint critical load is de- 
noted by Py (Case 0), the extra compressive force, in 
order to give them the same deflection w as the reduced 
plate, is equal to Py (Fig. 3d). Hence, it follows that 
the actual critical load of the complete sandwich plate 
is, according to Eq. (10), 


P. = Pot P, = Pot (Pi + Peo')-! (11) 


If at the critical load the stress in the faces is in the 
plastic region, this has to be taken into account by 
multiplying Py) and P; by a reduction factor 7. For 
this purpose the safe assumption is made that the en- 
tire force P,, is carried by the faces, so that 7 refers to 
the buckling stress o,, = P,,/2h. It can be found by 
trial and error from the appropriate diagram (Fig. 4), 
as will be explained in an example at the end of this 
paper. 

Consider, furthermore, a sandwich plate subjected 
to compression and supported at one or both of the un- 
loaded edges. Here, it is difficult to compare internal 
and external moments, and it is more convenient to 
compare the restraining forces R and the deflecting 
forces D acting on an element // dx dy (Fig. 5a). As- 
sume again a reduced plate of which the rigidity is split 
into Cases 1 and 2. Assume first a deformation ac- 
cording to Case 1 alone, so that the critical thrust is 
P,. The deflecting force D,, caused by the transverse 
components of the compressive force P), is held in 
equilibrium by the restraining force Rj, given by the 
transverse shearing stresses. At a certain point and 
with a certain form of the middle plane, the deflecting 
force D, is proportional to P; and to w;, and, hence, it 


may be denoted by cP;w;. Consequently, R = D or 
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Comparing deflecting and restraining forces for an 
element of a sandwich plate 


Fic. 5 


R, = cP \w (12 


Again, let the shearing rigidity of the core decrease to 
its real value, by which the deflection increases by an 
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STABILITY OF 
amount w2, while the internal stresses on the element 
remain constant and, hence, also R;. Assuming the 
deflection surface w. to have the same form as w, and, 
hence, also as w = Ww, + Wp, it follows that the reduced 
critical thrust P, now must satisfy, for any point of the 
plate, the condition R = Dor 

cP\w, = cP,w (13) 
so that 
(14) 


w, = (P,/P;)w 
being identical to Eq. (6). 
Assuming, finally, Case 2 to occur (Fig. 5b), the 
critical thrust will be P,, so that with the deflection w» 
itfollows that R = D or 
R. = cPowr 


(15) 


Keeping the internal stresses on the element and, 
hence, R, constant and letting the rigidity of the faces 
decrease to its real value, the deflection increases by an 


amount w,. Itisnow R = Dor 


CP owe = C 7 Ww (16) 
so that 
w. = (P,/P.)w (17) 


Eqs. (14) and (17) are identical to Eqs. (6) and (9), so 
that they also lead to Eq. (10) and, finally, to Eq. (11). 

It will be clear that Eq. (11) also obtains if the plate is 
subjected to compression in two directions and to shear 
or a combination of these forces. With shear forces 
P,, alone, with a certain form of the deflection surface 
w, the deflecting force D for a certain element is pro- 
portioned to P,,, and to w, so that it can be expressed 
as CP,,w, which leads in the same way to Eq. (11). 
If at the same time compressive forces are present with 
a fixed ratio of these forces to P,,, again the force D is 
proportional to P,, and w, so that Eq. (11) also obtains. 

In this derivation it was assumed that the deflection 
surfaces had the same form, so that the ratios w,/w and 
w/w were the same for any point. This obtains in- 
deed for Cases A,, Ao, J, and J and for Case B of Table 
lif Po, P:, and P2 are computed for the same half wave 
length in the X-direction. In the latter case, this half 
wave length has to be chosen such as to make P,, a 
minimum. Hence, for these cases, Eq. (11) yields ac- 
curate results. 

In other cases, the surfaces w), we, and w have differ- 
ent forms, but, since Py, P;, and P» are always calculated 
for the same half wave length, the difference is never 
extravagant. Since, in addition, admission of diverging 
deflection surfaces means a lack of Eq. 
(11) will yield conservative results for these cases. 


restraint, 


Indeed, by allowing the plate to choose for Cases | 
and 2 that deflection surface which makes P; and P, 
a minimum, it is provided with more choice than it 
actually has, so that the computed critical thrust from 
Eq. (11) may be smaller but never larger than the actual 
one. Seide* computed the buckling strength of sand- 
wich plates with simply supported loaded and clamped 


SANDWICH 


4 
Ge 


PLATES 


unloaded edges, assuming the faces to act as membranes. 
The author has shown that, after taking account of the 
thickness and proper flexural rigidity of the faces, it 
follows from his results that, even for this very unfa- 
vorable case where the shapes of w; and w» differ consid- 
erably, Eq. (11) yields results that are only from 0 to 4 
per cent conservative.’ For pure shear,®~* although 
w, and w». have different forms, the discrepancies seem 
to be smaller than | per cent and not more than those 
obtaining for more laborious methods. 

Ihe author gives the above rather lengthy proof of 
his formula because it appears, from papers that refer 
to his method, that the more condensed version of it 
given previously® * ' was not fully understood. 

P., must be calculated from Eq. (11) for several half 
wave lengths Z in order to find its minimum value, so 
that it is necessary to assemble data for Po, P;, and P» 
as functions of 8 = L/b for several cases of loading and 
boundary conditions, as for example for Cases A—A in 
Table 1. 


THe CriticAL Loaps Py AND P; FOR CASES 0 AND 1 


For all cases represented in Table 1, the critical loads 


Py and P, may be expressed as 
P» = nkon*D , b? 


P, = nkor’D, b? 


(18) 
(19) 


in which D; is the joint flexural rigidity of the faces, 
being (Fig. 2) 


D,; = E,h®/6(1 — v?) (20) 


where /, and v are the elastic modulus and Poisson’s 
ratio of the faces. JD, is the flexural rigidity of the 


plate as a whole, so that 
D, = E,-h(t + h)?/201 — v?) 


It should have an additional term /,f°/12(1 — »,?), 


in which /, and yp, refer to the core, to take care of the 
flexural rigidity of the core, but this term is so small 
that it nearly always may be neglected. Values / and 
t are the thicknesses of faces and core, respectively 
Value 6 is the width of the plates, as indi- 
Except for Cases A.1, J, and J, the 
n is the 


(Fig. 2). 
cated in Table 1. 
plates are assumed to have infinite lengths. 
reduction coefficient for plasticity. 


Values k 


Values ky in Eqs. (18) and (19) depend on loading 
and boundary conditions and follow from the theory of 
elastic stability of homogeneous plates.''~'* The per- 
tinent values are given in Table 1 as functions of the 
ratio 8 between half wave length L and width b. They 
are given graphically in Fig. 6 for Cases G and //. 
For Cases A, J, and J/, the load P: is independent of the 
form of the deflection surface, so that the plate will 
buckle here in that form, which makes Py and P; mini- 


mum. Hence, for these cases, only minimum values 
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Fic. 6. Coefficients ko for Cases G and Has functions of 8 = L/b. 


of Py and P; are needed, so that here only minimum 
values of ko are given. More details may be found in 
reference L5. 

It is known" that, if minimum buckling stresses are 
compared, the interaction curve for shear stresses r,, 
and longitudinal normal stresses o,, as occur in Case K 
of Table 1, is the parabola 


(Try Tryer)* + a, ‘Cue = |] (22) 


in which 7,,,, and o,,, are the critical stresses if r,, and 
o, work alone. The author!’ found that the 
relation holds if all critical stresses refer to the same 


same 


wave length, which case has to be considered here. 
This is plausible, because it may be shown that this 
parabolic relation means that the so-called reduced 
principal compressive stress governs the buckling. 

If for a certain half wave length Z in X-direction the 
critical stress for pure compression (Case B) is o,,, and 
that for pure shea: (Case G) is 7,,;,, it follows that, 
since for pure shear the principal compressive stress is 
also Trycr, in shear a principal compression is (oz¢, 


Nc = nH = 0.91 


The subscript (min.) indicates that a has to be chosen 
such that both numerator and denominator obtain 
minimum values. Furthermore, it was shown” that 


(a-'{[A(1 + a) + 2(B + 2F)a*] B~* + 23Aa? + B+ 2F) + AB*})min. 
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Tryer) AS dangerous as in pure compression. If in com. 
bined loading by r,, and o, the shearing stresses are 
accordingly replaced by what will be called their re. 
duced value ty, = (Grer/Tryer)Tzy, also in shear buckling 
occurs at a reduced principal stress 


Pir = Tryr (Orer Trycr)Trycr = Orer 


Hence, it may be expected that also in combined shear 
(ty) and compression (¢,) buckling occurs if the s0o- 
called reduced principal compressive stress is equal to 


Orer OF if ‘ 

é 9\2 o1'/2 
Pir = (Or 2) + [(c, 2)? + Tese | _ Orcr 
or 
2 2 2 

Or Orer re Or 
= a. = Orer — 9 
4 : 2 


The latter equation is identical to Eq. (22). Denoting 


for Case K 
— > > _— )* 
Oz/Try = Fx re = 7 (23 


Eq. (22) yields for Case K 


Tryer Y Trycr Y . 
Trycr \- <i 9 + [ + (: x) | t (24) 


In reference 15, the author checked this equation by 
caleulating o,,, (Case B), Tryp (Case G), and (ry). 
The last value, referring to Case A, was calculated for 


( Try )k 


plastically equivalent stresses o, and 7,, (i.€., ¢, = Try X 

ls _ /z ies —s 9 — 
V3ory = V3). The accuracy of Eq. (24) was found 
to be satisfactory. From Eq. (24), ko for Case K may 
be expressed in the ky values for Cases B and G by the 
formula 


(ko)e Y j | ryt ) es 
(ko , = ko G = . l ‘ (2: } 
pilin, ( (eon 2 VT Likode 24S ' 


Hence, for Cases G, H/, and K, the critical loads Py 
and P, in Eqs. (18) and (19) refer to the shear forces 
P,, per unit length of the plate. 


The Reduction Coefficients for Plasticity 


Values 7 in Eqs. (18) and (19) represent the ratios 
between the plastic and elastic buckling stresses for 
the given half wave lengths LZ or ratios 8 = L/b. They 
follow from the author’s theory of the plastic stability 
of plates and shells.'7~'® He derived the coefficients 7 
For Cases A—F, J, and J, values 7 
are given in Table 1, while for Cases G and /7,"® 


in reference 15. 


(26) 


fa 1 (1 + q@?)28 2 + 2(3a? +1)+ B*] btn. 


for combined compression o, and shear 1,,, Case A, 
nx may be expressed in ng and ng for Cases B and G for 


the same ratio 8 = L/h and the same equivalent stress 
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Plastic reduction coefficients 7 for Cases .4;, Ao, J and J 





Fic. 7. 
as functions of the equivalent buckling stress 
op = (or? — oroy + oy*)'/? 
for several ratios B = L/bd. 


e e\ 1/ - 
Gp = (o,? + 3rzy*) Hence, nz refers to 0, = og and 


Ng tO Try = Op ‘n/3. 
nk = (ve + 3na)/(y* + 3) 


Also nx is given in Table 


nx is given by the relation 
(27) 


where y is given by Eq. (23). 
i. 

For Cases A—F, J, and J, where the coordinate axes 
coincide with the principal axes of the state of stress, 
values A, B, D, and F in the formulas for 7 are given 
by Eqs. (21) and (2la) of reference 19. For Cases 
A-F in these Eqs. (21) and (21a), the ratio 8 = p2/p; of 
the principal stresses (which has nothing to do with the 
ratio 8 = L/b in the present paper) has to be equated 


tozero. This yields 
A [1 + 301 + e)F,/E]/y 
B = 2[1 — (1 — 2v)E,/E]/y¥ 
D=4/y 
where 
y = 5 — 4v+ 3e — (1 — 2n)?E,/E 
F 1/(2 +- 2v + 3e) 


Here, /:, is the tangent modulus and e = (E£/EF,) — 1, 
where /, is the secant modulus. More directly, e = 


¢>/€, Where ¢, = € — €, and e, = a/F are the plastic 
and elastic strains at the buckling stress just before 
buckling. 

For Cases J and J, 8 = pe2/p; is equal to unity, by 


which, for these cases, in Table 1 


— ] + 3(1 + E/E 
“(2 4 Qv + Be) [1 + (1 — 2)E,/E] 


For Cases G and H, where the coordinate axes are 
in the planes of pure shear, values A, B, and F are given 
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by Eqs. (57) of reference 17 or Eqs. (67) of reference 


18. 
a 4(1 + e) 
; 4(1 — v2) + 4(2 — v)e + 3e? 
i i 2(2v + e) 
4(1 — v?) + 4(2 — vje + 3e? 
E,/E 


F =; —— 
3 — (1 — 2)E,/E 
In all these cases, values /, and e or E, in these equa- 
tions refer to a pure compressive stress equal to the 
equivalent buckling stress 


o,0, + 0,2 + 3rz,2)" (28) 


op = (,” — 
In Case B, for example, og = o,¢, and, in Case G, og = 
Tryer V 3. 

In Figs. 7-13, the coefficients 7 according to the for- 
mulas in Table | are given graphically as functions of 
this equivalent buckling stress og, for several ratios 8 = 
L/b for avional, a Swiss aluminum alloy, to which Figs. 
11-13 in reference 19 also apply. Similar graphs may 
be made for any other metal. Until these are available, 
a rough estimate may be made by an approximate de- 
rivation of the og, —-7 relation for a given metal from 
the given graphs for avional. This could be done by 
drawing first the -,/- and E,/E curves for the given 
metal, as sketched in Fig. 4 by the dashed curves. 
Then, for a given equivalent buckling stress (¢%)m 
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Fic. 9. Plastic reduction coefficients 4 for Case C 
functions of the equivalent buckling stress 
OB = Gx 


for several ratios 8 = L/b. 
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Fic. 11. Plastic reduction coefficients » for Case E as 
functions of the equivalent buckling stress 
CB = dz 


for several ratios 8B = L/b. 








Q0OO 0.25 050 0.75 1.0 r 


Fic. 10. Plastic reduction coefficients » for Case D 
functions of the equivalent buckling stress 
COB => Oz 


for several ratios B = L/b. 
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STABILITY OF 


of the metal, /,/# is read from Fig. 4. Subsequently, 
(on)a, belonging to the same ratio /,/E for avional, is 


is given by the dashed 


read from Fig. 8, where /,/ 
cuve. Then, /,/£ and n, belonging to (og), for avi- 
onal, are found from the /,/E curve in Fig. 7 and from 
the diagram for the pertinent case and for the chosen 
ratio 8 = L/b, respectively. The actual n value for 
the given metal is then found by assuming that, for 
on = (Cr) m, the ratio 

((E,/E) — n]/[(E,/E) — (E,/E)] 


ajo = (29) 


in Fig. 4 is the same as that found with (¢,), for avi- 


onal. 


Tue CriticaL Loaps P2 FOR CASE 2 
The critical loads P: for all cases of Table | may be 
expressed as 


P, = al(t + h)?/t]G, (30) 


where G, is the modulus of rigidity of the core. The 
pertinent coefficients a for various cases are given in 
Table 1. The author derived these coefficients in refer- 
ences 3 (Case A), 6 (Cases B—D and G-—I/), and 15 
(Cases A—-K). For Cases G, //, and A, the critical 
loads P2 from Eq. (30) refer to the shear forces P,,. 
THE CRITICAL Loaps P,, OF THE 
SANDWICH PLATES 


COMPUTATION OF 


A more direct expression for the critical load P,, is 
obtained by inserting Eqs. (18), (19), and (30) into 
Eq. (11), which yields 


Py = nkor*D, b? | 


where 
D, — D, fe gD, ' (31) 
4 le ) 
D, = E,h®/6(1 — v?) | : 
D, = E,h(t + h)?/2(1 — v?) 


g =a (nkor + a) 
(x2D,t)/[b°G,(t + h)?] 


r 


Since for Case A.1 values ko and mp for Po differ from 
values ky and m for P,; (Table 1), for this case Eq. (11) 
transforms to 

P= 7p,/e 
where 

D, = noD, + 0.91 moD, 

¢ = a/(0.91 mr + a) 

D,, D;, and r according to Eqs. (31) 
Values 7, ko, and @ in Eqs. (31) all refer to the same 
wave length Z and ratio 8 = L/b. If only buckling 
in the elastic range is considered, 7 = 1. 
P.. is calculated for some ratios 8 in order to find its 
However, in actual cases, 7 will prac- 


In this case, 


minimum value. 
tically always be found smaller than 0.1, in which case 
it is sufficiently accurate to calculate P,, for the ratio 


SA 
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8 which makes P;, and thus kj), a minimum. These 
values 6 and the pertinent ky) values may be found di- 
rectly from reference 14; for Cases G and //, they may 
be found from Fig. 6. For Case B, for example, 8 = 
l and ky) = 4, so that a = 2. For Cases G, H, and K, 
the critical loads P,, from Eqs. (31) refer to the critical 
shear forces P,, = 2hr;zy. 

With buckling in the plastic range, P,, has to be 
calculated for some ratios 8 = L/b. After finding a 
first approximation of the buckling load, as explained 
in the numerical example below, values n that belong 
to the stress o, from Eq. (28), which is equivalent to 
this approximated state of stress at buckling, are found 
from the pertinent graph. For Cases A-—F, for ex- 
ample, in Eq. (28), o, = P.,-/2h and o, = tr, = 0; 
for Case AK, rz, = P,,/2h = P,,/2h and, from Eq. (23), 
ytr,- After finding for which ratio 6 P,, is mini- 
mum, the assumed approximated load has to be ad- 
justed such that it checks with the thus computed more 


Ce = 


exact load. 


WRINKLING OF THE FACES 


Up to now, only buckling in relatively large half 
waves was considered, so that the strains of the core 
perpendicular to the middle plane could be neglected. 
This is, however, not allowed if buckling in waves of 
very small wave length, so-called wrinkling, is con- 


The experimental wrinkling stress’ was found 
3 


sidered. 
to be considerably lower than the theoretical one.*: 
Obviously, this is caused by the unreliable properties 
of the core, since in this case the critical stress depends 
on the local rigidity of the core, so that irregularities 
diminish directly the buckling stress. With buckling 
as a plate, with larger wave length, only the average 
rigidity plays a role. Therefore, in the elastic domain, 
the wrinkling stress o,,. found by Hoff and Mautner 
(reference 2, page 71) from tests should be adopted 


here, from which 
Porn = 2hoerw = h(E;EG,) (33) 


where £, is the elastic modulus of the core in a direction 
This formula 


perpendicular to the middle plane. 
may be applied to any state of buckling stress if o¢,, 
is considered as the principal compressive stress. 

In the plastic range, /, has to be replaced by nk. 
Since, here, practically cylindrical buckling occurs, 
for pure compression 7 is equal to its value for Case 
A.2 of Table 1, so that 


Pore = 2WNGerm = h(nEsEG,)"* = 
h(O.91AE,EG,)“ (34) 


With an arbitrary state of stress, wrinkling occurs in 
waves perpendicular to the maximum principal com- 
pressive stress, so that Eq. (34) applies for any state 
of stress if the X-axis is assumed in the direction of the 
Value A has to be cal- 


principal compressive stress. 
In the ratio 


culated from Eqs. (2la) of reference 19. 
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cae 2,24: 
; TABLE 2 4 for [ 
8 1 0.90 0.85 0.80 a Eq. 
(oer)a ———_—____—_—_————— 2,250 kg. per sq.cm. ———__—__—- a= 
n 0.87 0.85 0.84 0.83 - 
ko 4 4.04 4.104 4.20 the § 
a 2 1.81 1.7225 1.64 
7) 0.89 0.883 0.877 0.87 
D.(in.Ibs.) 23 , 660 23 , 500 23 ,350 3, 150 Shea 
(Per)r (Ibs. per in.) 1,300 1,275 1,271 1,275 A 
(ocr)r (Ibs. per sq.in.) 32,500 31,900 31,770 31,900 Al 
(Gcr)r (kg. per sq.cm.) 2,290 2,245 2 , 236 2,245 shou 
failu 
B = p2/pi in Eqs. (21a), p: and p: are here the principal P., = kor*D,/b*? = 1,476 lbs. per in. wre 
compressive stresses at the buckling load and p: < py. Fer = P,,/2h = 36,900 lbs. per sq.in. wae 
Furthermore, ,and e = (//E, — 1)inthese equations ee ; = ; " 
belong to the equivalent stress From Eq. (33), with £, = 330,000 Ibs. per sq.in., stret 
: siti Pow = h(E,EG,)'"* = 6,920 Ibs. per in. oo 
op = (p12 — pip2 + pr”) ”* ” (EEG) (0.4! 
: ‘ ; ane Hence, the governing critical load is 1,476 lbs. per i 
or, with reference to the coordinate axes given in Table — preneing cree Tend ts 1,676 Gen. per te 
1, too, from Eq. (28). For Case K, for example, Plastic Buckling APP 
pi,2 = (07/2) + [(0,?/4) + t2,7]” (35) If the critical stress o,, = 36,900 Ibs. per sq.in. is c 
For ‘Cases A-F in Table 1,'value » = 0.914 in Eq. above the proportional limit of the face material, the a 
(34) is given for avional by the curves for Case A.2 in actual buckling stress will be less than 36,900 Ibs. per al 
Fig. 7 or by the curves 6 = L/b = 0 for Cases B-Fin $4-™- but higher than the proportional (elastic) limit. ae 
Figs. 8-12, which are all similar. For Cases J and J, Assuming the faces to be made of avional, with an elas- nlat 
he a a ae ic limi ) . encom. = 98: — é 
it is given by the curve for these cases in Fig. 7. tic limit of 2,000 kg. per sq.cm. = 28,400 Ibs. per sq.in., “a 
P., has to be calculated from Eqs. (31) or (32), as 2” approximate value of the actual buckling stress ¢,, gs 
well as from Eq. (33) or (34). The smallest value P,,is ™#@Y be obtained as follows: From Eqs. (31), Per sas 
governing. o-, are approximately n(P.,)e or n(ocr),, where (o¢,)¢ is ; 
the elastic buckling stress of 36,900 Ibs. per sq.in. = sa 
on wep : neg 924 7 or 
NUMERICAL EXAMPLE 2600 kg. per sq.cm., so that o,, ‘= 2,600 kg. per sq. 
cm., where 7 belongs to og = o;,,; in Fig. 8. Assuming 
Given a metalite-type sandwich plate. The thick- also 8 = 1, from Fig. 8, on/n = ¢,,/n = 2,600 kg. per sq. 
ness of the aluminum faces is h = 0.02 in.; that of the om. if og = 2,250kg. per sq.cm. and 7 = 0.87, so that ap- T 
balsa core is t = 0.46 in. The plate is compressed in proximately o,, = 2,250 kg. per sq.cm. = 32,000 Ibs. pert 
longitudinal direction and is simply supported at its un- per sq.in. and P,, = (0.04) (32,000) Ibs. per in. = | whi 
loaded sides (Case B). It has a width b = 25 in. 1 280Ibs. per in. nell 
Find its critical load. ' ' 
More accurately, o,, is found by assuming this ap- 
Elastic Buckling proximate buckling stress (¢,,), = 2,250 kg. per sq.cm., 
; , . reading the accessory 7 from Fig. 8 for several 8 values ' 
In this case, 7 = 1. Furthermore, Poisson’s ratio . - 3 “ i = 
ata te ale is 0.3, E, = 10,500,000 Ibs. per sq.in and calculating the resulting P,, from Eqs. (31), de- wages 
- t , " yt , . - . “> ia . ¢ ‘ 9 . . . ™ . . " ig e 7 
and G, = 12,000 lbs. per sq.in. Hence, from Eqs. nanny (Perr ene See Ce SR turai 
; values are also given. Die 
(31), iet 
3 
ee ‘ se ae Hence, (P,,), is minimum for B = L/b = 0.85. To 
D, = E,h®/6(1 — v?) = 15.4 in. Ibs. (Perr ‘ B ‘ port 
: at fe oe , find the real value of P.,, which may be assumed also 
D, = E,h(t + h)?/2(1 — v?) = 26,600 in.Ibs. aie he : é Am; 
r = (x*Dd)/[b°G.(t + h)?] = 0.07 to be minimum for 8 = 0.85, the assumed and resulting ‘ 
en : , values, (¢¢,)¢ = 2,250 kg. per sq.cm. and (¢,,), = 2,236 Typ 
Since r < 0.1, it is sufficiently accurate (discrepancy kg. per sq.cm., should coincide. This appears to occur - 
about 1 per cent) to assume 6 = L/b equal to that for if (o.,)¢ = (oe), = 2,245 kg. per sq.cm. Hence, o,; = ee 
buckling of homogeneous plates. Since, from Table 1, 2,245 kg. per sq.cm. = 31,900 Ibs. per sq.in., and 194 
for Case B, P. = 2ho = 1,275 Ibs. per in. 6 
” , ; : . I-l 
ko = (1/8?) + 2+ 6? Obviously, the first found approximate value P,, = i oie 
me ; , — : Pe = 1,280 Jbs. per in. was sufficiently accurate in 7 
this is the ratio 6 for which ko becomes minimum, being nf er)e pe ae age Sart ' 
; this case. It underestimates P,, for 8 = 1, but, since Met 
8 = 1, from which kj) = 4. Hence, from Table 1, ; se wisi 194: 
value a = 1+ 6? = 2. This yields, from Eqs. (31) actually P,, for 8 < 1is smaller than for 8 = 1, it yields . 
a , : _ ’ a fairly accurate result. It would be worth-while to] ,... 
eg = a/(kor + a) = 0.878 check some cases in this respect, since this would highly 346 
D, = Dz + ¢eD, = 23,375 in.|lbs. simplify the calculation. ber, 
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Assuming P,, = 1,275 Ibs. per in. and og = og = 
9.245 kg. per sq.cm., 7 in Eq. (34) is found from Fig. 8 
for 8 = O, yielding » = 0.91A = 0.56. Hence, from 
Eq. (34), (Pew), = 5,700 Ibs. per in. Consequently, 
P.w Will be between 1,275 and 5,700 Ibs. per in., so that 


the governing critical load is P,, = 1,275 lbs. per in. 


Shear Strength 


Analogous to the design of lattice bars the plate 
should have a given shear strength in order to prevent 
failure by small irregularities. Requiring a shear 
strength of 3 per cent of P,,, this yields for the above 
cases of elastic or plastic buckling (0.03)(1,476) = 44 Ibs. 
per in. or (0.03)(1,275) = 38 Ibs. per in. With a shear 
strength for balsa wood of 7, = 175 lbs. per sq.in., the 
shear strength of the plate is® © S, = (¢ + h)r, = 
(0.48)(175) = 84 Ibs. per in. and is thus sufficient. 
OTHER 


THE SAME METHOD TO 


PROBLEMS 


APPLICATIONS OF 


The author’s method of split rigidities is more gen- 
eral than exposed here for the case of sandwich plates 
and may be applied to many other problems, as, for 
example, the buckling of latticed columns with batten 
plates or lacing,* '° rammed piles surrounded by soil,'° 
and torsional and flexural buckling of open sections.”° 
It was also used to establish simple formulas for the 
interaction of column and plate buckling in built-up 
columns of various cross sections, which will be re- 
ported in an N.A.C.A. paper. 
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SUMMARY 


By using the local Mach Number rather than the free-stream Mach Num- 
ber in the ordinary Prandtl-Glauert correction for subsonic two-dimensional 
isentropic gas flow, a simple new relation is obtained between the com 
pressible and the incompressible surface pressure coefficient The new 
relation includes thé effect of the specific heat ratio and gives a better agree- 
ment with experimental data than does the von Karman-Tsien relation. 


INTRODUCTION 


— PRANDTL-GLAUERT COMPRESSIBILITY CORRECTION for the 
local surface pressure coefficient is generally written (see 


reference 1, page 140) as 
Cy = Cn/V1 — Mat (1) 


Eq. (1) is derived from the linearized potential equation by neg- 
lecting all the square or higher order terms of the perturbation 
potential (see reference 1, page 125). The surface pressure co- 
efficient C, is then of the same order as the first power of the 
perturbation potential (see reference 1, page 126). 

In this form, Eq. (1) is rigorously applicable only to a very 
slender body having a cusped nose with no stagnation point. 
When applied to even relatively thin airfoils, it gives a reason- 
able prediction only for the portions near free-stream pressure 
and only for low values of the free-stream Mach Number. 

By an analytical continuation it is easily shown, however, 
that a better approximation for the local compressibility effect 
on any finite body is obtained by replacing the free-stream Mach 
Number 1/7. by the local Mach Number J/,; so that 


Cp = Cy/V1 — M1? (2) 


This essentially applies the Prandtl-Glauert linearized potential 
equation to the local flow field considered approximately uni- 
form. The same procedure, using the local Mach Number in 
the Prandtl-Glauert correction, was first suggested in reference 
2 and was later justified by the experimental results presented in 
reference 3 for the particular case of the subsonic flow field exist- 
ing behind a detached shock wave in transonic flow. 


THE NEw COMPRESSIBILITY CORRECTION 


The surface pressure coefficient may be written (see reference 


1, page 28) as 


Cp = ae 


—— 9) 
5 ol a ae (3) 
(1/2)po Unt yMat\pq 


where for isentropic flow (see reference 1, page 125) 


PL fl + [(y — 1)/2] Ma?) 7 — YD 
= : (4) 
bo 1+ [Cy —1/2]Mi? f 4 
Then substituting Eq. (4) into Eq. (3) and solving for the loca] 
Mach Number 1, we obtain the exact expression 


2 M «? (Y-D/y 2 
Yeemeey rer oe 
y-1 2 y¥-1 


Now, before substituting Eq. (5) into Eq. (2), we must be cog- 
nizant of the fact that only the first-order terms of C,, correspond- 
ing to the first order of the local perturbation potential, were re- 
tained in deriving Eq. (2) for the local flow field. Consequently 
Eq. (5) must be reduced to the first order in C, so that 


eee 
M,? = Mo? — (: +— Ua!) Mo? Coo (6) 


Then, substituting Eq. (6) into Eq. (2) and retaining only the 
first-order term of C, in the local flow field correction in the de- 


M,? = (11. + 


nominator, we obtain 


Cn 


;, M ~? Cs i dace bi :) 
V1l—-M.? / 9 ‘ cacti 
T Vi — Me? - 2 


Eq. (7) is a new relation for the two-dimensional compressible 
pressure coefficient and provides a better agreement with experi- 
mental data for airfoils than does the von Karman-Tsien rela- 


tion (see either reference 4 or reference 1, page 185) 


li+vVi1—-M,? < 


It is interesting to note that the term containing the specific 
heat ratio y is associated with the term 7.4, as would be ex- 
pected from a comparison of the von K4rman-Tsien straight- 
line pressure-density variation (corresponding to y = —1) with 
the exact isentropic variation (see reference 1, page 173). 
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Analogy 


Note on the Minimum Critical Reynolds 
Number and the Form Parameter* 


John C. Freeman, Jr. 
Department of Meteorology, The University of Chicago 
February 9, 1951 


» THE THEORY OF THE BOUNDARY LAYER of an incompressible 
fluid, the “form parameter” H = 6,/3,' where 6, is the displace- 
ment thickness and # is the momentum thickness, has been found 
useful as a parameter to characterize turbulent boundary layers’ 
and (empirically) as a parameter that determines R; (the mini- 
mum critical Reynolds Number based on the displacement thick- 
ness) for laminar boundary-layer profiles on a porous flat plate 


* This work was started at Brown University 
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In this note it will be shown 


with two different types of suction.? 
(m — 1) 


that, if sufficient suction of the type v(x) ~x 
on a flat plate with pressure gradient [such that the free-stream 
velocity u(x) = Uox™], Ri for many of the resulting profiles is 
The proof will be accomplished by 


>. . 
“is introduced 


essentially determined by H. 
showing that, following a suitable change of scale, all of the pro- 
files mentioned above can be closely approximated by members 
of the family of profiles computed by Schlichting and Bussmann* 
for the case for which the suction velocity u(x) ~ 1/Vx and for 
which there is no pressure gradient (i.e., m = 0). R, and H are 
both unchanged by the transformations involved, and thus the 
dependence of R,; on H is established. 


These results should be useful in the study of boundary-layer 
control along a curved surface. 

Suction is introduced on a flat plate with pressure gradient 
such that mwo(x) = Uox™ |with Hartree parameter B = 2m 
(m + 1)]so that the velocity into the plate is 


—c,VvUdm + 1) x™- 1/2 (1) 
If the stream function is assumed to be of the form 


¥(x, y) = V2vUox™ * l/(m + 1) f(s) (2) 


where 
= yV Uox™(m + 1)/2vx (3) 


then the Prandtl boundary-layer equation becomes the familiar 
Hartree equation 
f'' +Sf" = Bf? — 1) (4) 
with the new boundary conditions 
f'(0) =90, f(0) = G 
f'(o) =1 


This equation has been solved numerically by Hartree‘ for C, = 
(0 and several values of 6 and by Schlichting and Bussmann* 
with 8 = 0 and 1 for several values of C;. We are not interested 
here in exact solutions to Eq. (4) but rather in the approximate 


solutions obtained as follows: We set 


f+K=w (6) 
where 
K = B/f’"(0) (7) 
and neglect the right-hand side of the resulting equation, 
w'’’ ww'’ = By [w’’/w''(0)| + w’? — 1} (8) 


with boundary conditions 


w(0) = 0, 


w’'( o) 


(9) 


| 


This gives us 


w’’’ +ww’'’ = 0 (10) 


with boundary conditions above. 


This equation has been solved numerically by Schlichting and 
Bussmann? for several values of p from —*/. to +5. For values 


of p > 0, the quantity 


Bi [w’’/w’(0)] + w’? — 1} < w’’’/8 (11) 


for all s and thus can be neglected. This is sufficient to show that 
each one of all the classes of profiles obtained when 8 is varied 
between —0.198 and 2 is approximately the same as the c’ass for 
8 = 0, provided C; is such that p > 0 (i.e., the profiles depend 


essentially on p and not on £8). 


A similar but slightly less satisfactory demonstration that the 
profiles for arbitrary suction on a plate with pressure gradient 
depend only on pcan be given. 


x 
5 


Lin’s formula for the minimum value of R;! 
lows: 


is written as fol- 
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R, = 25w’’(0) (1 — w’)dz/|w'(2)]* (12) 
0 
where 2; is the value of z for which 
(13) 


a 2w'"(O)s | w'(s)w’’"(s) " 
rw (0) : -— 3 nee ; =z ().58 
w'(z) lw’’(s)]3 


Thus, it is obvious that R; is essentially independent of 6 and 


vo(x) and depends only on p. Since 
H = i,/0 = (1 — w’ )dz w'(1 — w’)dz (14) 
J 0 0 
for all of these profiles, it also depends essentially on p. Thus, 


the empirical results that H determines R, for any of the profiles 
with B = 0 (w(x) ~ 1/7 x|and p = 0 (see Fig. 18 of the Stability 
Monograph*) are given a theoretical basis and are seen to be 


only special cases of a more general result 
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The Influence of the Temperature Variation 
of the Specific Heats of Air in Shock-Wave 
Calculations 


J. Kestin 

Head, Department of Mechanical Engineering, Polish University 
College, London, England 

February 1, 1951 


T° A RECENT PAPER on the structure of shock waves, Meyer- 
hcff! discusses the influence of the temperature variation of 
the specific heats of air on shock-wave calculations. The method 
of calculation requires a great deal of numerical computation 
that can be avoided by an approximate, but sufficiently accu- 
rate, method. 

The shock transition is defined by the following three equa- 


tions: 
polly = pilly = Mm (1) 
po + muy = fi + mu, = mC (2) 
ho + (1/2)uo? = ky + (1/2)m,? = H (3) 


where h denotes enthalpy, H is the constant stagnation enthalpy, 
and the remaining symbols have the same meaning as in refer- 
These equations, together with the perfect gas equation 


p = RTp (4) 


ence 1. 


and the condition that for a perfect gas, the enthalpy is a function 
of temperature alone 

h = h(T) (5) 
with / given in tabulated form,? provide all the equations needed 
to calculate quantities with subscript 1 if those with subscript 
QO are given. Apart from the trivial solution 


Uo/U,y = Po pi = 7, T; =... 1 (6) 
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the set of equations leads to one solution for which the ratios in 
Eq. (6) differ from unity. 
If we analyze the dimensionless function 


oo 
u 
2= : = ok - + - 
* * * 
a pua a 
7 (7) 
R | T(h) 
Sayer -|+x 
a*?| x 
where x = u/a* and a* is the sonic stream velocity, determined 


solely by the value of H, we can easily prove that it has a mini- 
mum at 
(7a) 


“u=a* or x=! 


it being evident that it remains constant across a shock transi- 
tion. In Eq. (7) T(h) is the inverse of Eq. (5). Moreover, z is 
a function of a single variable, since, by Eq. (3), 


h = H — (1/2)u? _) (3a) 

= H — (1/2)a**x?\ se 

We shall now inquire into the conditions under which the very 
well known Prandtl’ equation 


Uo; = a*? or xox, = 1 (8) 


which is rigorous for constant specific heats applies to cases with 

variable specific heats. For this to be so, the relationship be- 

tween 7 and h must be such as to render the expression 

R {T[H — (1/2)a**x?] | 
7 ry r+ 2x 

a \ x f 


(7b) 
invariant with respect to the transformation of x to 1/x, in view 
of Eq. (8). 
we see that this condition can be fulfilled only if T is proportional 
to H — (1/2)a**x? (i.e., if T and hare linearly dependent). Thus, 
the validity of the Prandtl equation is limited to temperature 
ranges over which a linear relation between JT and h offers a good 
approximation to the tabulated function 4(7T). From Fig. 7, 
reference 1, we can see that this is the case over relatively large 
ranges of temperature, and, consequently, Eq. (8) can be expected 


Limiting our considerations to elementary functions, 
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to lead to sufficiently accurate results if the initial supersonic 
Mach Number and the stagnation enthalpy are not too large. 
The assumption of a linear relationship between enthalpy, }, 
and absolute temperature, 7’, is equivalent to the assumption of 
a constant value for cp, c,, and y, but by the use of Eq. (8) the 
cumbersome calculation of the appropriate mean value of y jg 
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Fic. C. Temperature ratio across normal shock for various 
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FIG 
natior 
2,000' 
heats, 


Frc 
varia! 
3,000 


avoid 
step 
Th 
spon 
wher« 
H, or 
of H 
the a 
locity 


from 
in ref 
4. 
Th 
of Ex 
le * 


obtai 


and . 





1personic 
arge, 
halpy, h, 
Iption of 
- (8) the 
> of ¥ is 


gin 
ae 


_ 


ous 


_= 


bs 





READERS’ 








Fic. D. Pressure ratio across normal shock for various stag 
nation temperatures. variable specific heats; (2) Ty = 
2,000°F. abs.; (3) Ts = 3,000°F. abs. constant specific 
heats, Y= 1.400. 
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2 3 + a 5s 
M. 
Fic. E. Mach Number behind normal shock wave 
variable specific heats; (2) Ty = 2,000°F. abs.; (3) Tu = 
3,000°F. abs. constant specific heats, y = 1.400. 


Furthermore, the approximation is used only in one 
in that involving the calculation of 1, from wp. 


avoided. 
step—i.e., 

The values of sonic stream temperatures, 7*, and the corre- 
sponding sonic stream velocities, a*, have been calculated else- 
where! for a number of stagnation temperatures and enthalpies 
H, on the basis of reference 2. Fig. A gives a plot of a* in terms 
of H and the corresponding stagnation temperature 7,,. With 
the aid of these values we can use Eq. (8) to calculate the ve- 
locity ratio 


Uy = u,/uo = a*?/u? = (a*/aq)?/ Mo? (9) 


from which all other quantities can be easily deduced, as shown 
in reference 1, particularly with the aid of the tables in reference 
4. 

The error can be assessed by calculating uo and uw, with the aid 
of Eq. (8) and by comparing the values of z: For example, for 
T,, = 3,000°F. abs. (H = 790.68 B.t.u. per Ib. °F.), 4, = 5, we 


obtain 


a* = 2,419 ft. per sec. 


uo = 5,748 ft. per sec. Zo = 2.46 
u, = 1,018 ft. per sec. a = 2.44 


and Az/z ~ 0.8 per cent, which is negligible 
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Figs. B, C, D, and E contain some results from which the in- 
fluence of the stagnation temperature can be assessed. The 
graphs also contain curves calculated with y = 1.400 for com- 


parison. 
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Convective Heat Transfer from Spheres in a 
Free-Molecule Flow 


F. M. Sauer 

Assistant Professor of Mechanical Engineering, University of 
California at Berkeley 

January 19, 1951 


lew PURPOSE OF THIS NOTE is to extend the heat-transfer 
calculations of Stalder and Jukoff! in the free-molecule flow 
régime of rarefied gas dynamics to include the case of heat trans 
fer from spheres. 

Using the conventional methods of kinetic theory, in reference 
1 is developed an analysis for the heat transfer from a plane sur 
face inclined at an angle to the flow velocity. If the surface of 
the sphere to be investigated is postulated to be isothermal, the 
results of reference 1 may be integrated over the sphere to yield 
the average heat-transfer rate. 

The convective heat flux from an elementary plane surface, 
dA, inclined at an angle @ with the flow velocity, is, for a mon 
atomic gas, given as* 


dq dA = 2ankTy — ak; (1) 


and for a diatomic gas (dumbbell type molecule) as 
dq/dA = 3ankTyp — af Ex + nkT;) (2) 
where 
E, = n|(mU?2/2) + gkT;| 
n = (NVm/2V/r)f{e-2? + Var x[1 + erf(x)]} 


* The symbols used are as far as possible those of reference 1 Eqs. (1) 


and (2) are the corrected form of equations presented in the reference 














Fic. 1. The coordinate system used in the analysis 
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Fic. 3. Convective heat transfer from spheres in a free-molecule 
flow. 
1 + 3/o\/m x[1 + erf(x) le?” 
¢ =1+ y, 
1+ Var x[l + erf(x) Je** 
x = (U/Vm) cos 6 


and from Fig. 1,dA = 2rr? sin @ dé. 
Substituting the above into Eq. (1) the heat flux (in dimen- 
sionless form) into a monatomic gas becomes 


1o/r dq Tp : 
= (4— — Met — 4 )e-** 
aNkTiVn (#) ( — ) r 


Tp 
(. - — 7M? -— 5) Vax|l + erf(x)| (3) 


The average convective heat transfer rate per unit area, assum- 
ing 7’, to be constant, is found by integrating Eq. (3) over the 
The result is 


8Q0/aNRT;U = g(s) (Tp/Ti) — f(s) 


sphere. 


(4) 


where 


l | 
g(z) = 4 E + - terfe (2) + 5— en] 


2 as 


, ' ‘ Res . as C“f(s) 
f(z) = (22? + 5)]} 1 + ~ terfe(z) | + (22? + 3) =r 


2? = (7/2) Mo? 
In the same manner the average heat-transfer rate per unit 
area into a diatomic gas is found using Eq. (2). This result is 


80/aNkRT:U = pz) (Tp/Ts) — oz) (5) 
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where 
_—_ l ; 
p(s) = 6) 14+ ~ terfc(s) + aes erf(z) 
- _ _. erf(z) 
g(z) = (22? + 7)} 1+ 2 terfc(z) | + (2s? + 5) - = 


The equilibrium temperature for the sphere neglecting radj. 
ation, Taw, may be found by setting the heat transfer equal to 


zero. For the monatomic gas, 

Taw/T; = f(z)/2(s) (6 
For the diatomic gas, 

Taw/T: = q(z)/p(z) (7 


The average convective heat-transfer coefficient for high-speed 
flow is defined as 

h = O/(Tp — Taw) (8 

Substituting Eqs. (6) and (8) into Eq. (4) and writing a dimen- 

sionless heat-transfer coefficient in terms of the Stanton Number 

gives, for the monatomic gas, 

Nu 


St = = 
Re-Pr 8y 


ay—1) , 
Q(z) (9 


Similarly, substitution of the corresponding equations for the 


diatomic gas yields 


St = [aly — 1)/8y]p(s) 10 


The results of the foregoing analysis are plotted in Figs. 2 and 


3. The ratio of the gas stream to stagnation temperature was 
calculated from the well-known equation 


T;s/T; = 1+ [(y — 1)/2]M? 
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Heat Transfer by Laminar Flow from a 
Rotating Plate 


Knox Millsaps and Karl Pohlhausen 
Office of Air Research, U.S.A.F. 
February 9, 1951 


—. EQUATION OF CONTINUITY, the Navier-Stokes equations, 
and the energy equation for the steady motion of an incom- 
pressible fluid when expressed in cylindrical coordinates (r, 3, 2 
and with velocity components (u“,, “9, uz) are, respectively, 


Ou, Uy 1O0uy | Ou: 
+ = () 
or r r ow Oz 
Ou, Uy OU, Uy? Ou, 1 Op 
Ur _ + uz =-— 
or r Ov r Os p or 
O7u, 1 Ou, Ur 1 O7u, 2 Ous O*u, 
v -* eS ee ae v 
or? r or ¥* r* Ov? r= Ov Os* 
ous uy ONg uyuy ous 1 Op 
Ur + + + Uz > = + 
or r Ov r Oz p rov 
O7u9 1 Ous uy 1 O'u9 2 Ou, O7u8 
a eS eau oa 
or* r or r* r* Ov? r? ov O2* 
Ou; uy OU, Ou; 1 Op 
Ur - + u, == 
or r ov 0: p Oz 
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oT oT _k 0? i 10T 
+ Uz _ + * + 
ov Oz Or? r or 
2T 2T y 
Lae. )+; + er? + 2es9*) 
r ov? Oz? 2 


The case of a rotating infinite plate (constant angular velocity 
w) heated at a constant temperature (f9) is characterized by the 


ug 


oT 
by | Ur = —- 
or r 


(rr? + e997 + 227 + 2¢r5? 


boundary conditions 
u(r, 3,0) = 0, uv, 0,0) = 0 
u(r, 3, ©) =0, ug(r, 0, ©) = O 


us(7, 0, O) = or 


The partial differential system may be separated and reduced? 
for laminar flow to ordinary nonlinear differential equations by 
assuming 

ur = rf(z), us = rg(s) 

uz, = h(z), pb = p(s) 

T = r*S(z) + g(z) + bo 


The reduced equations are, respectively, 


2f +h’ = ( 
P= te = we 
2fg + hg’ 

hh’ = —(1/p) p’ + vh" 
2c, fS + ¢,hS’ = (k/p)S" + of”? 
(k/p) (4S + q”) + v(4f? + 
where primes denote differentiation with respect to z. 


This system can be integrated numerically‘ or analogically. 
with a physical dis- 


” 


y” 
= ps 
44 


+g) 


crhq’ = 2h’?) 


A dimensionless representation, together 


cussion, will be given in a later paper. 
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Proper Use of the M.1.T. Tables for 
Supersonic Flow Past Inclined Cones 


Milton D. Van Dyke,* George B. W. Young, and Charles Siska 

The RAND Corporation, Santa Monica, Calif. 

February 1, 1951 

I APPLYING THE M.I.T. TABLEs'! to inclined cones, confusion can 
the results are 

than body axes. In the notation of the tables (which is used 

here throughout ), the velocity components, pressure, and density 


arise because tabulated in wind axes rather 


at any point are, to first order in e, 


u= t+ xecos ¢ (la) 
v =v + yecoso (1b) 
w= sesind (1c) 
Pp = p+ necos > (1d) 
p = pt tecos@ (le) 


The coefficients (#, 3, p, p, X, ¥, 5, n, €) are tabulated as functions 
of 6. Evaluating these expressions at 6 = 6, does not, in general, 
yield values of the flow quantities on the surface of the cone; it 
yields values at the position occupied by the cone before it was 


yawed (shown dotted in Fig. 1), which are of no interest. 


* Now at Ames Aeronautical Laboratory, Moffett Field, Calif. 
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Coordinate axes for inclined cone. 


Fic. 1. 


Values on the surface can be found only after the results have 
been transformed from wind to body axes. To first order in e, the 


flow quantities in body axes (subscript 7) are found to be 


ur =u +(x + v) ecosy (2a) 
_— Xu + v cot @) : r 
v = = — tujecosy (2b) 
WwW, = — yvese 6) e sin y (2c) 
a*({u + 3% cot 6) 
br =pttin pa = - e cos ¥ (2d) 
=F 
u + v cot 0 
Pr =pt+ pe - ). cos y (2e) 
a =e 


Here, the quantities (a, a, v, etc.) are to be evaluated from the 
tables using, for 6, the angle 6, measured from the body axis 
(Fig. 1). 


0, and (except for the trivial 
There- 


On the surface of the cone, y = 
case of v) Eqs. (2) reduce to Eqs. (1) with ¢ replaced by y. 
the conversion is unnecessary if values are required only on 
but this fact is not evi- 


fore, 
the surface and only to first order in e, 
dent until the conversion has been carried out. 

When terms in ¢«? are retained, the conversion must be made 
even for surface values and is rather complicated. 
the one quantity of greatest interest to aerodynamicists 
pressure on the surface—is found in terms of tabulated results as 


br/p = 1 + (n/p)e cos ¥ + [(po/p) + (yii?/2a*) + 


/s cot 6.(n p) |e? + \( po Pp) + (yi? 
'/s cot (n/p) |e cos 2¥ 


For example, 
the 


2a?) — 
(3) 


all quantities are evaluated from the tables at 6 = 4s. 
the expression is even more 


where 
Away from the surface, of course, 
complicated. 

As the cone rotates through an angle e, the shock wave remains 
a circular cone to first order, and rotates through an angle 6 
First-order values of flow quantities just behind the shock wave 
can be found from a slight modification of Eqs. (2). In the 
coefficients of «, each term but the first should be multiplied by 
a = 6/e; for example, at the shock wave, 


Ur = Ud + (x + ape cos ¥ (4) 


with (a, %, x) evaluated at 6 = 6,. The corresponding modifi- 
cation of terms in «2 would involve also the factors 8) and fy» of 
Report No. 5. 

It appears that the usefulness of the M.1I.T 
greatly increased by retabulating them in body coordinates. 
That is, the expression (x + 7) appearing in Eq. (2a) would be 
etc. The special values required for de- 


tables could be 


tabulated instead of x, 
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termining the flow just behind the shock wave could be included 
separately. 

The values of flow quantities calculated just behind the shock 
wave join smoothly to the remainder of the flow field only if ¢ is 
infinitesimal. Discrepancies arise if the results are applied as 
an approximation for finite values of « (as is done in starting a 
characteristic solution for an inclined ogival body). Prelim- 
inary study indicates that this defect can be eliminated by using 
relations such as Eq. (4) everywhere, with a, 8, and 6. varying 
through the flow field. 


a = 1 — [1 — (5/e)] (0 — 6.)/(@0 — 0) 


For example, to first order in e, let 


It should also be observed that the notation for drag coefficient 
is not uniform throughout the three volumes of tables. The 
symbol Ky has one meaning in Report No. 1, another meaning 
in Report No. 3 and most of Report No. 5, and a third meaning 
in Eq. (41) of Report No. 5. Interms of tabulated quantities, the 
conventional aerodynamic drag coefficient (referred to base area) 
is given to order e? by 


Cp = AKp + {Kpyn|AKp + (2/7 M*)] + (1/y M*)}e® (5) 


where A = 2 if Kp is taken from Report No. 1 and where A = 
8/m if Kp is taken from Report No. 3. 
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A Probe for Measuring Flow Inclination in a 
Supersonic Air Stream 


Robert D. Fusfeld 

> ao Laboratory, North American Aviation, Inc., Inglewood, 
Calif. 

February 5, 1951 


— USUAL DEVICE for measuring flow inclination in a super- 
sonic air stream makes use of a thin two-dimensional wedge 
or some modification of a wedge. Orifices placed on opposite 
sides of the wedge experience a pressure differential that serves as 
an indication of the inclination of the air stream. 
the Mach Number becomes large, the pressure coefficient differ- 
Other diffi- 
culties are that wedge-type devices are not especially easy to 


However, as 
ential due to the inclination angle becomes small. 


constsuct, they have a relatively large size, there is a tendency to 
deflect under air loads, etc. 

Other instruments sometimes used to measure flow inclination 
are free floating vanes, hot-wire probes, and blunt bodies with 
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FIG. 2 FROBE CAL/BRAT/ION 


pressure orifices near the leading edge. The last devices, similar 


to typical subsonic yaw heads, show most promise. The probe 
described here is a form of blunt body specially designed for use 
in supersonic flows. 

The probe is shown in Fig. 1. It consists of four small-diam- 
eter tubes soldered together with the leading edges cut off at some 
A conical nose shape is the simplest 


A probe 


large angle with the axis. 
to make and is one that gives satisfactory results. 
of this design can easily be fabricated in large numbers for use 
in rakes, etc. By making the instrument as small as possible, the 
disturbance introduced into the air stream can be minimized 
It is, therefore, suitable for use in close proximity to small models 
in wind-tunnel testing. Although designed specifically for use 
at supersonic speeds, the probe can also be used at subsonic and 
transonic speeds. 

Sample calibration curves for inclination angles in the plane of 
one pair of tubes are shown in Fig. 2. The pressure differentials 
at M = 2.87 are about two and one-half times as great as would 
be obtained on a thin wedge. At M = 1.86, the pressures are 


approximately the same as for a wedge. 


Formulas for Bending with Axial Tension 


Alfred S. Niles 


Professor of Aeronautical Engineering, Aeronautics Laboratory, 
Stanford University, Stanford, Calif. 


February 21, 1951 


4 pe GENERAL DIFFERENTIAL EQUATION for a beam-column with 
axial compression is 


(d?M/dx?) + (PM/EI) = $(w) (1 
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where M is the bending moment, P is the axial load, E is the 
modulus of elasticity, J is the moment of inertia of the section, 
and ¢(w) is a function of the transverse loading that is inde- 
pendent of P. P is considered positive, and EI/P may be re- 
placed by j*, also a positive quantity. The general solution of 


Eq. (1) is then 
M = C sin (x/j) + C. cos (x/j) + f(w) (2) 


When the axial load is tensile, in order to be able to continue to 
consider EI/P a positive quantity, Eq. (1) is usually replaced 
by 


(d?M/dx*?) — (PM/EI) = ¢(w) (3) 
which has the general solution 
M = T, sinh (x/j) + T>2 cosh (x/j) + g(w) (4) 


Formulas for C;, C2, and f(w) for various types of transverse 
load are presented in Table 14:3 of Airplane Structures, Vol. II, 
by Niles and Newell,’ and parallel formulas for what are here 
termed 7), 72, and g(w) are to be found in Table 14:5 of the same 
volume. Comparison of the tables shows that wherever W, w, 
or j appear in a formula for axial compression, they appear in the 
same place in the conjugate formula for axial tension. Also, 
whenever a trigonometric function appears in a compression 
formula, the corresponding hyperbolic function appears in the 
same place in the conjugate formula for axial tension. Some- 
times the signs in the conjugate formulas are the same. How- 
ever, sometimes they are reversed, and there is no obvious rule 
by which the signs in a formula for axial tension can be deter- 
mined from those in the conjugate formula for axial compression. 
The purpose of this note is to demonstrate a rule by which this 


can be accomplished. 


Eq. (1) would apply to the beam-column with axial tension if P 
were assumed to have a negative value. Then, if j? were taken as 
EI divided by the absolute value of P, it would be a_ pos- 
itive quantity with a real square root, but Eq. (1) would be- 


come 
d?M/dx? + [M/(ij)?] = o(w) (5) 

and its general solution would be 
M = K, sin (x/1j) + Ke cos (x/1j) + g(w) (6) 


where A,, As, and g(w) for any given type of transverse load 
would have the same forms as C, C:, and f(w), except that 7 
would be consistently replaced by 77. 

But, from the theory of complex numbers, 
cos (x/i7) = cosh (x/j) 
-itanh (x/j) (7) 


sin (x/ij) = —i sinh (x/j) 


tan (x/1j) = 
Substituting I°q. (7) into Eq. (6), the latter takes the from 
M = Ty, sinh (x/j) + T>2 cosh (x/j) + g(w (4) 
where 7; = —iK, and 72 = Ko. 

Thus we arrive at the following rules for developing the form- 
ulas for a given transverse load and axial tension from the conju- 
gate formula for axial compression: 

(a) Substitute 77 for 7 wherever the latter appears as a dimen- 
sion, as in Wj or wj?. 

(b) For each trigonometric function, substitute the correspond- 
ing hyperbolic function, multiplying the resulting hyperbolic 
sines and tangents by —i and the hyperbolic cosines and secants 
by + 1. 


For example, the formula for the bending moment on a pin- 
ended beam-column of length, Z, due to an axial compression, P, 
and a transverse concentrated load, W, acting at a distance, a, 
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from the left end, is, for x > a (see Case 4, Table 14:3, Airplane 
Structures, Vol. II, pages 92-93), 
Wj sin (a/j) sin (x/j) een : : 
M= : — Wj sin (a/j) cos (x/j) (8) 
tan (L/j) 
Making the specified substitutions, the conjugate formula for 
axial tension will be 
W1ij(—1z) sinh (a/j)(—7) sinh (x//) 
—i tanh (L/j) 
Wij(—i) sinh (a/j) cosh (x/j) 
Wj sinh (a/j) sinh (x/j) vay ; ; 
= . — Wj sinh (a/j) cosh (x/j) (9) 
tanh (L/7) 
Eq. (9) agrees with Case 4 of Table 14:5.! 
For a transverse load varying uniformly from zero at x = 0 tow 
at x = L, the formula for bending moment due to axial compres- 
sion is (Case 5, Table 14:3!) 


M = —[wj? sin (x/j)/sin (L/j)] + (wj?x/L) (10) 


Making the specified substitutions, the formula for axial tension 


1s 
ul w(7j)?(—7) sinh (x/j) wij) *x 
: —i sinh (L/j) s 
= [wj* sinh (x/j)/sinh (L/j)] — (wj?x/L) (11) 


For a uniformly distributed load w and end moments J, and 
M2, the formula for axial compression is (Case 3, Table 14:3!) 


Mz — wj? — (M, — wj?) cos (L/7) 


M= sin (x/j) + 


sin (1/7) 
(M, — wy?) cos (x/j) + wy? (12) 
Making the specified substitutions, the formula for axial tension 
is 
Me — w(ij)? — (M, — wi?j?) cosh (L/7) 


M = ar . 
—i sinh (L/j) 


(—7) sinh (x/j7) + (Mi, 


— wi*j?) cosh (x/j) + wi?7? 
M2 + wj? — (M, + wj?) cosh (L/j) . : 
= : . > sinh (x/7) + 
sinh (L/j) 
(M, + wy?) cosh (x/j) — wj? (13) 


It may be noted that in the first example the signs of the terms 
in the two formulas were the same; in the second, they were re 
versed: and in the third, some were reversed while others re 
mained unchanged. 

All of the expressions of Table 14: 5,' including the terms for the 
Three Moment Equation, can be obtained in this manner from 
those of Table 14:3. 
T; in this note), the results must be multiplied by —/ to take 
account of the substitution of —7 sinh (x/j) for the sin (x/j) by 
which the C, of Table 14:3! would be multiplied in a formula for 


In obtaining the expressions for C; (called 


bending moment. 


The formulas for the total shear on a beam-column can be 
treated in the same manner as those for bending moment. The 
formulas for slope and deflection under axial tension can be de- 
veloped by using a negative value for P in those for axial com 
pression, and vice versa. 

In making this study, a typographical error was found in Case 
10 of Table 14:5! (axial tension and a clockwise couple applied in 
the span). The constant C, for x < a should be preceded by a 
minus sign. 
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Simplified Compressible Laminar Boundary- 
Layer Theory* 


Arthur N. Tifford 

Department of Aeronautical Engineering, The Ohio State 
University, Columbus 

January 29, 1951 


2 ancsnig RECENT YEARS several excellent mathematical anal- 
yses of the compressible laminar boundary layer have 
been published. Unfortunately, however, the essential unity of 
the analytical results at Mach Numbers as high as 5 (under some 
circumstances) has apparently not been noted by others. It is 
the purpose of this communication to emphasize that unity. 

The writer was first concerned with the effect of variable fluid 
properties upon the surface characteristics of the boundary layer 
when doing the research reported in reference 1. Busemann? and 
Eckert and Drewitz* had already shown that, in the case of a 
compressible gas flowing along a flat plate, the effect is very small 
when the Mach Number of the flow is not too large—e.g., 3 per 
cent at Mach Number = 2. The writer concluded on the basis 
of qualitative reasoning that it must also be small in the general 
case of a surface of arbitrary shape (except in the neighborhood 
of a shock wave) if the ‘‘local’” Mach Number of the external flow 
is somewhat lower. 

More recently, Howarth‘ has studied the effect of compressi- 
bility upon the boundary layer at Mach Numbers near unity 
when the local free-stream velocity increases or decreases linearly 
with the distance along the surface. Let us reconsider his re- 
sults. 

In the case of the linearly increasing local free-stream velocity, 
U = Bx, he found the local skin friction coefficient, 7/povo' 28/2 


* Work supported by the Office of Air Research, U.S.A.F. 
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or 7/p)U? (vo/Ux)'/*, to decrease from 1.23 to 0.87 as the local 
Mach Number increases from 0 to Bx/ayp = 1—i.e., local Mach 
Number based on free-stream static temperature increases from 
0 to 1.12. The quantity, v'/28-*/%x-1(0u/0Y)o, on the Other 
hand, increases only 5 per cent, from 1.23 to 1.29, over the same 
range of Mach Numbers. However, the latter quantity is “‘idep. 
tically”’ the local skin friction coefficient based upon the loca| 
fluid properties at the wall, r/pyU?(vw/Ux)' 2 The expectation 
of the writer, that regardless of pressure gradient only very small 
surface compressibility effects (properly defined) exist in this 
Mach Number range, is thus supported. f 

Let us turn now to the case of the linearly decreasing loca] 
free-stream velocity,t U = Uy — Uix. Howarth has found that, 
at the point at which 8(U,;x/U») = 0.1, the friction coefficient, 


Arx'/?/povo'/? Uo'/*? or 417/po Uo? (vo/Uox)'/2, is 1.228 when M, = 
1 as compared with 1.226 in the incompressikle flow case—excel- 


lent agreement. When we define the local friction coefficient jp 
terms of the local free-stream velocity and the local fluid proper 
ties at the wall—as before—excellent agreement (i.e., better than 
1 per cent) with the incompressible flow case is still obtained. 

It does not seem out of place to remark here that in a Readers’ 
Forum note,> Chapman and Rubesin’s results for the noniso 
thermal flat plate® have been shown to be essentially the same as 
those obtainable on the basis of fluid properties being ‘‘effectively’ 
constant at the local wall value. The argument is probably more 
easily followed when the local surface shearing stress and the local 
rate of heat transfer are examined, rather than their average 


+ It is interesting to replot the boundary-layer velocity profiles of Fig. 5 
reference 4, as in Fig. A The merging of the profiles associated with dif 
ferent local Mach Numbers is due solely to the use of the local fluid kine 
matic viscosity at the wall, instead of its stagnation value, in the nondimen 
sional boundary-layer thickness coordinate 

} The subscript zero refers in this case to conditions in the free stream at 
the leading edge rather than to conditions at a stagnation point. 
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value over the plate. For example, the local surface shearing 


stress has been found in reference 6 to be given by Eq. (57a). 

tw = [f"(0)/4]UeCwV po Vana /Cx (57a) 
At low free-stream temperatures, Cy and C are identical. Eq. 
(57a) thus reduces to the solution obtained on the basis of the 


fluid properties being constant at the local wall value. 


Tw = [f"(0) 1)V pw U ony x 


At high free-stream temperatures, however, the definitions of 
C» and C allow them to be different in general, the former being 
appropriate to the local surface fluid physical properties; the 
latter, to the average surface fluid physical properties. Consider, 
though, the strong asymptotic nature of the boundary layer upon 
which essentially all boundary-layer analysis leans. Increasingly 
crude assumptions may be made at points lying successively 
further ahead of the point at which the boundary-layer character- 
e.g., the conditions at the leading edge of a 
Conversely, increasingly pre- 


istics are sought 
flat plate are completely ignored. 
cise assumptions must be made at points lying successively closer 
to the point of interest. It is to be expected therefore that the 
best value for C is essentially the local Cy, rather than an average 
Cw. On this basis the results of the two analyses become identi- 
eal “under all conditions.” 

It is not the intent of this note to de-emphasize the importance 


of compressible laminar boundary-layer analyses in which fluid 
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properties are allowed to vary. They are needed for a fuller 
understanding of the details of the flow, for the evaluation of sur- 
face phenomena at very high Mach Numbers and at certain com- 
binations of high local Mach Numbers and surface pressure gra- 
dient, etc. It is rather the intent to emphasize that, in many 
current applications, the characteristics of the laminar boundary 
layer at any point on a surface may be determined within engi- 
neering precision on the basis of the simpler analyses extant in 
the literature, assuming fluid properties throughout the boundary 
layer to be the same as at that particular point. 


REFERENCES 


| Tifford, A. N., The Thermodynamics of the Laminar Boundary Layer of 
a Heated Body in a High-Speed Gas Flow Field, Journal of the Aeronautical 
Sciences, Vol. 12, No. 2, pp. 241-251, April, 1945 

2 Busemann, A., Gas-stromung mit laminaren Grenzschicht entlang einer 
Platte, ZAMM, Vol. 15, No. 23, pp. 22-25, 1935 

3 Eckert, E., and Drewitz, O., The Heat Transfer to a Plate in Flow at High 
Speed, Forschung a.d. Gebiete d. Ingenieurwesens, Vol. II, No. 3, pp. 116 
122, May-June, 1940; also N.A.C.A. T.M. No. 1045, May, 1943. 
, Concerning the Effect of Compressibility on Laminar Bound- 
of the Royal Society of London, 


4 Howarth, L 
ary Layers and Their Separation, Proc 
Series A, Vol. 194, No. 1036, pp. 16-42, July, 1948. 

5 Tifford, A. N., On the Surface Effects of a Compressible Laminar Boundary 
Layer, Readers’ Forum, Journal of the Aeronautical Sciences, Vol. 17, No 
3, pp. 187-188, March, 1950 

6 Chapman, D. R., and Rubesin, M 
files in the Compressible Laminar Boundary Layer with Arbitrary Distribution 
16, No. 9, 


W., Temperature and Velocity Pro 


of Surface Temperature, Journal of the Aeronautical Sciences, Vol 
pp. 547-565, September, 1949 








Notices should be sent directly to: 





Changes of Address 


Since the Post Office Department does not as a rule forward magazines to forwarding addresses, 
it is important that the Institute be notified of changes in address 30 days in advance of publishing 


date to ensure receipt of every issue of the JOURNAL and REVIEW. 


Institute of the Aeronautical Sciences, Inc. 
2 East 64th Street, New York 21, New York 




















INSTITUTE of tke AERONAUTICAL SCIENCES 


1932— 


1951 


Since its incorporation in October, 1932, as a scientific membership society to advance the art and science 
of aeronautics, the Institute of the Aeronautical Sciences, in over eighteen years of service to aviation, has 
become the representative technical society for the aviation industry and profession. It is the only organiza- 
tion in the United States which brings together all of the sciences and branches of engineering and tech- 


nology which are applied to. aeronautics. 


Publications 


Papers presented at meetings or submitted to the Edi- 
torial Board are published in the JourNat or THE AERO- 
NAUTICAL ScigNcEs or the ABRONAUTICAL ENGINEERING 
Review. The Journat prints in each issue several full- 
length scientific papers on new research and developments 
in various fields applied to aviation. The Revisew keeps 
members and subscribers up to date on aeronautical news 
and literature through reviews of new books and periodical 
articles, government publications and trade literature and 
sear papers on applied engineering and aircraft pro- 

uction. 

The AgronauticaL ENGINEERING CatTatoc, published 
yearly, is a guide to sources and specifications of materials, 
parts, and accessories used in the design and production of 
aircraft, aircraft engines, and parts. 


Research Facilities 


Through the Aeronautical Archives of the Institute, 
members have for their use the most complete aeronautical 
library and reference research facilities maintained by any 
technical society. At the New York building of the Insti- 
tute, The W. A. M. Burden Library of books, reports and 
periodicals, the Aeronautical Index of subject files, bio- 
—— files and bibliography, and the historical and art 
collections in the Archives are available for study. The 
Pacific Aeronautical Library at 7660 Beverly Boulevard, 
Los Angeles, maintains a reference collection and reading 
room for members and loans books to company libraries 
in the area. Through The Paul Kollsman Lending Li- 
brary, aeronautical books are lent by mail without charge 
to members anywhere in the continental United States. 
In addition to a complete collection of standard aero- 
nautical reference works, the Kollsman Lending Library 
has available for loan new aviation books as they are re- 
viewed in the AERONAUTICAL ENGINEERING Review. 


Membership 


In addition to membership in professional societies 
representing their particular fields, specialists applying 
their knowledge and experience to aeronautics find mem- 
bership in the Institute of great value in the contacts made 
with others engaged in aeronautical work and in the broad 
coverage of all aeronautical problems furnished through its 
meetings and publications. All applications for member- 
ship are carefully reviewed by the a Committee, 
which recommends the grade of membership for which the 


applicant is found eligible. There are advanced grades of 
membership—Fellow and Honorary Fellow—to which 
members may be elected as their experience warrants. 
Thus, a member receives due recognition—more valuable 
because it is given by his colleagues—for his contributions 
to the progress of aeronautics. New applicants may be 
elected by the Membership Committee to the grade of 
Associate Fellow if they have been engaged for not less 
than ten years in the practice of the aeronautical sciences 
and shall have been in responsible charge of important 
scientific or engineering work, or shall have otherwise 
made outstanding contributions to the aeronautical sci- 
ences; to the grade of MEMBER if they have been en- 
gaged in professional work for at least eight years and have 
acquired a recognized standing in engineering, design, or 
other special work applied to aeronautics; to the grade of 
Associate Member if they have acquired a recognized 
standing in an administrative capacity in the aviation in- 
dustry; to the grade of Technical Member if they are recent 
engineering school graduates or are engaged in technical 
aeronautical work. Student Members are admitted by 
application through Student Branches organized at their 
schools. Acronautical companies, as well as individuals, 
are affiliated with, and participate in, the support of the 
Institute through Corporate Membership. 


Dues 


An entrance fee of $10 is required of all new members 
who apply for admission to a grade of membership to be 
specified by the Membership Committee (Corporate Mem- 
bers, Foreign Members, and Student Members being ex- 
cepted). 

Annual membership dues include a subscription to the 
AERONAUTICAL ENGINEERING Review. Dues for the vari- 
ous grades of membership are as follows: Fellows— 
$16.50; Associate Fellows—$16.50; MEMBERS—$13.50; 
Associate Members—$13.50; Technical Members—(26 
years of age and over) $13.50, (under 26) $6.50. 

The regular subscription price for the JouRNAL OF THE 
AERONAUTICAL SciENcEs and the AERONAUTICAL ENGINEER- 
ING Review is $12 and $3.00 per year, respectively. Mem- 
bers may subscribe to the Journat at $4.00 per year. 

Those who are engaged in aeronautical work or who 
have an interest in any technical phase of aviation will find 
the services of the Institute to be a necessary aid in keeping 
abreast of developments in these times of such rapid prog- 
ress. An application form and further information about 
membership can be obtained from an Institute member or 
by writing to the Secretary. 


INSTITUTE OF THE AERONAUTICAL SCIENCES, INC. 


2 EAST 64TH STREET 


NEW YORK 21, N.Y. 














